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ON POSSIBLE REALIZATIONS OF MACH’S PROGRAM! 
F. A. KAEMPFFER 


ABSTRACT 


By ‘‘Mach’s Program” is meant the intention to understand all inertial 
effects as being caused by gravitational interaction. Within the framework of a 
simplified theory of gravitation, based upon ideas by Hund and Sciama, it is easy 
to invent primitive cosmological models for which Mach’s Program is realized. 
In particular, it is possible to understand the inertia of a linearly accelerated body 
as being due to radiative gravitational interaction with all other masses in the 
universe. If it is true that anti-matter has negative gravitational mass, as has 
been suggested by Morrison and Gold, then one may look upon Newton’s second 
law as a manifestation of the preponderance of matter over anti-matter in our 
universe. 


1. INTRODUCTION 


Newton’s second law is certainly one of the most obscure of all the un- 
derstandable relations underlying our description of the physical world in 
which we find ourselves. Anyone who has ever tried to explain this law to a 
person who insisted on asking questions will know the difficulty of giving good 
reasons for the two faets embodied in it, namely, that production of an accelera- 
tion, rather than some other derivative of the position, requires a force exerted 
on the body to be accelerated, and that this so-called law is true only if the 
acceleration to be produced takes place with respect to certain mysterious 
coordinate systems, called ‘inertial frames’. 

It is curious that the existence of inertial frames is usually demonstrated for 
the kinematically more complicated case of central acceleration, like in 
Newton’s famous rotating bucket experiment, instead of for the simple case 
of linear acceleration. 

For later reference let us state here an experimental definition of what is 
meant by an inertial frame in the case of rotational motion. Suppose two 
observers O and O’ are fixed to their respective coordinate systems S and 5S’, 
which are in rotation with respect to each other around a common axis (see 
Fig. 1). Then the motion of O’ in S can be described by angular velocity wo, 
position r, and velocity Vv = (w Xr). Elementary kinematics tells us then that 


1Manuscript received September 20, 1957. 
Contribution from the Department of Physics, The University of British Columbia, 
Vancouver 8, B.C. 
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the acceleration of O’ in S consists of two terms, Coriolis acceleration c = 


2(#Xv) and centrifugal acceleration f = —(@X(@Xr)), which add up to 
the central acceleration 

(1) a = 2(wXv)—(#X(oXr)) = (woX(eXr)). 

Similarly, O has in S’ angular velocity w’ = —a, position r’, velocity v’ = 
(w’ Xr’) = —v, and therefore a central acceleration 

(1’) a = (ao Xa xP). 


The remarkable experimental fact is that there exist observers who can 
maintain central accelerations, as described above, without requiring any 
forces to do so. The coordinate system tied to such an observer will then be 
called an inertial frame. Only if the rotational motion of an observer takes 
place with respect to an inertial frame will the production of the corresponding 
central acceleration require a force F according to Newton’s second law F = 
ma. 





Fic. 1. 


The only clue we have regarding the reason for the existence of inertial 
frames is the astronomical observation that all inertial frames coincide with 
those frames in which the observable universe as a whole appears to be in a 
state of unaccelerated motion. Newton himself apparently did not recognize 
this coincidence as a link, and it remains a psychological puzzle why so few 
physicists have been attracted by the almost obvious suggestion made by 
Berkeley and, more emphatically, by Mach, who stated clearly that the dis- 
tribution of matter in the universe ought to be considered as the cause of the 
inertia of all bodies embedded in it. 

Imagine that throughout human history our earth, by some freak of climate, 
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had been covered with thick clouds, so that yesterday we were without astro- 
nomical experience but in possession of the laws of mechanics, including know- 
ledge of the disturbing and quite mysterious existence of inertial frames, which 
are not tied to any obviously simple state of motion on this earth. Imagine 
further that last night another freak of climate had dissipated for a few hours 
all clouds and for the first time we were seeing the stars around us. A few simple 
measurements, quickly done, would have revealed that the system of stars 
and the inertial frames coincide. It is hard not to believe that even the most 
soberminded of witnesses would have agreed that here might be a possibility 
to at last understand what so far had to be accepted as one of those facts. 

Faced thus with the suggestion of a causal relationship between the universal 
distribution of matter and the inertia of a body immersed in it, one is led to 
search for an interaction that may be made responsible for this conjectured 
state of affairs. Since gravitation is apparently the only interaction that indis- 
criminately affects all bodies having inertia, regardless of whether the bodies 
have or have not additional properties like mesonic charge, electric charge, 
Fermi charge, etc., it is extremely tempting to adopt and to attempt a realiza- 
tion of Mach’s Program. Adopting Mach’s Program means to declare one’s 
intention to understand all inertial effects in terms of the gravitational 
interaction between the masses in the universe and the body that exhibits 
inertia. Realizing Mach’s Program means to derive Newton's second law from 
a suitable theory of gravitation. 


2. THEORY OF GRAVITATION 

It has been known for a long time that the very existence of Coriolis and 
centrifugal forces cannot be understood from the point of view of Mach’s 
Program within the framework of a scalar field theory of gravitation, and it was 
precisely this insufficiency of Newton's theory of gravitation that was one of 
the motivations for Einstein’s attempt to describe gravitation in terms of a 
tensor field. Unfortunately, it turned out that the various versions of Einstein's 
theory do not by themselves constitute realizations of Mach’s Program, and it 
seems that to date not even any special cosmological model has been found, 
formulated in terms of general relativity theory, that fully realizes this program 
(see Bondi 1952). 

Therefore, for the present discussion, let us cast aside all those cosmological 
models that require Riemann geometry, and let us simplify the description of 
the probably rather complicated phenomenon of gravitation, by devising a 
vector field theory analogous to electromagnetic theory, which combines ideas 
due to Hund (1948) and Sciama (1953, 1957) with a formal procedure by 
Gupta (1954), and which enables one to demonstrate how Mach's Program may 
be realized in principle within primitive Euclidean models of the universe. 

In analogy to the electromagnetic field, described by a polar electric vector 
E and an axial magnetic vector H, which act upon any charge e with the force 


F = e[E+(1/c)(vXH)], 


it will now be assumed that it is possible to describe the gravitational field in 
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terms of a polar vector f and an axial vector g, which act upon any gravitational 
mass m with the force 


(2) F = m(f+(1/c)(vXg)]. 


The fields f and g are supposed to be produced by gravitating masses and 
currents, analogous to the production of the fields E and H by electric charges 
and currents, but with the important additional requirement that all forms of 
energy, including the field energy of f and g itself, are considered as sources of 
gravitation. This means that those terms in the field equations which corre- 
spond to the source terms in Maxwell's equations are from the outset to be 
considered as functions of the fields f and g themselves, and it becomes neces- 
sary to use a description by successive approximations, similar to the method 
employed by Gupta (1954). 

It is clear that in zero approximation, i.e. if one neglects both the effects of 
g and the gravitational effect of the field energy itself, the description must 
revert to Newton’s theory of gravitation, where there is just a polar vector 
f‘” which satisfies for given gravitational mass density p the equation 


(3) div f = —4ryp (y = gravitational constant) 


and which acts upon any gravitational mass m according to Newton's law of 
gravitation with the force 


io") Fyewton = mf, 


If in first approximation one does take into account the effect of g, but 
considers only fields that are sufficiently weak, so that the gravitational effect 
of the field energy itself may still be neglected, one will have a description in 
terms of f‘) and g“’, which is the complete analogue to Maxwell's equations, 
namely 


a) div f” = —4rvyp, curl g = a ypv +--—--. 


The equations enabling one to consider f and g as derivable from four- 
potentials, 


(4) div g = 0, curl f = —- =, 


can be required to hold in all approximations. One may thus introduce, as is 
usual in electromagnetic theory, potentials by 
1 0A 


(5) ¢ = curlA, f = —grad as 


It turns out that this first approximation is already sufficient for the deriva- 
tion of Newton's second law in the case of a body accelerated linearly, and for 
the understanding of Coriolis force in the case of a body accelerated centrally 
with respect to the universe. However, in order to understand the centrifugal 
force experienced by a body rotating with respect to the universe one has to go 
to the next approximation, in which the gravitational contribution from the 
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energy of the fields in first approximation, namely u“) = —(1/8my)(f?+ 
¢@‘?), is added to the gravitational mass density p by putting p" = p+ 
(1/c?)u™, so that, for example, the equation governing the production of f by 
gravitating sources reads in second approximation 


(3) div £2 = —4ry[o+(1/c)u], 


with a similar equation, containing the gravitational contribution from the 
energy flux density vector s“) = —(c/4ry)(f" Xg™), namely (1/c?)s™, 
added to the current in the equation governing the production of g. 

One may thus proceed by successive approximations to a more and more 
accurate formal description of gravitational phenomena. This procedure can 
be written down more concisely in terms of the Lagrangian formalism. One 
starts in first approximation with a Lagrangian density 


a. dA, OA, 
2Y Ox, OX; 


+ (0) 


(6?) L” = — 4rA,j; 


with a given, field independent, current four-vector j°°, so that the corre- 
sponding field equations are 


» 6L‘)/6A, = 0, ie. OA, = 4eryji”. 


~J 


( 


The equations in second approximation are obtained by adding to the source 
j\° the gravitational contribution from the energy momentum density, namely 
(1/c?)T“”, so that 


v0? 
(7°) OA, = 4ayg = 4ay[j0+(1/2)TO)] 
where 
a) » ee dA, a + Lt 


dx9 0(0A,/8x,) 


These equations (7°) are in turn derivable from a Lagrangian density L@.so 
1 gran 

that they become identical with 6L°)/64, = 0. In general, the equations in 

nth approximation are identical with 


(7) BLO) 8A, = (Aey/2)TE™. 


It is an interesting and open question whether one can derive expressions for 
the field equations in the limit »—© in closed form. The pursuit of this 
problem is outside the scope of the present discussion. 


3. NEWTON'S SECOND LAW FOR LINEAR ACCELERATION 


Turning our attention now to the equations (7?) describing the field in first 
approximation, let us consider the case of a body of gravitational mass m 
accelerated linearly with respect to the other masses constituting the universe. 
From the point of view of an observer tied to m, a mass element pd V at distance 
r will be moving with velocity v and linear acceleration V (see Fig. 2). According 
to the field equations this element pd V will produce at distance r a gravitational 
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pdV 





3 total 


Fic. 2. 


field, analogous to the static plus radiation field in electrodynamics, which, 
taking only terms of lowest order in v/c, has the form 


|r “a rt hiex (xv) fod 
(9) 
qd) _ 1 (1) 
ig = vif Xr). 


This field will then, according to equation (2), exert a force, again taking 
consistently only terms of lowest order in v/c, 


(10) F = mf), 


Now let us replace the actual universe by a primitive model, consisting of a 
sphere of radius R and constant mass density p, with m being located at the 
center. Then the total force exerted on m will be obtained by the integration 


(11) | ee m f£dV 


which is elementary and yields the interesting result that the total force is 
parallel and proportional to the acceleration V, 


(12) Frota = (4r-ypR?/3c?)mv 


with all other components averaging to zero. It is clear that this force will be 
experienced by m as a drag, proportional to V, and will have to be compensated 
by a force F’ equal and opposite to it if the state of acceleration WW = —v is 
to be maintained. Thus Newton’s second law in the case of linear acceleration 
can be understood as being due to the radiative gravitational interaction of 
the universe with the mass to be accelerated. 

In order to obtain numerical agreement with experience one requires in the 
approximation considered here the relation, obtained by inspection from (12), 


(13) 4rypR?/3c? = 1, 
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which is familiar from cosmology (see Bondi 1952) and which seems in reason- 
able order-of-magnitude agreement with the observational data available for 
p and R in this model. 


4. NEWTON’S SECOND LAW FOR ROTATIONAL MOTION 


Consider next a body of gravitational mass m moving with respect to the 
inertial frame on a circular path of radius r. Seen from an observer tied to m 
the universe will be revolving around an axis going through r = 0 with angular 
velocity w. The inertial forces experienced by m are then Coriolis force 
—2m(vXw) and centrifugal force —m(wX(@Xr)). If the fields f and g pro- 
duced by the rotating universe are considered as the cause of these forces 
according to equation (2), they must have the form 


(14a) f = oF, 
(14)) & = —2cw. 


The field equations in first order (3) can be solved easily in this case, if 
again a model universe of radius R and constant density p is taken. By sym- 
metry f“ will vanish at the center, but for 8“ one obtains, in complete analogy 
with the magnetic field produced at the center of a rotating homogeneously 
charged sphere, 

(15) g°) = —(2rypR?/3c*) 2cw. 
By comparison with (140) it is thus found that in order to yield the observed 
Coriolis force the model under consideration must be subject to the relation 


(16) 2rypk?/3c? = 1, 


which up to a factor 2 is identical with the order-of-magnitude relation (13) 
obtained in the preceding section from quite different considerations. 
It was Hund (1948) who first observed that the very large negative energy 


density corresponding to §""’, 
(17) u") = —(1/8ry)g)? = —we?/2ry, 


is sufficient to account for the centrifugal force by virtue of the equations in 
second order (3°). Indeed, wu“ corresponds to a negative mass density 


(18) p) = —w*/2ry 

which by gravitational repulsion, as demanded by the field equation 
(19) div £2 = —4ryp™ = 2’, 

gives rise to the observed centrifugal force 

(20) f® = wr. 


(Note that, because of cylindrical symmetry, div r = 2.) 
It should be stressed that from the point of view of this model the inertial 
forces are caused mainly by the distant masses, since the solutions of the field 
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equations always contain an integral [*rdr. It is instructive to compute the 
¢-field produced by the rotation of a single galaxy near its center. Let R’, 
p’, and w’ be radius, mean density, and angular velocity with respect to the 
inertial frame of a typical galaxy. Then one should expect near its center an 
additional g-field of the order 

(21) g’ = yp'R”°ow'/c 


corresponding to a negative mass density 
(22) pS —yp?Ra!2/c4. 


Using the reasonable values p’ ~ 107?! g./cm.’, R’ ~ 10” cm., w’ ~ 107" sec.—! 
one finds 
(23) —p = 10-* g./cm.?, 


which appears to be quite beyond observation at the present time. In fact, 
such a negative mass density would cause gravitationally an effect equal to the 
inertial forces observed on a rotating system revolving with angular velocity 
Q, where 

(24) Y = —2ryp = 10-* sec.~?. 


The corresponding effective centrifugal acceleration even near the rim of the 
galaxy would still have the exceedingly small value 


(25) (?R’ = 10-7! cm. sec.~, 


which is 11 orders of magnitude smaller than the ordinary centrifugal accelera- 
tion due to the rotation with angular velocity w’. 


5. REMARKS ON THE PRINCIPLE OF EQUIVALENCE 


In the preceding discussion it has been assumed that all forms of matter and 
energy are sources of a universal gravitational field, which is attractive for 
resting matter. This amounts to assuming the validity of the so-called prin- 
ciple of equivalence, which states that any form of matter which exhibits 
inertia is at the same time source of a gravitational field of a strength pro- 
portional to its inertia and may therefore be assigned a gravitational mass, 
which is always taken to be positive and proportional to its inertia. In par- 
ticular, any piece of anti-matter in our universe has been assumed to have 
positive gravitational mass, because there is no experimental evidence to date 
which would indicate any difference between the inertias experienced by a 
particle and by its anti-particle. Within the framework of the model theory 
discussed above it is thus in principle possible to assume that any fraction of 
the universe consists of anti-matter without changing any of the conclusions 
that have been reached so far. 

Any tampering with the principle of equivalence would require a drastic 
change in the field equations to make realization of Mach’s Program feasible, 
and would entail a radically different approach. Recently Morrison and Gold 
(1956) have put forward the interesting hypothesis that anti-matter has 
negative gravitational mass, and have envisaged a universe in which the 
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number of galaxies made up out of anti-matter is equal to the number of 
galaxies made up out of matter. It is clear that one cannot realize Mach’s 
Program on the basis of the field theory of gravitation used here if the model 
of Morrison and Gold is adopted, since, except for tiny local fluctuations of the 
order of magnitude of the contribution to the inertia by a single galaxy, there 
will be cancellation of the gravitational effects of matter and anti-matter con- 
tributed by the distant galaxies whose influence dominates local inertial 
behavior. One may thus say, on the basis of the gravitational theory used here, 
that if anti-matter has indeed negative gravitational mass, then the very 
existence of inertia and the facts incorporated in Newton’s second law are 
manifestations of the preponderance of matter over anti-matter in our universe. 

Morrison (1957) has suggested that in addition to the positive gravitational 
“rest”? mass of matter and the negative gravitational “‘rest’’ mass of anti- 
matter there is a third kind of gravitational mass which has the same sign 
for both matter and anti-matter and which is ‘‘kinetic’’ in origin so that it 
vanishes for resting matter or anti-matter. Upon adoption of Mach’s Program 
one may then attempt to account for the inertia of both matter and anti- 
matter by interaction with the ‘‘kinetic’’ mass of the universe, to which the 
universal recessive motion of distant galaxies may conceivably make a domi- 
nant contribution. One can, of course, not rule out such possibilities, but it 
appears that to this end a theory of gravitation is required which leaves the 
realm of invariance under Lorentz transformations, and for which therefore 
the laws of conservation of energy, momentum, and angular momentum are 
put into question as far as the gravitational field is concerned. 
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HORIZONTAL MOVEMENTS OF AURORA' 


J. S. Kim ann B. W. CURRIE 


ABSTRACT 


Measurements on the drifts of auroral forms at three stations to the south 
of the auroral zone in West-Central Canada failed to show evidence of a 
motion due to the earth’s rotation relative to fixed excitation patterns in space. 
The distribution and the magnitude of the speed of auroral structures parallel 
and normal to the geomagnetic meridians are substantially the same as for 
non-luminous ionic irregularities observed by radio methods. The speed of an 
auroral arc or band parallel to the geomagnetic meridian is apparently constant. 
Speeds increase with geomagnetic activity, particularly in an east-west direc- 
tion. There is no definite indication of a reversal of the east-west motions 
close to local midnight. Speeds and directions show no characteristic variations 
with the time of night. 


INTRODUCTION 


Numerous radio investigations in recent years have demonstrated the 
existence and movement of ionic irregularities in the upper atmosphere 
(Ratcliffe 1955). The origins of these irregularities are still uncertain. Some are 
unquestionably due to the impacts of extraterrestrial matter such as meteors, 
particles ejected from the sun, and galactic dust attracted towards the sun. 
Some may originate from turbulence—the turbulence causing electronic 
distributions that are normally uniform to become irregular, or producing 
irregular distributions of air density that become apparent when the air is 
ionized. Some may be due to electrical discharge processes within the atmos- 
phere such as discharges upward into the ionosphere from thunderstorms. 
The causes of the observed movements are even less certain. They may result 
from the movement of the surrounding neutral air through meteorological 
processes, or the action of electric fields in combination with the steady 
magnetic field of the earth, or wave-motions in the distributions of electrons, 
or a combination of these causes. For these reasons, it is of interest to investi- 
gate the horizontal motions of aurora to see what features, if any, they may 
have in common with the motions of ionic irregularities that are observed by 
radio methods. 

The movements of auroral ionization, as observed visually, may be divided 
into two types. One is generally very rapid and takes place within the body 
of the display. The other is relatively slow and involves displacements of 
displays as a whole. The first is often more apparent than real, auroral struc- 
tures appearing and disappearing in rapid succession and thus giving the 
impression of waves of luminosity travelling either horizontally or toward 
the observer’s zenith. The second is evident most frequently as a displace- 
ment of an arc or band normal to its lateral extent. To the south of the zone 
of maximum auroral occurrence, an arc or band will often appear in the 


‘Manuscript received October 7, 1957. 
Contribution from the Institute of Upper Atmospheric Physics, University of Saskatchewan, 
Saskatoon, Sask., Canada. 
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observer's northern sky and then gradually move southward without any 
marked change in its configuration. This may be followed by a succession of 
slow retreats and advances. Sometimes the arc or band disappears as it 
approaches the observer's zenith, and is followed by the appearance of a 
new arc or band in the north that moves southward. Similar motions are 
observed to the north of the auroral zone, the primary movement now appear- 
ing to be northward. Motions parallel to the lateral extent of a display also 
occur, but these are difficult to detect. Detection depends on the presence of 
persistent convolutions and structures that cannot be confused with the 
more transient features that characterize the first type of movement. Slow 
rotations of a display, clockwise or counterclockwise, are also occasionally 
observed. These again may be more apparent than real, variations in the 
lateral direction of an arc as it moves eastward or westward being interpreted 
as rotations. It is movements of the second type that are considered in this 
paper. 

Quantitative investigations of these movements have only recently become 
possible through the development of all-sky cameras (Meinel and Schulte 
1953; Stoffregen 1955; Meek 1955; Park 1957). The sky is imaged in a convex 
or parabolic mirror with its axis vertical. This image is photographed (usually 
once every minute) by a movie camera placed directly above the mirror. A 
clock or a digital counter controlled by a clock is located in the field of view 
of the movie camera to give the time. Because exposure times are seldom long 
enough to record the star background, altitudes and azimuths of points on a 
photograph are obtained from suitably-located markers or lamps placed in 
the field of view of the mirror and photographed along with the sky. The 
markers sometimes consist of a hemispherical grid of wires placed symmetric- 
ally about the mirror. Often they consist of nothing more than the rods 
supporting the camera and short crosspieces attached to them at heights 
corresponding to selected angular altitudes as imaged by the camera. For the 
accurate determination of angular altitudes and azimuths of auroral points, 
the photographs must be projected onto a grid of concentric circles and radial 
lines. The grid is constructed usually by reference to star photographs taken 
with the camera so as to eliminate errors due to aberrations in the optical 
systems. The images of the markers serve to adjust the magnification of the 
projector, and to locate the projected image properly on the grid. 

The movements of an auroral display are seen readily by running the film 
from an all-sky camera at a moderate rate through a projector. This is par- 
ticularly necessary for east-west motions. An observer looking directly at 
an auroral display seldom recognizes these motions, possibly because persistent 
features that serve as indicators of the motions are confused with the more 
transient features of the display. Because of the relatively long exposure 
times required (20 seconds for some types of cameras), the transient features 
either do not appear or are smeared by the auroral motions into a picture 
that shows only the limiting edges of the auroral display. 

Certain assumptions are necessary in order to measure horizontal displace- 
ments and the rates at which they take place. One is that the lower edge of 
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arcs and bands occurs at a constant height above the earth. Numerous paral- 
lactic measurements of auroral heights indicate that this is generally the case, 
and that a height of 100 km. is a reasonable one to assume (Currie 1955; 
Jensen and Currie 1953). Another is that the components of the displace- 
ments most likely to be of value for comparative purposes are ones parallel 
and normal respectively to the geomagnetic meridians, since auroral character- 
istics such as the isolines of occurrence and the lateral extent in an east-west 
sense coincide more closely with the geomagnetic than the geographic lines 
of latitude (Jones and Currie 1941). Only displacements of arcs and bands 
parallel to a meridian, either geographic or geomagnetic, can be measured 
with any considerable frequency; and to a first approximation these are normal 
to the arcs and bands in West-Central Canada. 


OBSERVATIONAL DETAILS 


Motions parallel to the geomagnetic meridians were measured on all-sky 
photographs taken at Saskatoon, Flin Flon, and Uranium City during January, 
February, and March, 1956. The geographic latitude and longitude (¢, A), 
the geomagnetic latitude and longitude (#, A), and the angle between the 
geographic and the geomagnetic meridian (Ww) for each station are given in 
Table I. Uranium City is on the southern edge of the zone of maximum 
auroral occurrence for the Northern Hemisphere. The photographs were 


TABLE I 
GEOGRAPHIC LATITUDE AND LONGITUDE (@, A), GEOMAGNETIC LATITUDE 
AND LONGITUDE (#, A), AND THE ANGLE WY BETWEEN THE GEOGRAPHIC 
AND GEOMAGNETIC MERIDIANS FOR SASKATOON, FLIN  FLON, 
AND URANIUM CITY 


d, r, ?, A, Y, 


Station IN. oW. IN. °W. fe 
Saskatoon 52.1 106.6 60.5 48.9 13.8 
Flin Flon 54.8 101.8 63.8 43.5 13.6 

59.6 108.8 67.7 56:6 19.1 


Uranium City 





taken on 16-mm. film at 1-minute intervals with about 20-second exposure 
times. The period covered by the photographs was one with a predominance’ 
of arc and band-like structures. Comparison of the relative occurrences of 
auroras at each station (number of nights with aurora expressed as a per- 
centage of the number of nights when aurora could have been photographed 
if present) showed that the data were approximately representative of the 
region and of the time of year. Auroral occurrence increased with geomagnetic 
latitude at about the expected rate (Vestine 1944). Auroral occurrence at all 
three stations was substantially greater for January than for February, con- 
trary to long-time averages for regions to the south of the auroral zone. 
However, auroral occurrence increased rapidly during March in accordance 
with the approach to the usual spring maximum. 

Motions normal to the geomagnetic meridians were measured on all-sky 
photographs taken at Flin Flon throughout 1956. Owing apparently to in- 
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creased solar activity, the character of the displays changed during the later 
part of 1956, structures suitable for detecting lateral motions occurring more 
frequently. The additional months of data provided sufficient measurements 
to be of some value for statistical purposes. 

The photographs were projected onto a special grid for each station. This 
consisted of a system of geographic latitude and longitude lines, centered on 
the station and extending outward to cover the region over which auroras 
were photographed. The lines were spaced according to the assumption that 
the lower edge of all auroral forms is at a height of 100 km. above the earth. 
The latitude and longitude of the point directly below any selected point on 
the lower edge of an auroral form could be read directly from the grid. Absolute 
displacements could then be computed from the latitude and longitude readings 
corresponding to the positions of the selected auroral point on successive 
photographs. Actually, the grids were used largely to check on the direction 
and uniformity of motion of the auroral forms. A line corresponding to the 
geomagnetic meridian for the station was placed on its grid. It was marked 
with a scale of distances from the station. Measurements of displacements 
parallel to the geomagnetic meridians were made only on arcs and bands that 
were above the geomagnetic meridian for the station and that moved so that 
their positions on successive photographs were substantially parallel to each 
other. The grid as a whole was necessary to check this point. Measurements 
of displacements normal to the geomagnetic meridians were more complicated, 
since computations of the displacements had to be made from the latitude 
and longitude readings. 

Displacements were ordinarily measured for 4-minute intervals, and in 
most cases on auroral structures at least 20 degrees above the horizon. While 
shorter time intervals and lower altitudes could have been used, the errors 
would have been appreciable because of the small size of the photographs 
and the rapid decrease in magnification of horizontal distances toward the 
edges of the photographs. 

RESULTS 
1. Motions Parallel to the Geomagnetic Meridians 

Several general features of the motions of arcs and bands parallel to the 
geomagnetic meridians, and hence approximately normal to their lateral 
extent, are indicated by the measurements. The speed of each arc or band 
is constant to within the accuracy with which the measurements could be 
made. An are or band is rarely stationary, although the direction of motion 
of some reverses. Motions in a northward or southward sense are not related 
to the time of night. The higher speeds, especially in a southward sense, 
occur generally at times of considerable magnetic disturbance. 

The deduction concerning the constancy of the speed of particular arcs 
and bands parallel to the geomagnetic meridian was based on measurements 
for 77 independent arcs and bands as photographed from Saskatoon, 587 
from Flin Flon, and 502 from Uranium City. Of these, 49, 335, and 328 respec- 
tively were moving southward. For all except the fastest moving forms, the 
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speeds were measured for;at least two successive 4-minute intervals; and in 
numerous cases for five or'more. One persistent arc moved slowly enough for 
19 such measurements to be made; another for 15 measurements. 

About one-third of the arcs and bands reversed their direction of motion. 
Whether or not they came to rest and remained so for a short period of time 
before starting in the opposite direction is uncertain. A gradual slowing before 
and speeding up after the reversal was not observed, but this may have 
been because of the l-minute intervals between photographs. Very occasion- 
ally a series of photographs showed an arc that was apparently stationary, 
but these were for low altitudes where small displacements were difficult 
to detect. 

The percentage occurrences of arcs moving southward, northward, and 
reversing direction at cach of the three stations for the period before midnight 
were approximately the same as for the period after midnight. There was 
some indication that auroras during the hour or so immediately after sunset 
and again before sunrise have motions that are predominantly southward. 
The data were insufficient to be certain about this point. The speeds also 
showed no significant differences for different periods of the night, although 
they were usually greater at times close to sunset and sunrise than at other 
times. 

The largest speeds invariably occurred at times with much magnetic dis- 
turbance. A correlation of speeds with magnetic disturbance indices would 
have required considerably more data, since the number of arcs and bands 
suitable for measuring displacements parallel to the geomagnetic meridians 
decreases as the magnetic activity increases. Rayed and other structures 
which change their form rapidly with time begin to predominate, and the 
number of reliable measurements of displacements that can be made gets 
small. 

Table II gives a tabulation of the observed speeds parallel to the geomag- 
netic meridians, grouped according to 3-km./min. intervals of speed. Because 
of the small number of cases with speeds above 39 km./min., these are shown 
as a single group. The highest observed speed was about 50 km./min., once 
southward at Saskatoon and twice northward at Flin Flon. It is evident from 
the table that the number of cases with southward motion exceeds the number 
with northward motion for all ranges of speed except the one for speeds less 
than 3 km./min. Obviously, there is a resultant displacement of auroral 
ionization southward. The average southward speeds of arcs and bands at 
Saskatoon, Flin Flon, and Uranium City from these observations are re- 
spectively 10.6, 9.4, and 10.4 km./min., and average northward speeds, 
8.5, 8.0, and 9.3 km./min. Averages computed using the 4-minute values, 
irrespective of the arcs or bands on which they were made, give approximately 
the same speeds, namely 9.7, 8.0, and 9.3 km./min. southward, and 8.8, 
6.3, and 8.9 km./min. northward. Apparently there is no systematic variation 
of speeds with latitude. Applying the usual statistical tests, the differences 
between the southward and northward speeds are significant for each station. 
The differences are appreciably smaller than would be expected for an excita- 
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TABLE II 


SPEEDS OF AURORAL ARCS AND BANDS PARALLEL TO THE GEOMAGNETIC MERIDIANS AS OBSERVED 
AT SASKATOON, FLIN FLON, AND URANIUM CITY 


The observations are grouped according to 3-km./min. ranges of speed. Each row includes 
occurrences of the highest indicated speed for the range but not the lowest. The first quantity 
in the columns under the names of the stations gives the number of 4-minute occurrences of 
the indicated speed; the second quantity, which is in brackets, the number of independent 

arcs and bands with the indicated speed. 








Saskatoon Flin Flon Uranium City 
Speed, —_—__—-—-— —_-—__— $$$ —~ 
km./min. South North South North South North 
0-3 3(1) 11(2) 149(45) 188(46) 78(20) 51(17) 
3-6 47(11) 11(4) 235 (82) 150(71) 187(62) 83 (44) 
6-9 18(7) 15(9) 134(61) 98(49) 142(63) 77(33) 
9-12 28(9) 11(4) 95(51) 66(38) 141(66) 63(33) 
12-15 21(10) 2(1) 77 (37) 31(21) 102(43) 24(16) 
15-18 5(3) 0 42(21) 16(8) 73(36) 14(10) 
18-21 7(2) 3(2) 18(11) 17(9) 25(16) 17(12) 
21-24 2(2) 2(2) 18(11) 6(4) 17(10) 8(3) 
24-27 2(1) 1(1) 13(9) 2(2) 2(2) 3(2) 
27-30 2(2) 1(1) 3(1) 2(1) 4(4) 9(3) 
30-33 1(1) 5(4) 1(1) 
33-36 2(2) 1(1) 
36-39 1(1) 2(2) 
39 1(1) 2(1) 3(2) 3(2) 2(1) 
792(335) 580(252) 778(328) 350(174) 


Totals 136(49) 60(28) 


tion pattern fixed in space relative to the earth. For arcs and bands lying along 
the geomagnetic parallels of latitude, the expected speeds southward due to 
the earth’s rotation are respectively 4.1, 3.8, and 4.6 km./min. at Saskatoon, 
Flin Flon, and Uranium City. If the random portions of the speeds southward 
and northward have the same average values for a large number of observa- 
tions the resultant drift southward along the geomagnetic meridian at each 
station is obtained by dividing the observed differences by two. These are 
about 1.0, 0.7, and 0.6 km./min. respectively. 


2. Motions Normal to the Geomagnetic Meridians 

Several general features of the motions of auroral structures normal to the 
geomagnetic meridians were indicated by the measurements. These are less 
certain than ones for the motions parallel to the geomagnetic meridians be- 
cause of the limited number of cases. When structures suitable for measuring 
absolute displacement occur, the motions are predominantly in an east-west 
sense and occur at much greater speeds than for the south—north motions. 
Motions in an east-west sense are rarely detected except at times of consider- 
able magnetic disturbance; and the more disturbed the magnetic field, the 
larger are the speeds. Speeds are apparently independent of the time of night, 
and on an average are the same for eastward and westward motions. Westward 
motions occur more frequently before midnight and eastward motions after 
midnight. The speed of particular structures appears to remain constant; 
but this is difficult to verify because of their relatively short durations. 

Measurements were made on 90 independent auroral structures, 35 exhibit- 
ing eastward motion and 55 westward motion. The average duration of each 
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structure was about 5 minutes. The average speed eastward of the structures 
was 28.9 km./min., and westward, 28.7 km./min. Average speeds computed 
from the total displacements and the total durations were 26.3 km./min. 
eastward and 26.4 km./min. westward. Again the observations show no evi- 
dence of an excitation pattern fixed in space relative to the earth. 

Table III shows a tabulation of the average departures from normal of 
the magnetic field intensity at Flin Flon corresponding to various ranges of 
speed. The magnetic values were computed from the average values of AX, 
AY, and AZ (X, Y, and Z are the components northward, eastward, and down- 
ward of the earth’s field as recorded by a three-component magnetic vario- 


TABLE Ill 


AVERAGE DEPARTURES FROM NORMAL OF THE MAGNETIC FIELD INTENSITY AT FLIN FLON 
CORRESPONDING TO VARIOUS RANGES OF SPEEDS OF AURORAL STRUCTURES NORMAL TO THE 
GEOMAGNETIC MERIDIANS 





Ranges Average Average departure 
of speed, Number speed, field intensity, 
km./min. of cases km./min. gammas 

0-10 5 5.4 24 

10-20 21 14.5 102 

20-30 19 25.4 142 

30-40 17 34.1 145 

40-50 8 46.6 182 

50-60 8 53.7 173 

60-70 2 64.6 125 


70-80 1 71.0 113 


meter), for the time intervals during which the displacements occurred. 
Nine of the 90 cases are not included because of incomplete magnetic records. 
Fig. 1 shows a plot of the average speeds as abscissae against the average 
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Fic. 1. Plot of the average departures from normal of the magnetic field intensity during 
the east-west motions of auroras for 10 km./min. ranges of speed. The vertical lines show 
the extreme variation of the magnetic values included in each average value. 
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departures of the magnetic field intensity as ordinates. Although the variance 
is large, there is unquestionably an increase of speed as the magnetic distur- 
bance increases. A relationship between the direction of motion and the magni- 
tude of the disturbance was not evident. 

Table IV shows the numbers of occurrences of eastward and westward 
motions and the corresponding average speeds for different hours of the night. 
Westward motions were characteristic of the early morning hours, eastward 
motions of the late morning hours, and motions in both directions of about 
a 4-hour period centered on midnight. The speeds do not indicate a gradual 
slowing down of the westward motions prior to a change to eastward motions, 
as might be expected if motions in the evening hours are generally opposite 
to those for the morning hours. 


TABLE IV 


THE DIURNAL OCCURRENCES AND SPEEDS OF 
EAST-WEST MOTIONS OF AURORAL STRUCTURES 
AT FLIN FLON 


Occurrences Speeds, km./min. 
Time, —_—__—__—_——— - 
M.S.T. East West East West 


18-19 


1 23.0 
19-20 l 34.3 
20-21 4 26.0 
21-22 8 38.7 
22-23 4 9 22.6 23.4 
23-2: 6 15 32.3 30.6 
00-01 13 10 27.0 29.8 
01-02 5 7 35.3 20.0 
02-03 4 28.3 
03-04 3 35.4 


04-05 


DISCUSSION 

Auroras are assumed generally to be excited by charged particles that 
travel in streams from the sun and penetrate into the earth’s atmosphere. 
Because of the interaction between the charged particles and the earth's 
magnetic field, the latitudes where auroras occur depend on the speed of the 
particles. Variations in speeds of the particles in the solar streams, and the 
distortions of the earth’s magnetic field by the ring current that is assumed to 
develop outside the earth and in its equatorial plane during geomagnetic 
storms, are presumably responsible for the sporadic occurrence of auroras 
over regions in high latitudes, and the apparently random ways in which 
displays move. The only systematic motions that should be observed are ones 
superimposed on these due to the earth’s rotation relative to the streams of 
exciting particles. The failure of these observations to show any definite 
indication of a motion due to the earth’s rotation confirms the earlier and 
less extensive observations on east-west motions at Yerkes Observatory 
(Meinel and Schulte 1953) and at Saskatoon (Meek 1954). It is unlikely 
that the small resultant drift southward of the auroral ionization along the 
geomagnetic meridians is due to sampling and errors in measuring the speeds. 
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The average orientation of arcs and bands in West-Central Canada (Jensen 
and Currie 1953) is such as to have given appreciably greater values of the 
southward drift, if it had been due to the earth’s rotation. 

It is of some interest to compare the motions of the auroral luminosities 
with those of the non-luminous ionic irregularities detected in the E and F 
regions by radio methods. The speeds of auroral forms and structures are of 
the same order of magnitude as the speeds of the ionic irregularities (Chapman 
1953; Briggs and Spencer 1954; Toman 1955; Maxwell 1955; Dagg 1957; 
Harang and Pedersen 1957). Histograms of the frequencies of occurrence of 
the speeds tabulated in Table II (see Fig. 2) would approximate in many 
cases to the corresponding histograms for the ionic irregularities, for example, 


OCCURRENCE 


PERCENTAGE 





o 3 6 9 ? 15 182i 24° er = 30.3 
SPEED — KM/MIN 
Fic. 2. Histograms of the percentage occurrences of 3-km./min. ranges of speed of north- 
ward and southward moving auroras. The larger speeds are more frequent for southward 
moving auroras. 


Fig. 16 by Toman for F-region drifts. The directions of the E and F region 
drifts, observed by radio methods, are usually more regular and change 
systematically with the time of day and the season. This may be because the 
motions of the ionic irregularities can be followed continuously, and some 
of the radio methods permit more exact determinations of direction of motion. 
The ionic irregularities also experience sudden changes of direction of motion, 
for example, Fig. 1 by Briggs and Spencer. Recent observations on E-layer 
drifts at Kjeller, Norway, by Harang and Pedersen, indicate a prevailing 
motion toward the southeast in winter, the southward component of which is 
about the same as for the drift of the auroral arcs and bands toward the south. 
Since aurora occurs predominantly in the £& region, it might be anticipated 
that the aurora would be subject to any relatively-constant drifts characteristic 
of the £& region. On the other hand, if aurora is due to a discharge mechanism 
due to electrical fields that develop between the / and F regions, the auroral 
drifts might well be influenced by drifts in the / region. While it is difficult 
to see how the discharge fields could develop because of the high electrical 
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conductivity along the lines of the earth’s magnetic field, meridional winds 
deduced for the F region are southward to the south of latitude 60° N. (approxi- 
mately the position of the auroral zone in West-Central Canada), and north- 
ward to the north of latitude 60° N. (Pant 1956). Visual observations made to 
the north of the auroral zone during 1932-33 at Chesterfield, Coppermine, 
and Cape Hope’s Advance indicate that an investigation of auroral drifts, 
similar to this one, would show a northward drift over latitudes to the north 
of the auroral zone. While confirmation of the visual observations by all-sky 
camera studies could not be taken as positive proof of a control on auroral 
drifts by processes taking place in the F region, it would suggest a relation- 
ship between winds of meteorological origin and the average drifts of auroras. 

The increase in speed of the non-luminous ionic irregularities with in- 
creasing Magnetic activity is now a well-established phenomenon, especially 
for the F region. For example, see Fig. 2 by Maxwell and Fig. 1 by Dagg for 
F-region drift speeds plotted against geomagnetic K indices. Since the auroral 
speeds (especially of structures which usually extend upward beyond the 
E, region) also show an increase with geomagnetic activity, the possibility 
exists that the same mechanism controls the speeds of both the ionic irregu- 
larities and the aurora. 

The observations by Meinel and Schulte at Yerkes indicated that the motion 
was westward in the evening, reversed during the midnight hours, and was 
eastward in the morning. Meek’s observations at Saskatoon showed eastward 
motions taking place at various times throughout the night, westward motions 
up to 2:00 a.m., and no evidence of a reversal. In general, the Flin Flon 
observations agree with the ones made by Meek. While only westward motions 
were observed in the early evening and eastward motions in the late morning, 
motions in one direction were relatively as frequent as in the other for the 
four hours centered on midnight and showed no evidence of a slowing down 
close to midnight. The failure of the Saskatoon and Flin Flon observations 
to agree with the ones at Yerkes is disappointing, since investigations by 
Chapman, Vestine, Martyn, and others of geomagnetic disturbance features 
indicate the possibility of current systems flowing at ionospheric levels and 
reversing direction near midnight, the current systems being induced by winds 
at ionospheric levels through a dynamo action. 

Three distinct types of clouds of abnormal ionization occur in auroral 
and polar regions as observed by ionosonde equipment (Hagg and Hanson 
1954). Two of these, clouds of sporadic E and clouds accompanied by radio 
blackouts, appear to be associated with aurora. The first has speeds and speed 
distributions comparable to those for aurora, but occurs more frequently at 
Resolute Bay (74.7° N.) than at Baker Lake (64.3° N.). The frequency of 
occurrence of aurora at the former is very much less than at the latter station. 
The diurnal distribution of the first also differs notably from that for aurora. 
The second type is observed at stations close to the auroral zone, and has a 
diurnal distribution comparable to that for aurora. It has average speeds of 
75 km./min., and maximum speeds of the order of 250 km./min. Radio black- 
outs are often an aftermath of such occurrences, and these may be wide- 
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spread and persist for many hours. The second type is apparently associated 
with the penetration of ionizing particles into the atmosphere, and is usually 
accompanied by auroral displays, but the speed, the extent of the clouds, 
and the persistence indicate that it is not aurora. Again it is interesting to 
speculate whether or not widespread ionization can occur at and below auroral 
levels with visual aurora taking place intermittently only at particular loca- 
tions. 

It is believed that these observations on auroral motions, and the discus- 
sion of their possible relationships to other phenomena involving ionic irregu- 
larities in the upper atmosphere, are sufficient to show that studies of auroral 
motions by sequences of all-sky camera photographs are worth while. Eighteen 
of these cameras, including nine 35-mm. cameras designed by the National 
Research Council that make possible more accurate measurements on aurora 
than the usual 16-mm. camera, are operating in Canada during the Interna- 
tional Geophysical Year. While they are intended primarily to record the 
occurrence of aurora and the development of auroral displays, an investigation 
of auroral movements both inside and outside the zone of maximum auroral 
occurrence from their photographs may contribute significantly toward an 
explanation of the process by which aurora is excited. 
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ROTATIONAL SPECTRA OF RbF BY THE ELECTRIC RESONANCE 
METHOD! 


H. Lew, D. Morris,” F. E. GEIGER, JR.,? AND J. T. E1sinGER‘ 


ABSTRACT 
Transitions between the J = 0 and J = 1 rotational states of RbF have 
been measured by means of the molecular beam electric resonance method. 
The following rotational constants have been determined (all frequencies in 
Mc. /sec.): 


Rb®F Rb*F 
Yu = B. 6315.543+0.009 6289 .002+0.010 
—-Yn=a, 45 .638+0.017 45 .350+0.017 
Yo = Ye 0.094+0.007 0.094+0.007 
B, 6315 .639+0.036 6289 .098+0 .036 
The quadrupole interaction constants —eqQ/h in the J = 1 state are found to 
be 
Rb®F v=0 70.340+0.030 Mc./sec. 
v=1 69 .543+0.060 Mc./sec. 
Rb*F v=0 34.024+0.040 Mc./sec. 


eqQ/eqQ* = 2.06730. 0030 


The equilibrium internuclear distance is r, = (2.26554+0.00005) X10-* cm. 
The electric dipole moment of Rb*®F in the v = 0 state is u = (8.80+0.10) 
X107-'8 e.s.u. The mass ratio of the Rb isotopes is W5;/Ms; = 0.9770148+ 


0.0000052. 
I. INTRODUCTION 


Over the last few years, pure rotational spectra in almost all the alkali 
halides have been obtained either by microwave absorption spectroscopy 
at high temperatures (Honig, Mandel, Stitch, and Townes 1954) or by the 
molecular beam electric resonance method (e.g. Lee, Fabricand, Carlson, 
and Rabi 1953). From the extensive data general trends in such quantities 
as interatomic distances and electric dipole moments have been observed 
and empirical relations derived for them. As a step towards completing the 
picture for the alkali halides, it was decided to study RbF in the electric 
resonance apparatus of this laboratory. This molecule was one of the few 
that had been left undone by the microwave spectroscopists (Honig et a/. 1954) 
because of experimental difficulties. The other alkali halides yet to be studied 
for rotational spectra are LiF, KF, NaF, and LiCl. Rubidium fluoride was 
chosen for the present study because the frequencies required would be the 
lowest among the molecules mentioned and because the necessary microwave 
equipment for the estimated transition frequencies was already on hand. 
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Rubidium fluoride was first studied by the electric resonance method by 
Hughes and Grabner (1950), who observed transitions between Stark sub- 
levels belonging to a given rotational level. They obtained values for the 
quadrupole interaction constant egQg» for different rotational states and 
different vibrational states and values for the so-called Jy.J interaction 
constant where /, is the spin of the fluorine nucleus. The present measure- 
ments, in addition to yielding quadrupole interaction constants which are 
in agreement with those of Hughes and Grabner, yield values for the molecular 
rotational constants B, and a, and a value for the electric dipole moment u. 
The J.J interaction, however, was not resolved. 

Prior to a search for transitions an estimate of B, was made through the 
use of internuclear distances tabulated by Honig ef a/. (1954). These authors 
give a calculated internuclear distance for RbF of 2.246 A with a probable 
accuracy of about one per cent. Thus B, should be around 6450+125 Mc./sec. 
or the frequency of a (J = 1)—(J =0) transition about 12900+250 
Mc. /sec. 


Il. THEORY 

For the analysis of the rotational spectra observed in the present experi- 
ment, aside from the hyperfine structure, Dunham's (1932) expressions for 
the energy levels of a vibrating rotator are used. In order to avoid excessive 
reference to the original or other papers, some of the appropriate relations 
are reproduced here. In terms of Dunham’s coefficients the frequencies of a 
(J = 1) <—(J = 0) transition for the first four vibrational states are (Lee 
et al, 1953) 
fo = 2Vut4¥oot+ Yut(1 /2) Yo, 
fi = 2VYu+4V¥o+3¥ut (9/2) ¥ a, 
fe = 2Vot+4 Voot+5 Vit (25/2) Yn, 
fs = 2¥out+4¥ot+7 Yu+ (49/2) Ya. 


The various Y,; are related to more familiar molecular constants as follows: 


(1) 


You = B,, Yo = —D,, Yu = —@e, Vu = Ye- 


Since the determinant of the coefficients in (1) is of rank 3, only three coeffi- 
cients may be determined even though more than three f, may be available. 
The smallest of the coefficients, Yo, is usually calculated from the relation 
Yo = —4B,3/w,? provided w, is known from other work. 

Before Eq. 1 can be used to treat the experimental data, the hyperfine 
structure of the J = 1 state must be taken into account. The theory for this 
has been given by Hughes and Grabner (1950). These authors found that, 
in their study of transitions between various sublevels belonging to the J = 1 
state, it was sufficient to assume two internal interactions in the molecule: 
one an interaction between the electric quadrupole moment of the Rb nucleus 
and the gradient of the electric field acting on it and the other a magnetic 
interaction between the fluorine nuclear magnetic moment and the field 
arising from the rotation of the molecule (the J.J or spin-rotation interaction). 
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These two interactions split the J = 1 level into six components as shown 
in Figs. 1 and 2. In the figures, the various quantum numbers have the follow- 
ing significance: F; = 1,;+J where J, is the spin of the Rb nucleus and 





F F F, F 
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Fic. 1. Energy levels of Rb®F. Fic. 2. Energy levels of Rb*F. 


F = I,+F, where J; is the spin of the F nucleus. There are theoretically six 
observable transitions between J = 1 and J = 0 but, since the J.J splitting 
is not resolved in the present experiment, only three distinct lines are observed 
and these are labelled A, B, and C. The J.J splittings are of the order of 
10 ke./sec. (Hughes and Grabner 1950). The patterns shown in Figs. 1 and 
2 are, of course, for zero external electric field. 

The electric dipole moment u of the molecule may be determined by observ- 
ing the shift in a transition frequency on the application of a static electric 
field. For RbF the most convenient line for this purpose is line C. This line 
involves the levels (J = 0) and (J = 1, F; = %). By Eq. 2 of Hughes (1947) 
the shift in the J = 0 level for weak electric fields E is simply 


(2) AW = — p?E?/B,. 


By Eq. 13 of Fano (1948) or similar expressions elsewhere (e.g. Townes and 
Schawlow (1955), Eqs. 10-31) the upper level of line C splits into two com- 
ponents which are separated from the field-free level by 


1 Ee 
‘ va heen 
(3) AW (+35 B. 
where the upper sign applies for the substates |M,y,| = } and the lower 
sign for |My,| = 3. At sufficiently small fields these two Stark levels are un- 


resolved and the equal and opposite shifts of the two substates result only 
in a broadening of the line. Hence the frequency shift in line C is simply 
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that given by Eq. 2. This very simple relationship between the electric dipole 
moment and easily measured quantities makes the line C particularly suitable 
for the determination of the dipole moment. 


III. APPARATUS 


The apparatus is of conventional design, basically that described by Hughes 
(1947) but in the modification first described by Lee et al. (1953). A diagram 
of the apparatus giving some of the important dimensions is shown in Fig. 
3. The vacuum tank is of non-magnetic stainless steel plate, 0.475 cm. thick, 
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Fic. 3. Schematic plan view of apparatus. 


rolled into a cylinder and welded. It is divided into the usual three chambers 
and evacuated by three oil-diffusion pumps. 

The oven is a vertically mounted cylindrical steel tube closed at the bottom 
and provided with a plug at the top. It has an inside diameter of 1.59 cm. 
and a wall thickness of 0.025 cm. The total depth of the oven is 5 cm. and 
the center of the slit is located 3.2 cm. from the bottom. The width and height 
of the slit are respectively 0.0076 cm. and 1.27 cm. The heating elements 
for the oven are molybdenum coils enclosed in thin ceramic tubes mounted 
parallel to the axis of the oven and arranged in a circle around it. This 
assembly is surrounded by two radiation shields, the outer one of which is 
water cooled. 

The brass pole-pieces of the A and B fields are of identical cross section 
but of different lengths. The cross section is as shown in Fig. 4. The various 
design constants are given in Table I. The beam position is exactly equi- 
distant from the two outer faces of the pole pieces so that the latter may be 
reversed in their mount and hence the gradient reversed without a need for 
realignment of the system. 
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~~] 
or 





W = 3.754 cm (B field) 
= 3.147 cm (A field) 


Fic. 4. Pole pieces of A and B fields, sectional view. 


TABLE I 
DESIGN CONSTANTS OF DEFLECTING FIELDS 
(A and B fields of identical cross section) 


Separation of equivalent parallel wires 2a = 1.000 cm. 
Radius of convex piece R; = 0.508 cm. 

Radius of concave piece R2 = 0.564 cm. 

Separation between centers of curvature d = 0.345 cm. 
Beam position = 0.1 cm. from apex of convex piece 
Field intensity/voltage = 2.888 (volt /cm.) /volt 
Gradient /field = 1.981 cm.~? 


The C field consists of two parallel plates, one of which (Item 1 in Fig. 5) 
is of aluminum deposited on glass and the other (Item 2) of brass. The brass 
piece contains a slit 0.038 cm. high and 4.75 cm. long in the direction of the 
beam. Microwave energy is fed into the C region by way of this slit. The 


OO) @ 


(b) 


Fic. 5. C field. 
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slit couples to an external waveguide through a 55° bend and a taper section. 
The bend was necessitated by an obstacle in the side of the vacuum tank 
directly opposite the C field. Although the brass plate was ground optically 
flat prior to assembly it is not known how flat it remained after assembly. 
The microwave energy in the waveguide is in the TE; mode. 

The use of an aluminized glass plate as one of the surfaces of the C field 
was decided on because of the ease with which such a plate could be made 
flat and the ease of having the active surface divided into two mutually 
insulated regions as shown in the lower half of Fig. 5. The aluminum surface 
was merely scored with a pointed instrument. The purpose of this division 
was to permit low frequencies (a few megacycles) to be applied to the C 
region for the observation of AM, transitions. The large contact potential 
between aluminum and brass, however, had been inadvertently overlooked 
and was not recognized until late in the measurements. Fortunately this 
potential introduced a negligible correction in the measurements at zero field, 
as will be shown. 

The aluminized surface was coated with a thin film of silicon dioxide to 
protect it from abrasion. Quartz spacers between the plates provide a spacing 
of 0.69065 cm. 

The collimator consists of two quartz semicylindrical jaws mounted at a 
fixed separation of 0.152 cm. between apexes. The radius of curvature of the 
jaws is 0.318 cm. By rotating the unit about an axis coincident with the 
mid-line of the opening between the jaws the effective slit width as presented 
to the beam can be varied from 0 to 0.152 cm. The usual width for the present 
experiments was 0.0076 cm. 

The detector is a tungsten wire 0.0127 cm. in diameter with a slit in front 
of it of 0.010 cm. width. The ions that are formed by surface ionization are 
pulled off at right angles to the incident beam to a mass spectrometer. The 
ion optical system for this extraction is identical with that used by Braun- 
stein and Trischka (1955). With this ion gun, although 80% of Rb ions formed 
on the tungsten wire could be pulled out of the gun, only 20% reached the 
collector of the mass spectrometer. The usual energy of the ions was about 
1000 ev. No thorough attempt has been made to determine the exact cause 
for the loss of ions. 

At the collector end of the mass spectrometer the ion current may be col- 
lected directly in a cage and measured on an electrometer circuit or it may be 
amplified by means of an electron multiplier and then measured on an electro- 
meter circuit. The latter arrangement was the one usually used in the experi- 
ments. 

The obstacle wire in the middle of the B field is a nylon filament of 0.051 
cm. diameter. It is coated with colloidal graphite and maintained at a potential 
midway between those of the B field pole pieces. 


IV. ELECTRONIC EQUIPMENT 


In order to measure transitions to an accuracy of 1 part in 10° the micro- 
wave signal must be stable to about 1 part in 10’ and its frequency purity 
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should be no worse than | in 10°. A system fulfilling these requirements is 
shown in Fig. 6. Although this system is essentially the same as that used by 
Lee et al. (1953) and later by Stroke, Jaccarino, Edmonds, and Weiss (1957), 
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Fic. 6. The microwave system. 


it is considered desirable to describe it here both for the sake of completeness 
and for the purpose of mentioning some details which may be of interest to 
anyone wishing to set up a similar system. The source of microwaves is a 
Varian X-12 klystron immersed in oil. For greater ease of understanding 
the operation of the system the klystron is shown generating a definite fre- 
quency of 12574.214 Mc./sec., which is actually the frequency of one of the 
lines observed. Part of the energy is fed to the C field through a series of 
attenuators, as shown. Part of the energy is fed into a hybrid tee where 
mixing with two reference frequencies takes place. The reference frequencies 
are obtained by crystal rectification of 270 Mc./sec. To obtain adequate 
power in the region of the 47th harmonic to permit operation of the negative 
feedback loop, a crystal current of about 10 ma. is necessary. The 270 Mc./sec. 
standard is obtained by mutiplications of a 1 Mc./sec. signal from a Primary 
Frequency Standard maintained by the Division of Applied Physics of this 
institution. The multiplier circuit is that given by Hedrick (1953). This 
circuit has provision for introducing sidebands into the output at 5, 10, 30, 
or 90 Mc./sec. from the 270 Mc./sec. Hence at the output of the crystal 
rectifier one can find the frequencies 12,600 Mc./sec. and 12,610 Mc./sec. 
With the former frequency the klystron signal produces a beat of 25.786 
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Mc./sec. and with the latter a beat of 35.786 Mc./sec. The 25.786 Mc./sec. 
signal is used simply for frequency measurement. That is, this beat frequency 
is measured precisely to yield the klystron frequency precisely. The 35.786 
Mc./sec. signal is used for frequency control of the klystron. The FM receiver 
is tuned to this frequency and if the beat tends to drift from this frequency 
as a result of a drift in the klystron frequency an error voltage is produced 
which can be applied to the klystron reflector through the reflector control 
circuit to counteract the drift. The reflector control circuit is given in detail 
in Figure 8 of Stroke et a/. (1957). With this circuit and a Hallicrafter SX-62 
FM receiver it was found necessary to connect a 0.5 uf. condenser across the 
output of the FM discriminator in order to prevent the feedback system from 
oscillating. 

The feedback system not only stabilizes the klystron frequency but also 
reduces its frequency spread. Without the feedback, the output of the Varian 
X-12 has a frequency spread of about 100 kc./sec.; with the feedback the 
spread is reduced to 15 kc./sec. The 100 kce./sec. spread cannot be ascribed 
to ripple in the power supply voltages; the heater of the klystron is operated 
from a storage battery and the 3 millivolt ripple in the beam and reflector 
voltages of the FXR Model Z 815B power supply is not sufficient to produce 
this spread. Furthermore the fundamental frequency of this noise is not 60 
c./sec. but is variable between 65 and 70 c./sec. 

It is difficult to determine just how stable a signal is which has a noise 
bandwidth of 15 kc./sec., but repeated measurements of a given signal have 
been found to yield results which agree to +1 ke./sec. Thus the stability is 
at least 1 part in 10’. 

The 100 ke./sec. Primary Frequency Standard from which our 1 Mc./sec. 
standard is derived is always within 3 parts in 10’ of the nominal frequency 
and is known to better than 1 in 10°. 


V. EXPERIMENT 

Because RbF is highly deliquescent it was found necessary to dry it in 
an auxiliary vacuum system at 800°-900° K. before loading it in the oven. 
Even with this preliminary drying the substance was found to have a strong 
tendency to clog the slit, especially after the oven had been raised to the usual 
operating temperature of about 1030° K. The only way the writer found of 
avoiding the clogging was to load the oven very lightly. The usual charge 
was about 0.2 g. This was sufficient for several days’ operation. 

The beam intensities used in the measurements were about +X 10~-"° amp. 
or 400,000 cm. of galvanometer deflection (1 cm. = 10-'® amp.). This was 
the intensity as measured at the tungsten wire detector and not at the end 
of the mass spectrometer. This intensity is considerably less than those used 
by others in similar experiments (cf. Lee et a/. 1953; Fabricand et al. 1953; 
Trischka and Braunstein 1954). Attempts to increase the beam intensity 
were not very successful. Furthermore the attempts suggested that a thorough 
investigation into the reasons for the low intensities might be a long drawn- 
out affair. Hence, since transitions were already observable with the available 
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intensities, it was decided to proceed with the measurements and see how far 
one could get. It turned out that one could get most of the information desired 
but that the measurements would have been easier and the accuracy greater 
if greater intensities had been available. 

Since the operating conditions of the tungsten wire have a direct bearing 
on the detected beam intensities, it is considered desirable to record them here. 
The tungsten wire was of commercial quality* 0.0127 cm. in diameter. Before 
being used for beam detection it was heated for several hours at around 
890° C. in its normal mount in vacuum with a current of 1.0 amp. The above 
temperature and the 1500° C. mentioned below are temperatures as measured 
directly by means of an optical pyrometer and not corrected for emissivity. 
This gentle baking out removed most of the alkali impurities in the tungsten, 
as indicated by the ion current coming out of it, but did not remove more 
refractory impurities such as oxides. Then just prior to being used for beam 
detection the wire was flashed for a few seconds at 1500° C. with a current 
of 1.7 amp. and immediately returned to 890° C. at 1.0 amp. Immediately 
after flashing, the wire often showed a very high sensitivity, i.e. one in which 
over 1.2 million cm. of RbF was detected. The sensitivity usually dropped 
very rapidly after flashing and levelled off in about 20 minutes to a value 
roughly one-third the initial value. It was at this lower sensitivity that the 
experiment was done.f The flashing had to be repeated occasionally. 

In the measurement of a line it was found difficult to plot a line shape 
because of fluctuations in the refocused beam. There were quiescent periods 
between fluctuations, however, which were of sufficient duration to permit 
a reasonably reliable determination of the peak of the transition. Hence all 
the results of the present work were obtained by repeated settings of the micro- 
wave oscillator on the peak of the line. Under the best of conditions and for 
the strongest lines, the signal-to-noise ratio was about 20 to 1. The line widths 
were estimated to be approximately 50 ke. /sec. On the basis of the uncertainty 
principle the expected line width is 10 ke./sec. The large observed width 
probably arose partly from the inhomogeneity of the microwave signal and 
partly from the fact that perhaps more than the optimum power was required 
for the transitions. 

VI. RESULTS AND DISCUSSION 
1. Observed Lines, Quadrupole Interaction 

In Rb*F, rotational transitions belonging to three vibrational states 
v = 0, 1, 2 were observed, whereas in Rb*’F, lines belonging to the v = 0 
state only were observed. These are all listed in Table II, without any cor- 
rection for the contact potential of the C field. It is shown below that the 
correction amounts to only 400 c./sec. or well below the uncertainties of the 


measurements. 


*“Elmet’’ brand by North American Philips Co., Lewiston, Maine. 

tExperiments done after the measurements suggest that a thorough aging of the tungsten 
followed by oxidation may yield a detecting surface of constantly high sensitivity (Datz 
and Taylor 1956). 
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TABLE II 
OBSERVED FREQUENCIES (Mc./SEC.) 
(JJ =1)<—(VJ =0) 


F, of J =1 
Line level v=0 v=1 y=2 

Rb*®F 

( 3/2 12595.316+0.005 12504. 307+0.010 

B 7/2 12588 .985+0.005 12498 .054+0.010 

\ 5/2 12574.214+0.005 12483 .444+0.010 12393 .058+0.010 
Rb*F 

Cc 1/2 12541.181+0.010 

B 5/2 12534.381+0.010 

\ 


3/2 12525 .870+0.010 


The quadrupole interaction constants in those cases in which all three of 
the hyperfine structure lines have been measured are given in Table III. 
For comparison the corresponding values found by Hughes and Grabner 
(1950) are also given. It is seen that these agree very well with our somewhat 
more accurate values. 


TABLE III 
VALUES OF —egQ/h FOR THE ROTATIONAL STATE J = 1 IN Mc./SEC. 








Molecule Vibrational state Present work Hughes and Grabner 
Rb*®F 0 70.340+0.0380 70.31+40.10 
l 69 .548+0.060 69.54+0.10 


Rb*F 0 34.024+0.040 34.00+0.06 


eqQ® /eqQ® = 2.0673 0.0030. 


From the ratio of the quadrupole interactions in the zeroth vibrational 
state in the two isotopic species we find for the ratio of the nuclear quadrupole 
moments of the Rb isotopes 


Q*/Q% = 2.0673-+0.0030. 


This agrees with the more accurate value of 2.0669+0.0005 found by Trischka 
and Braunstein (1954) in their study of RbCl. 


B. Rotational Constants 

The positions of the rotational levels in the absence of hyperfine structure 
are easily determined in accordance with the patterns shown in Figs. 1 and 2. 
This can be done for the cases v = 0 and 1 in Rb*®F and v = 0 in Rb*’F. For 
the vibrational state v = 2 of Rb*®F, however, since only one of the three 
hyperfine levels has been measured and the quadrupole interaction constant 
is not determinable, it is necessary to use Hughes and Grabner’s value for 
this interaction constant, 68.71+0.10 Mc./sec. The rotational transition 
frequencies unperturbed by h.f.s. are then as shown in Table IV. 





LEW ET AL.: RbF ROTATIONAL SPECTRA 181 


TABLE IV 
CALCULATED ROTATIONAL TRANSITION FREQUENCIES 
IN THE ABSENCE OF HYPERFINE STRUCTURE (Mc. /SEC.) 
(J =1)—(VJ =0) 


0 12585 .468+0.010 


Rb®F v= 
v= 1 12494 .571+0.020 
v=2 12404 .052+0.026 
Rb*’'F v= 


0 12532 .674+0.020 





The various rotational constants calculated from the above frequencies 
are shown in Table V. The calculations required a knowledge of Yio = w, 
and Ya = —w,x,. From a study of the ultraviolet absorption spectra of 


TABLE V 
ROTATIONAL CONSTANTS (MC./SEC.) 


Rb®F Rb*"F 


Yn = B. 6315 .543+0.009 6289 .002+0.010 
— Yn = ae 45.6388+0.017 45 .350+0.017 
Yn = Ye 0.094+0.007 0.09440. 007 
Be 6315.689+0.036 6289 .098+0. 036 





alkali halides, Barrow and Caunt (1953) give for RbF, w, = 390 cm.—! and 
w-X. = 1.3 cm.—!. Hence we find Yo. = —4B,3/w?2 = —0.007 Mc./sec. and, 
by formulae conveniently listed in Lee et a/. (1953), 


B, = Yut(0.096+0.035) Mc./sec. 


The probable error in the added term above is based on an estimated 5% 
probable error in the spectroscopic value of w, and w,x,. Thus the quoted 
errors in B, arise almost entirely from the errors in the vibrational constants. 

In the case of Rb*’F, where only one vibrational state has been measured, 
it was necessary to assume a,” = (us5/us7)??a,. and Yo? = Yo, where 
the »’s are the reduced masses of the molecules. 


C. Mass Ratio 
The ratio of the masses of the Rb isotopes may be calculated from the rela- 
tion (Fabricand et al. 1953) 


Mss _ (Yor / Yor”) (Mr /Msz) 


Mex 1+(My/Ms7)—(Yo:/Yor”)’ 
where Vy is the mass of the fluorine atom. Strictly speaking, in the above 
formula, the ratio of the Yo:’s should be replaced by the ratio of B,’s. However, 
these two ratios agree to within 1 part in 10’ or well within the accuracy of 
our experiment. We find 


Mss/Msz = 0.9770148+0.0000052. 


For comparison, this and other measurements of the same ratio are shown 
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in Table VI. It seems the spectroscopic values are consistently lower than the 
mass spectroscopic value even though all the values overlap within the stated 


errors. 
TABLE VI 
MEASUREMENTS OF THE MASS RATIO OF THE Rb ISOTOPES 
Method Molecule Ratio (Mg5/ M37) Ref.* 

Mass spectroscopy - 0.9770191+0 0000022 (a) 
Microwave absorption RbI 0.9770177 +0 .0000045 (b) 

RbBr 0.9770146+0 . 0000055 (b) 
Molecular beams RbCl 0.9770168 +0 .0000045 (c) 


RbF 0.9770148+0 . 0000052 (d) 


*(a) Collins, Nier, and Johnson (1951). (6) Honig et a/. (1954). (c) Trischka and Braun- 
stein (1954). (d) This work. 


D. Internuclear Distance 
The equilibrium internuclear distance of RbF calculated from B, is 


re = (2.26554+0.00005) X LO-§ cm., 


where the uncertainty arises entirely from those in Planck’s constant and 
the atomic mass unit (Cohen, Dumond, Layton, and Rollet 1955). This 
internuclear distance is somewhat larger than either of the values 2.225 and 
2.246 given by Honig eft al. (1954) on the basis of the simple additivity of 
ionic radii. However, the agreement is within the 1 or 2% accuracy claimed 
for the latter. 

Honig et al. also give an empirical formula for internuclear distances in- 
volving the polarizabilities of the ions concerned and the electronegativities 
of the atoms. Using polarizabilities given by Rittner (1951) and electronega- 
tivities given by Pauling (1940) we obtain a value of 2.2545 A. This is only 
1% lower than the experimental value. This agreement suggests that the 
empirical formula is a good one to use to predict the B, values of alkali halides 
which have not yet been measured. 


E.. Electric Dipole Moment 

For the determination of the electric dipole moment of RbF in the v = 0 
state the Stark effect of the 12595.316 Mc./sec. line was studied. This line 
was chosen for reasons given in Section II. Static voltages of up to 20 volts 
were applied to the C field in which the gap was 0.69065 cm. The shift in the 
frequency of the line was measured with the static voltage first in one direction 
and then in the other. Typical shifts are, for V. = +20.00 volts, Af = (0.461 
+0.006) Mc./sec. and for V. = —20.00 volts, Af = (0.407+0.010) Mc./sec. 
The difference in the frequency shifts is of course due to the contact potential 
of the C field. From the measurements it is calculated that the contact potential 
of the aluminized surface coated with silicon dioxide relative to the brass 
surface is +(0.61+0.05) volts. The electric dipole moment of Rb“F in the 
v = 0 state is then found to be, by Eq. 2, 


nu = (8.80+0.10) X 10-"8 e.s.u. 
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Rittner gives an expression for the dipole moment of an alkali halide in- 
volving the polarizabilities of the ions. This relation yields a value of 7.96 X 10~"* 
e.s.u. The calculated value is thus lower than the experimental value, a 
relationship which is consistent with what is observed in all alkali halides for 
which the dipole moments are known. 

It is easily seen that the contact potential of 0.61 volt introduces a negli- 
gible error in the measurements of the lines at “zero” C field. The energies 
of the Stark levels of a molecule in the 'Y state do not depend, to first order, 
on the electric field when this field is small. Among the levels with which we 
are concerned the one with the highest second-order dependence on the external 
field is the J = 0 level. At a field of 0.883 volt/cm. corresponding to a potential 
of 0.61 volt, the shift in the J = 0 level of RbF is only 0.0004 Mc./sec., 
which is well below the limit of detection in the present experiment. 
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MASS SPECTROMETRIC ANALYSIS OF SOME 
HYDROGEN OXIDES. [ 


LARKIN KERWIN AND MAURICE COTTIN? 


ABSTRACT 


Typical mass spectrograms for H2O, D2O, H202, and D202 are presented. Both 
positive and negative ions were observed. The variations of the abundances of the 
spectrogram components with pressure are given. Some features of the experi- 
mental data are discussed. 


1. INTRODUCTION 

The mass spectrometric investigation of the hydrogen oxides presents 
considerable experimental difficulty. Water has been studied by several 
workers, including Smyth and Mueller (1933), and Mann, Hustrulid, and 
Tate (1940) in which reference is made to other papers. Hydrogen peroxide 
has recently been studied by Robertson (1952). No systematic investigations 
of heavy water and deuterium peroxide have appeared in the literature, 
although Washburn, Berry, and Hall (1953) have made a careful study of 
the precautions to be taken when measuring the small quantities of D,O 
present in natural water. Certain hydrogen—oxygen ions and free radicals have 
been the subject of close scrutiny, and papers by Foner and Hudson (1956) 
and the works cited therein are typical of results on such species as HO: and 
OH. Negative ions have received much less attention than positive ions. 

We have begun a mass spectrometric study of water, heavy water, hydrogen 
peroxide, and deuterium peroxide. In this paper, we report the mass spectro- 
grams obtained and the variation of the ion abundances as a function of 
pressure. The second part of the investigation dealing with appearance poten- 
tials and dissociation energies will be reported later. 


2. EXPERIMENTAL 

It is usual to improve the ion yield of a mass spectrometer source by causing 
the ionizing electrons to move in spiral paths, which effect is obtained by 
placing the source in a magnetic field. When the spectrometer fields are 
reversed for analyzing negative ions, the source field must be reversed as 
well, and even this does not usually eliminate the need for refocusing the 
ion beam. In addition, reversal of the source field inevitably affects the 
efficiency of the source, making comparisons difficult. In modifying our 
inflection spectrometer for studying negative ions the source magnets were 
removed. The inherent high intensity of the instrument (Kerwin 1950) 
resulted in the beam intensity still being sufficient to yield useful informa- 
tion. It was not necessary to reset the focusing controls, and positive and 
~ 1Manuscript received October 1, 1957. 
Contribution from the Département de Physique, Université Laval, Québec, Que. 

2On leave from L’Institut du Radium, Paris. 
Can. J. Phys. Vol. 36 (1958) 
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negative ions could be analyzed in rapid succession while source and sample 
characteristics remained constant. 

The nature of the samples analyzed required meticulous avoidance of 
impurities in the pyrex sampling system. This was kept as short as possible. 
After repeated cleaning with fuming H2SO,, the glass was heated to the 
softening point. This last step was essential in preventing decomposition of 
the peroxides. The samples were introduced directly from the pyrex line into 
the ionizing electron beam. Two stopcocks in the line were greased with Dow 
Corning silicone grease. This, and the conditioning of the spectrometer source 
by flushing with deuterium oxide vapor before runs, minimized deuterium 
exchange with hydrogen from adsorbed water in the sample line during runs. 
These precautions are similar to those recommended by Washburn et a!. 

All of the samples were in liquid form. Vapor pressure was controlled in 
part by varying the sample line leak diameter, and in part by varying the 
temperature of the liquid samples. These were immersed in dry ice — acetone 
baths, melting ice, or mixtures of the two. Sometimes reproducible runs could 
be obtained at intermediate pressures by allowing the sample to warm slowly 
from one type of bath to another. The pressures reported were measured by 
several gauges, and represent the pressure in the ionizing region to within a 
probable factor of two. 

The spectrometer was usually operated with a 100 ev. ionizing current of 
10 va. Positive and negative ions were detected with a vibrating reed electro- 
meter. At 10-‘ mm. Hg, the main spectrum components provided ion currents 
of about 10~'° amp., the smaller components currents of about 10-" amp. 

For a typical analysis, several runs were made on each of several successive 
days. Many such runs on the peroxides were rejected because slight decom- 
position from day to day altered the spectrogram significantly. The sample 
line would then be cleaned as described until the spectrogram could be repeated 
over a several-day period. The results presented here are composites of two 
series made at 3-month intervals, each series being composed of several runs 
considered acceptable as described. The spectra are considered typical of 
ordinary instrumental variations, all acceptable runs made over a period of 
months fitting them to within about 10%. 


3. MASS SPECTROGRAMS 

The first four figures show typical mass spectrograms of the species investi- 
gated. Their liquid purity was: water (distilled), heavy water (about 80% D), 
hydrogen peroxide (99.2% peroxide), deuterium peroxide (over 99% peroxide, 
and 99% D). These samples were graciously furnished by Professor Paul-A. 
Giguére. 

The percentage abundance is referred to the total of all ion currents ob- 
served. The shaded portion of the abundance bars represents the variation in 
abundance with pressure from 10~! mm. Hg to 10-* mm. Hg. Therefore typi- 
cal spectrograms obtained within this pressure range should fall within 
the shaded portion. The corresponding abundances for water reported by 
Mann, Hustrulid, and Tate are indicated. The expansion of the negative 
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ion abundance scale is to be noted. With each figure are listed the ions con- 
sidered as corresponding to the various masses, in order of increasing mass. 
Except where noted, molecules containing O-18 were not considered as con- 
tributing significantly to peaks. 
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Fic. 1. Mass spectrogram of water. Ions considered observed: H+, H2*, O+, OH*, H20*, 
H;0+, and H,O* (isotopic O-18); H~, O-, and OH-. 

Fic. 2. Mass spectrogram of heavy water. Ions considered observed: H*+, H2++D*, DH*, 
D,t, O+, OH+, H2O++DO+, DHOt, D.0*+, DzHOt, and D;0*; H-, D-, and O-. 

Fic. 3. Mass spectrogram of hydrogen peroxide. Ions considered observed: H*, H.*, O*, 
OH+, H,O+, H;0+, H2O* (isotopic O-18), O2+, HO.*, and H:0.+; O-, and OH-. 

Fic. 4. Mass spectrogram of deuterium peroxide. Ions considered observed: H*, He++D*, 
O+, OH+, H.O++OD+, DHO*, D.O+, DzHO*, D;0*, O.+, HO2*, DO.*, DHO2*, and D202*; 
O-, and DO-. 


4. ABUNDANCE VARIATIONS WITH PRESSURE 
Figs. 5 and 6 show the variation of the major ion abundances as a function 
of pressure. In Fig. 5 are grouped the light and heavy waters, in Fig. 6, the 
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light and heavy peroxides. Thus in Fig. 5, the curve labelled H* also repre- 
sents Dt, that labelled H;O0+ also represents D;0+, D,HO+, and DH,O+, and 
so on for the other species. 

The calculations involved for separating the effects of these species are 
simple but lengthy. As an example, from a spectrogram of heavy water such 
as is given in Fig. 2, peaks 21 and 22 may be ascribed to D,HO* and D,O*, 
since no other hydrogen—oxygen ions may be considered at these masses. We 
let the deuterium abundance equal d, and the hydrogen abundance (100—d). 
Now the D,HO?t ion may be present as DDHO, DHDO, or HDDO. Accord- 
ingly its abundance will be proportional to 3d?(100—d). The abundance of 
the D,0+* ion will be proportional to d*. We therefore have: 

3d°(100—d) _ abundance mass 21 








a abundance mass 22 


From this equation d may be calculated to somewhat low precision, e.g. 
about 80%. However, the figures obtained serve to establish that peak 19 
must be caused effectively by DHOt, since H;0*+ makes a negligible con- 
tribution, while peak 20 is caused by D,O+*, the form DH,O* being eliminated 
for the same reason. The value of d may then be established with more precision 
from these two peaks by the appropriate equation. For the 10-‘mm. Hg 
values illustrated in Fig. 2, d = 79.9%. Finally, the contributions of the 
various ions to each mass peak may be calculated. 

In this way it was shown that in general no isotopic effect was observed in 
the abundance-pressure relationships. Some particular effects that seem 
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Fic. 5. Ion abundance as a function of pressure, water and heavy water. 
Fic. 6. lon abundance as a function of pressure, hydrogen peroxide and deuterium peroxide. 
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established will be mentioned later. Except for these, the curves of Figs. 5 
and 6 represent the behavior of the isotopic variations of a given ion. The 
sum of the component abundances at a given pressure does not always equal 
exactly 100% because the curves are averages of many runs. The reproduction 
of the curves for the waters is somewhat better than that for the peroxides. 


5. DISCUSSION 
(a) Water 

The type and abundance of the ions observed (Fig. 1) are essentially the 
same as those found by Mann et a/., with the exception of the low-abundance 
OH-~- which they did not observe. Massey (1950) reported that Sloan and 
Love have observed OH™~ in the residual gas background of their spectro- 
meter. Smyth and Mueller also observed this ion. We found no HOz*, even 
when oxygen was mixed with water vapor in the source, although an abundance 
of 0.01% would have been detected. 

At low pressures the decomposition of water into OH* and Ht? is relatively 
greater than at high pressures (Fig. 5) and it was verified that this was not 
a result of any variations in the ionizing current. The hydroxonium ion H;O+ 
increases in abundance as the square of the pressure (thus its percentage 
abundance increases linearly) as was observed by Smyth and Mueller, and 
Mann et al. In a few experiments at very high pressure (filament life was 
drastically curtailed) it was possible to obtain hydroxonium abundances as 
high as 100% of the parent H,O+ peak. This supports the view of Lefort 
(1954) that H,0* would be one of the most abundant ions in radiolized liquid 
water. The variation of the ion abundance with the square of the pressure 
indicates formation by secondary reaction (Smyth and Mueller): 


H.O+ 4 H,O — H;0+ + OH. 


(b) Heavy Water 

The D,O spectrogram (Fig. 2) indicated essentially the same species of 
ions as were formed from water. When the various isotopic versions are 
separated out as described in Section 4, the calculated mass spectrogram of 
pure heavy water may be obtained, and is shown in Fig. 7. 

This is seen to be similar to Fig. 1 for HO, with the exception of OD~ at 
mass 18, whose abundance was less than the limit of sensitivity of the apparatus, 
as may be calculated from the D content of the sample. 

The variation in ion abundance with pressure was essentially the same as 
for water (Fig. 5) with the exception of the following small effects. The 
D,O* concentration in the ‘pure’ heavy water spectrogram was about one-half 
that of the H,O* concentration in the ordinary water spectrogram. This was 
true over the pressure range described. A second systematic effect indicated 
that the D.O* ion abundance was about 5% lower than the corresponding 
H.O* abundance, with the DO* abundance being correspondingly higher than 
the HO* abundance. These effects are indicated by dotted lines in Fig. 5. 


(c) Hydrogen Peroxide 
The higher-purity peroxide furnished by Professor Giguére gave about 96% 
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peroxide in the vapor phase. The relative abundance of the ions (Fig. 3) 
agrees roughly with those reported by Robertson, who did not give his com- 
plete spectrum in his paper, so that a peak-by-peak comparison is not possible. 


Calculated spectrogram 
7} Pure 0,0 
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Fic. 7. Calculated mass spectrogram, pure D.O. 


Our spectrum for the O.+, HO.*+, H,O.* group (32, 33, 34) corresponds to his 
but is about seven times more abundant compared with the H.O* peak, 
indicating less spontaneous (non-electron-impact) dissociation in our instru- 
ment. Robertson states that he obtained much dissociation and oxidation of 
metal in his source. 

The HO ,* ion appears in relatively high concentration— about 10°% of the 
H,O,.* ion. Foner and Hudson reported the ion as being 1°% of the peroxide 
ion abundance when using a discharge source with water vapor mixed with 
the peroxide. 

No negative hydrogen ions were found, doubtless because of lack of sensi- 
tivity. The OH™> ion was about twice as high as the O~ ion, while in the 
H,O-D.0 analyses, OH~ was only about 1° of the O~ peak. This supports 
the generally-held view that the peroxide is most prone to division into 
2(OH) as the primary step in its dissociation. 

The pressure-abundance curves (Fig. 6) are essentially linear, except for 
the decrease in H,O+ and OH* abundance with pressure. The H,O0* ion 
increases with the square of the pressure as before. 


(d) Deuterium Peroxide 

The ions considered observed are listed under Fig. 4, from which the very 
low hydrogen content of the sample is apparent. Noteworthy are the peaks 
D.HO+ and D,O+ at 21 and 22, and the HO,* and DO,* peaks at 33 and 
34. When the various ion species are separated out of Fig. 4 as indicated in 
Section 4, the calculated spectrogram for ‘pure’ deuterium peroxide may be 
obtained, and this is shown in Fig. 8. 
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Once again, no H~ ion was observed, because of lack of sensitivity. OD- 
was as relatively abundant as the corresponding OH™ in the peroxide spec- 
trum. 


Calculated spectrogram 
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Fic. 8. Calculated mass spectrogram, pure D2Qx. 


The variation in ion abundance with pressure was essentially the same as 
for hydrogen peroxide. The isotope effect on hydroxonium formation was 
noted, but the effect on the D,.O+ and DO?t peaks was absent, to within the 
precision of the experiment. 


6. CONCENTRATION OF O-18 

The D.O, used in these experiments is manufactured in an electrodeless 
discharge (Giguére, Secco, and Eaton 1953). The microwave spectrum of 
this heavy peroxide produced a great array of lines (Massey 1952; see also 
Massey 1954) and it has been suggested that these might have been caused 
by a concentration of the oxygen-18 isotope in the discharge. A careful 
examination of the mass spectrograms indicates that the concentration, if 
present, does not exceed a factor of three, i.e. the abundance of O-18 would 
be less than 0.7% of the O-16 isotope abundance. 
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ON VOLTERRA’S DISLOCATIONS IN A 
SEMI-INFINITE ELASTIC MEDIUM! 


J. A. STEKETEE 


ABSTRACT 


In this paper a Green's function method is developed to deal with the problem 
of a Volterra dislocation in a semi-infinite elastic medium in such a way that the 
boundary surface of the medium remains free from stresses. (A Volterra disloca- 
tion is here defined as a surface across which the displacement components show 
a discontinuity of the type Au = U+ Q Xr, where U and Q are constant vectors.) 
It is found that the general problem requires the construction of six sets of Green's 
functions. The method for the construction is outlined and applied to one of the 
six sets, which is of the type of two double forces with moments in a plane parallel 
with the boundary. The displacement field thus generated is computed. Several of 
the results obtained are believed to be of geophysical interest, but a more detailed 
discussion of these applications is postponed to a further communication which 
is being prepared. 


1. INTRODUCTION 

In the orogenetic considerations of J. T. Wilson and A. E. Scheidegger an 
important part is played by island arcs and arcuate structures in mountain 
ranges. It is supposed that these arcuate structures are the surface representa- 
tions of fracture zones in the orogenetic shell. These fracture zones, along which 
the foci of earthquakes are to be found, have been described on different 
occasions by the above authors (Wilson 1954; Scheidegger 1953) and others 
(compare Umbgrove 1954). It has been suggested by the author (Steketee 
1955) that Volterra’s theory of dislocations might be the proper tool for a more 
quantitative description of these fracture zones.? This suggestion was not 
made without reason. If the mechanism involved in earthquakes and the 
fracture zones is indeed one of fracture, there will exist discontinuities in the 
displacement components across the fractured surfaces. As dislocation theory 
may be described as that part of the theory of elasticity dealing with surfaces 
across which the displacement components are discontinuous, the suggestion 
seems reasonable. 

Another argument is that in the last 20 years dislocation theory has become 
of major importance in solid state physics, and it should therefore not be ex- 
cluded that also in geophysics it may be of some importance. Finally it was 
realized that the mathematical structure of dislocation theory is analogous to 
certain parts of potential theory in the sense that the first theory deals with the 
Navier equations and the second with the Laplace equation. In view of the 
many applications of potential theory in physics it was again felt that disloca- 


tion theory might be of some use in geophysics. 

‘Manuscript rece ‘ed October 21. 1957. 

Contribution from the Department of Applied Mathematics, University of Toronto, 
Pcronto 5, Ontario 

*Similar suggestions have been made at other places and similar times; the referee kindly 
informed me that G. W. Housner applied dislocations to earthquake faulting in 1953 (Housnet 
1955), while for example Dr. Kevlis-Borok in Moscow and Dr. Knopoff in California were 
probably working on similar problems at the time 
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Rochester (1956) was the first one to obtain practical results with this theory 
when he worked out an example referring to the San Andreas fault in Cali- 
fornia. His results seemed to justify further work along these lines. 

As usual in mathematical physics, it is necessary for simplicity’s sake to 
omit certain features from reality and to retain only those which are expected 
to be of principal importance. We omit therefore the curvature of the earth, 
its gravity, temperature, magnetism, and non-homogeneity, and shall deal 
with a semi-infinite medium, homogeneous and isotropic, assuming further the 
laws of the classical linear elasticity theory. 

The usual boundary conditions for dealing with the earth's problems require 
that the surface S of the elastic medium (the earth) shall be free from forces. 

To solve this boundary-value problem for the case that a dislocation is 
present in the medium, Volterra (1907) has given the following method: First 
compute the stresses and displacements in an infinite elastic medium due to the 
dislocation, in particular compute the forces on S. Then consider the ordinary 
boundary-value problem for the elastic body with surface S, where the surface 
is subjected to forces equal and opposite to the ones obtained before. Super- 
position of the two solutions will give the required result with a dislocation in 
the elastic body and no normal or tangential forces on its surface. 

It isclear that with this method for each new type of dislocation a new boun- 
dary-value problem has to be considered. Although the solution of this 
boundary-value problem does not offer fundamental difficulties in the case of 
the semi-infinite medium, the computations one has to face are complex so that 
for example in Rochester’s work the ordinary boundary-value problem could 
only be solved approximately. 

It will be the purpose of the following investigation to develop a Green's 
function method for the calculation of dislocations in a semi-infinite medium. 
The construction of these Green's functions is simple compared with the 
problems faced when following the method of Volterra. Once these functions 
are obtained, the problem of an arbitrary dislocation in the semi-infinite 
medium is reduced to the computation of certain surface integrals. If approxi- 
mations have to be made, it can be done in these integrals in the final stage 
of the problem. 

While the work contained in this paper was undertaken with the geophysical 
problem in mind, we shall be largely concerned here with the elasticity problem. 
It is expected to give some geophysical applications in a later communication. 

In Sections 2 and 8 we deal with the theory of dislocations in an infinite 
elastic medium. In Section 4 the boundary-value problem is discussed and the 
method for the construction of the Green’s functions is outlined. In Sections 5 
and 6 some parts of the construction are dealt with in greater detail. Finally 
the method is applied in Section 7 to one of the six groups of Green's functions. 

2. VOLTERRA’S THEORY OF DISLOCATIONS 

Let O be the origin of a Cartesian frame of reference in an infinite elastic 
medium, xv, the Cartesian coordinates (¢ = 1, 2,3), and e; a unit vector in the 
positive v,-direction. A force F = Fe, at O, compensated by continuously dis- 
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tributed forces at infinity to establish equilibrium, generates a displacement 
field u*(P, O) which is determined by the well-known Somigliana tensor 


. 
4 


(2:1) ui (P, O) = es (6 x —_ ar ix). 
Sau 


In this relation 6, is the Kronecker delta, a = (A+yu)/(A+2u), \ and uw are 
Lamé’s constants, 7 is the distance from P(x,) to O. Further 7, = dr/dx; and 
the summation convention applies. If the force is located in “Oly,) we have 
r? = (x;—y,)? and dr/dx; = —dr/dy;. It will be convenient to assume the 
constant F/82p equal to one. The dimension of this constant is [L?], so that 
the formula (2.1) gives indeed a displacement with the proper dimension. 

The stresses due to the displacement field (2.1) are easily computed from the 
general relation 


(2:3) Tis = NOG Uy tu (Ui, j+U;, i). 
We find 


A 


. Xx , Xe 
(2.3) Tis(P, O) = 2u(1—a)( 6, a a8 € ted ; 5 é 
r r 
Further we denote the components of the force per unit area on a surface 
element, following Sokolnikoff and Specht (1946), as 


Kk 
(2.4) qT = Thy. 0; 


where the v,’s are the direction cosines of the normal to the surface element. 
A Volterra dislocation i is defined as a surface 2 in the elastic medium across 


which there is a discontinuity Au, in the displacement components of the type 


) Au, =u; —u;, = Ui+Qx,, 


(3:5) 
\ Q,, = —Q),. 


The formula (2.5) in which U; and Q,;, are constants is the well-known Wein- 
garten relation which states that the discontinuity Au, should be of the type 
of a rigid body displacement, thereby maintaining continuity of the compo- 
nents of stress and strain across L. It is also clear that the edge of J will bea 
singularity as the displacement is essentially indeterminate. 

The displacement field in an infinite elastic medium due to a dislocation of 
type (2.5) is then determined by the relation of Volterra (1907) (compare 


Bergman and Schiffer 1953, pp. 219-222): 


, l a 
(2.6) Ux(y.) = Au, 7, dd, 
Sue e 
z 
Pi 
where 7°; is the 7-component of the force per unit area on a surface element 
d> containing the point P(x,) due to a force 8mpe, in Q(y,). The dimension of 
the constant in front of the integral is [M L~! T~?] so that the formula gives 
the proper dimension for u,. 
Once the surface 2 is given, the dislocation is essentially determined by the 
six constants U’; and Q;,;. We write therefore also: 
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97 ) 5 U; ( k P d> 
(2.4) Ux(y¥r = Sau ; 7i3( ’ Q) vida 
z= 
Qi; ' " af i ey v 
—— (xyri(P, Q)—xit(P, Q)}v, dd, 


where Q,,; takes only the values Q)2, Qe3, Q3:. Following Volterra (1907) and 
Love (1944) we call each of the six integrals in (2.7) an elementary dislocation 
(‘‘Distorsion élémentaire’’). 

It is clear from (2.6) and (2.7) that the computation of the displacement 
component u,(Q) is performed as follows: A force 8xye, is put in Q, and the 
stresses rx ,(P, Q) which this force generates are cena in the points P(x,) 
on Y. In waaticalen the components of the force on = are computed. After multi- 
plication with prescribed weights of magnitude Au; these forces are integrated 
over ~ to give the displacement component in Q due to the dislocation on 2. 
We shall later see that a different interpretation is possible which will turn 


out to be more practical. 


3. FURTHER DEVELOPMENTS, SYMMETRY PROPERTIES 
Once the displacement field u,(Q) has been obtained, the stresses can be 
found by differentiating to y, according to (2.2). If the point Q(y,) is not on &, 
these differentiations can be performed in general under the integration sign 


where the variables x,;— y, appear in 74,;. We obtain then for points not on ¥ 
: ) U, CC px: — kl Gt - 
(3.1) T%(Q) = Sau Giyv,d> = as Gin —X iG jn) ¥nd = 
where 
(3.2) Gi(P. Q) = Net tury +745) 
and 

«8 k 

OT Or k 
(3.3) nie = — 2 = — gs, 

‘ Oy, Ox, ag 


Writing out the expression (3.2) for the case that Q is the origin, we obtain 


Bide: du (1—a) 
r> r 


(3.4) —Gii(P,O) = 4y?(2—3a) (55: +5510} 


\ re 
— 12y*(1—a) yee 14d, (+ 600m SEs 
+3u° *(1- 20) a! B+ Tha 43, 2) 


It is immediately noticed that this expression does not change if: 


(i) ¢ and j are interchanged, 
(ii) & and / are interchanged, 
(iii) 27 and &/ are interchanged. 
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These symmetry properties reduce the number of different components of 
G*'(P, Q) from 81 to 21. So far we have considered 74,(P, Q) as the stress 
tensor in P generated in an infinite elastic medium by a force 8rpe, in Q. For 
what follows it is necessary to consider 7*,(P, Q) in a different way. 

Inspection of (2.7) and (3.1) shows that the subscripts in “ and Gr are 
determined by the given dislocation. In the simplest possible case, an element- 
ary dislocation with only one of the U,’s different from zero, the three Q;,’s 
equal to zero, and & part of a plane perpendicular to a coordinate axis, there 
is only one set of subscripts 77. We also observe that the displacement and 
stress components in Q have subscripts which correspond with the superscripts 
in r*, and G*/. This suggests that 74, and G*/ with fixed 7 and 7 should be con- 
sidered respectively as displacement and stress components, giving altogether 
six different displacement and stress fields. To make this dimensionally correct 
we should consider (U;/u)r{, as a displacement field, but for simplicity of 
argument we shall omit this refinement, and refer to 7‘, as displacement field. 

That each of the displacement fields 7 (P, Q) considered either as function 
of Q(y,) with P(x,) fixed, or as function of P(x;) with Q(y,) fixed, satisfies 
the Navier equations 


(3.5) (A+u)u; ji+uuij, = 0 


everywhere, except when P and Q coincide, follows immediately from the 
symmetry in the suffixes of the Somigliana tensor, the linearity of (2.2), the 
homogeneity and linearity of (3.5), and the symmetry of the variables x ,;— y;. 

To clear up further the nature of the displacement field ri (P, Q), we turn 
back to the displacement field (2.1). This may be considered for fixed 7 and 
variable & as the displacement field in Q due to a force 8zye; in P. This expres- 
sion (2.1) and other displacement fields u*" are usually referred to as nuclei 
of strain (Love 1944). For 7 = / we have a double force without moment; for 
i ~/ it is a double force with moment; if we have 7 = 7 and sum over the 
range (1 = 1, 2, 3) it is called a center of compression. 

It follows that 74 ,(P, Q), which may be written 


(3.6) Tis(P, Q) = —Vb iui —w(ur?+u5*), 


may be considered as the displacement field due to some combinations of nuclei 
of strain; in the case 7 = j this is a center of compression and an additional 
double force without moment; in the case 7 ¥ j, it is a combination of two 
double forces with moment. 

The displacement u,(Q) in (2.6) and (2.7) can therefore be considered as 
being generated by a distribution of these combinations of strain nuclei over 
the surface ¥ of the dislocation. 

4. THE BOUNDARY VALUE PROBLEM, THE GREEN’S FUNCTIONS 
The relations and properties discussed so far all hold in an infinite elastic 


medium. We turn now to the case of the semi-infinite elastic medium. Let the 
surface S of the semi-infinite medium be the plane x; = 0 where the positive 
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X3-axis penetrates into the medium. The boundary conditions we have to 
satisfy are that for each point Q(y1, y2, 0), 


(4.1) 731(Q) = 732(Q) = 733(Q) = 0, 


while we shall in general require also that the stresses and displacements are 
continuous and differentiable everywhere, with the possible exception of 2, 
and moreover vanish at infinity. 

We have already seen that the expressions 7‘ (P, Q) and G%/(P, Q) for fixed 
i and j can be interpreted as displacements and stresses while it is noticed 
that the conditions (4.1) affect the functions GY Gy and G¥‘. It is then clear 
that we have solved the problem of an arbitrary dislocation in a semi-infinite 
elastic medium if we can construct for each of the six combinations 7j functions 
wep wi, Pep and ay which satisfy the following conditions: 

. k , Bip. eee 
(1 ) Wij(Xny Vn) = T14(Xu— Va) £0is (Kas Yn)» 

‘ kl 
(4.2) ty (Xn» Vn) = Gij(Xn—Vn) +H; (Xns Vn): 
(i1) vi, and iy considered as function of y, and the superscripts form a dis- 

placement and stress field regular throughout the semi-infinite medium 


(including the surface 2) so that 


(4.3) Hi} = No tu (vi; +07;). 
(iii) For points Q(y1, v2, 0) we have 
(4.4) Fij(%ny Yn) = Gij(Xn— In) FH (Xny Jn) =O — (U = 1, 2,3). 


The displacement and stress field in the semi-infinite elastic medium due to 
an arbitrary dislocation, while satisfying the boundary conditions (4.1), is then 
determined from the relations 


i U; JJ y Qis jy - k aa v 
Ux(Vn) a Sen . was + Sau (X; Win —Xi Wjn)vnd =, 


Trv(Yn) = HO F* jv, +e “fl (x, Fi'—x, F’})v,d>. 


We shall now show that the functions vf and H‘‘ are obtained by superposition 


(4.5) 


of two solutions: 
(i) an image solution such that the boundary conditions: 


29 


> +31 132 
(4.6) Fy; = Fi; = 0 
are satisfied for (v1, ve, 0), 
(ii) the solution of an ordinary boundary value problem where the surface 


of the semi-infinite medium is subjected to a normal load distribution 
and which takes care of the last boundary condition 


(4.7) Fi; = 0 


in the points (V1, Ve, 0). 
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We have already seen in Section 3 that G{\(P, Q) can be considered as the 
stress in Q generated by certain combinations of double forces in P on 2. 
It is also clear from (3.4) that in case Gj} and Gi are odd (even) in 
x3—ys, the function G%4 is even (odd) in x3;—y3. By putting the same (opposite) 
combinations of double forces in the image point P’(x1, x2, —x3) of P with 
respect to the plane x; = 0, we obtain by superposition of the original and 
image fields a stress field which satisfies the boundary conditions (4.5), while 
no singularities are added in the half-space where the dislocation is situated. 
It is clear that the normal stress is doubled for x; = 0. The solution of the 
boundary value problem where the plane x; = 0 is subjected to normal loads 
equal and opposite to the ones obtained from the double forces in P and P’ 
has to be superposed to satisfy the condition (4.7). 

While this finishes in outline the construction of the Green’s functions, 
several details have to be added. 

Inspection of some of the preceding expressions and papers dealing with 
problems similar to ours (Mindlin 1936; Westergaard 1952) shows that several 
expressions may become clumsy. Following these writers we shall make use of 
Galerkin vectors (Galerkin 1930; Papcovitch 1932). While these vectors do 
not assist a great deal in the actual solution of the problem, they enable one to 
state the results as compactly as possible. As we define the Galerkin vector in a 
slightly different way from some former authors and to make the treatment 
self-contained, we shall discuss these points in the next section. 

The second part of the construction of the Green’s functions is the solution 
of the semi-infinite elastic medium with a normal load, usually called after 
Boussinesq. This problem has received considerable attention in the past (see 
p. 169 of Love 1944). We shall use here a method of Fourier transforms due to 
Trefftz (1935) but modified for normal loads only and put in terms of Galerkin 
vectors by Westergaard (1952). This problem will be discussed in Section 6. 


5. GALERKIN VECTORS 

In analogy with a strain potential, where a displacement field is found by 
taking the gradient of a scalar function, the Galerkin vector I’; is defined as a 
vector which gives a displacement field by differentiating according to the 
Laplace operator and the grad div. We have the relation 
(5.1) Up = Ti nc—ale ci, 
where a will be chosen in a convenient way. If the Galerkin vector is known 
for a problem, the problem is solved, as we have immediately for the strains 
e;,; and stresses 7,,, 
Cis =  ( Ds. eel 5, t0 ad ety 
Tis = ACI —a)6 ijl cute r, etl; ikk) — 2paly kije 
Substitution of (5.1) in the Navier equations (3.5) and selecting a = (A+ y) 
(A+2u) leaves us with the equation 


(5.3) D4 5 9xk = 0. 


This shows that each component of the Galerkin vector is biharmonic. 
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With this definition it is clear from (2.1) that the solution of the infinite 
elastic medium with a force 8zye, in O is given by the Galerkin vector F = re,, 
while the Somigliana tensor is obtained from the Galerkin tensor 


(5.4) li = oir = dar. 


It is also easily found that the Galerkin vector for the displacement f | 
r,(P, Q) with i and j fixed is given by 


(5.5) Ty = Tij(P, Q) = Muraetu(Sar,j+5p7,;). 
If the differentiations to obtain the displacement field are performed to the 
coordinates y, in Q, this expression has to be replaced by its opposite value as 
T, = —r%., 
6. THE BOUSSINESQ PROBLEM 

Let the boundary of the semi-infinite medium be the plane x; = 0 where the 
positive x;-axis points into the medium. According to Westergaard (1952) the 
solution of this problem is determined by a Galerkin vector 


(6.1) lr = ¢3 r3. 


Let the double Fourier transform of I’; be denoted as f;, then we have by 
definition 


T'3(x1, Xo, X3) = = ff P3(Ai, do, x3)exp[1(Aiwi + Aexe) |dAidX2, 


(6.2) fe 
P3(A1, Ag, ¥3) = 32 {J I'3 (x1, Xe, Xs)exp[—2(Aivi + Awe) |dxidxe. 


As T; is biharmonic, we find for 3, using the operational rules, the differential 
equation 


(6.3) (d? dx 3? — B*)? Fs = 0 
with the general solution 
(6.4) T,; = (4+ BBx;3)e"*3+ (C+ DBx;)e**8, 


where we have denoted 8B = ++/(A1°?+A,.?). Owing to the boundary conditions 
at infinity C = D = 0 and we are left with 


(6.5) Ps(Ai, Ao, X3) = (4 + BBx3)e7%*3, 
From the general relation (5.1) and the operational rules we find easily the 


transforms of the displacement and stress components. In particular for 7s; 


and 732 we have 


731 = iud,82{2B(2a—1)e**#—2al’'s}, 
T3223 = ipd28? {2B (2a— 1)e~83 — 2a; }. 


(6.6) 
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As we have only a normal stress at x; = 0, the shearing stresses and their 
transforms have to vanish at x; = 0. This gives us: 


(6.7) A = (2—1/a)B. 


The remaining problem is to find B. Substitution of (6.7) in (6.5) and working 
out the displacement and stress components gives us: 


Tr; = B(2—1/a+ 6x3)e783, 


@, = tad,8B(1—1/at+Bx3)e4"3, 

iy = tadBB(1—1/a+Bx3)e-973, 

i; = —a8*B(1/a+Bxs)e**, 

711 = 26B {wads?+d(1 —a) Ao? — pads?Bx | e-F3, 
(6.8) Foo = 2BB} pads? +A(1—a)A1?— wars?Bx; | e478, 

733 = 2ua8*B(1+Bx3)e—%3, 

To; = —2ipad.B*x;Be*3, 

731 = —2ipad6*x;Be-*s, 

Fi2 = —2uadr,dBB(1—1/a+B8x3)e*3. 


The boundary condition which determines B is the normal load distribution 
on x3 = 0, that is 


(6.9) 733(X1, Xe, 0) = P(X1, X2). 


If the Fourier transform of p is denoted as p, we find B from the transform 
of equation (6.9) at x; = 0. We find easily 


l P(Ay, Az) 


6.10 , Ae) = 
(6.1 B(Ai, A Sua B? 


The remaining problem is the inversion. In this problem the use of the Galerkin 
vector may be advantageous as the inversion of the Galerkin vector may be 
easier than the inversion of the individual expressions (6.8). The displacements 
and stresses can then be found by differentiations once the Galerkin vector 
is known 

This finishes the general theory. In the following section we illustrate the 
method by computation of the displacement field 7/.(?, Q). To prevent con 
fusion we shall adhere strictly to the notation developed in the preceding 


sections, although this is not necessary 


7. THE DISPLACEMENT FIELD 7}.(P, Q) 


k 
\s we have seen in the preceding sections ryo(P?, Q) is essentially the dis 


placement field in VQ venerated by two double forces with moments in ?. We 


choose for /? the coordinates (0,0, ¢) while 2 has the yeneral coordinates 
‘ — vk 
yi(l 1, 2, 3). It follows from (5.5) that the Galerkin vector Py.(P, Q) for 
this case is vyiven by 
a , ‘ ai* ‘ » 3 
J Mie(/ VY) poy, Kk po», 1 


‘ ’ 
where Kk byt t yet t (ys o)* 4? Writing out the displacements and stresses 








ir 


ng 


on 


m 


Is 


Ve 





problem ts 


in Q generated by these double forces in P we find, either by straight differen- 
tiations of (7.1) or from the general relations (2.3) and (3.4), 


1 
T12 


2 
71 


2 


v31 


Giz 


vl2 
G2 


We notice from these relations that the shearing stresses Gy}: and 
in ys3—c. Superposition of the stress field generated by the same double forces 
in P’(0, 0, —c) cancels the shearing stresses at the boundary y; 
doubling the normal stresses G}3. 

The Galerkin vector for the original double forces tn 
P’(0, 0, —c) is therefore given by 


t Mart We? (Vs 
It follows also that the normal surtace load tor 


vi" t Vo" 
It is interesting Co notice Chat the expression 


so that we have tnmedtately 
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= 2u(1—a) 


1%*(1—ay 


Gu" (2a— 1) 


ve 2 V1 
——— (Wau * 





‘ 1 = vm 
tuo (1—a)s3+ by (2a—-1 ‘1 


riz are odd 


he remaining Boussinesg 
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This property also makes it possible to reduce the exponential Fourier trans- 
forms to sine transforms. It is then easily found that 


48° 


(7.6) P(Ay, Ae) = a “tC —a)T,—5ac’T >} ; 


where 7, and 7. denote the integrals 


T, = SJ sin (A1y1)sin (Agy2)dyidy2, 

(7 - ) 00 
T- ff 21%? sin (Ary)sin (Asy2)dyid ys. 

00 


To evaluate these integrals we make use of Basset’s integral for the modified 
Bessel functions of the second kind (Watson 1944, p. 172) and some recurrence 
relations for those functions and the cosine transforms (Magnus, Oberhettin- 
ger, and Tricomi 1954): 


Pie 2 Pie) 35 a 2 2 

J Ko[Aiv (y2 +c )] cos(Azy2)dy2 = 9 Vim i Fy expel eV (Ar + Az2 )], 

oe *. 
K 1[Arv/ (y2" te" )] oA a a 2, \ 2 
vos +2) cos(Asy2)dy2 = 2) vc exPl cV (Ar +A2°)]. 

0 
In this way we finally obtain 
(7.9) pra, Ao) = wand 2 one", 
where we have denoted B = ++/(A;"+Az,’). 

Substitution of (7.9) in (6.10) gives 

a 4u AiA2of l—a = 
(7.10) B(Ay, A2) = _— B whe( 18 oe) - 


while the transform of the Galerkin vector is obtained from (6.8) as follows: 


(7.11) r; = es 


a 


Hs) Ose) Ga 1) gtld—a)(yte) ac] —aeyst fe Bluate) 


It is again noticed that [; is odd in A; and d2 so that the exponential inversion 
integral can be reduced to sine transforms. Working this out we obtain for the 
Galerkin vector of the Boussinesq problem 


Qu S( nee 


aT 


(7.12) T3=+ oF tielievks)~ellle—eend ag 


In this expression J;, /2, and 7; denote the integrals 
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qi 


—4 i) Odds” 12 exp[—8(ys-+c) sin (Ay)sin (Agys), 


Ce he=<4 f f ddxdX2 rs exp|—8(ys-+c) sin (Axy,)sin (Axys), 


=4 | ddydds 3 exp[—8(ys-+)Jsin (Avy, ein (gy). 


00 


I; 


To check the results of these integrations we note that 


a a a a 
(7.14) ay, 2») = — To, ay, 2?) = —Ts, (2 ia oa —I3. 


It follows that J,, Z2, and J; are harmonic functions of y; so that (7.12) 
biharmonic. 
To evaluate the integrals (7.13) we substitute new variables \; = 6 cos ¢, 
= 6 sin ¢, y1 = r cos 0, yo = r sin 6. The integrations to 6 are then easily 
performed as they are essentially Laplace transforms, while the integrations 
to ¢ are elementary but lengthy. Slipping over the details we obtain finally: 


I, = +2 mete -S? + ps! 


(7.15) I, = oe \20- S- ay 
I, = +2029?) —2432_ a 


where we have denoted p = y3+c, ae as before S = {yi2+y2?+ (ys+e)?}?. 
Substitution of (7.15) into (7.12) gives us finally the Galerkin vector 

a atas 4u vivo) 1—a 2 (l—a)(2a—1) (1—a)’ . ; 

(7.16), Ts; = ae ee p- oe . pS+acS 

F cp’ 


rab” (1a)! E E+ at ag 


” 
P cp 
2 € / 
—2ac +3a—~ —2a 
S 
To simplify our further work we make the assumption A = y, which is not 
unusual in applications of elasticity to geophysics. It follows that a = 3. To 
. re k ‘i s “ ‘i 
obtain the displacement components 712 we have to differentiate the Galerkin 
7 k . e 
vector I'y2, which we have now found to be 


Tie = ~u( 2422 ) ’ 


—_ 


S 


(7.17) T= -s( 24%), 
ri, = ee 13" — hr — pS+4cS—8e" +1906 
f prc ob i pie?’ p'cl 
es “s-* gtd Ss 
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Differentiation of (7.17) according to the formula (5.1) gives us the dis- 
placement components: 


1 y 
T12 =a 
r 


\ 


\ 


2 
(7.18) rip = ae A+"r By, 
r 


wha seeBy 
J 
i. 





; 2 (p—97-\ (pb —2 
Tie = p 2epeh gIRe StS) Ip 22) 4449 
an R 
166-1 s,s S4 10l +88; 
fs ee bt. wet 
+6455 80% cs — 408 +405, 
while 
— 2p_4 (b= 2c)", 2 (p—2e)"_ =A 
Am Gig ks ee 1p+s S+3 
ee 29. pe pic pic” 
+e = a ri gi tee 
Saat ee g'h=2 C)" 4 1gP=20)"_ s(0=20)" tgp tos 


ra rae 


satis PO OE 
36 5 40 Ss +405 4265+ 120 $3 


2.2 6 5 42 7 6 2 
pe p pe pe pe pe 
— 120°33- — 8S5 — 120°a5 + 120*25- + 40°37 — 40*a7- 
fg een eo 
One more differentiation will give the stresses. So far it has only been 
checked that Gi; = 0 for ys; = 0. The formulae (7.18) and (7.19) simplify 
a great deal if we restrict ourselves to the displacement components at y; = 0. 
We find: 


1 yo) Vi l 
T1i2 = walt r By ‘7 
7.20 = PN a+tet 
(7.20) T12 = WA): pis 
3 yiye) c 9 < ye’ 
Ti2 = wy 4-18 +26-5— 12-5 ‘ 
r a. # p 
1 = —4e4+2p+2~, 
p 
a °° 
B = 16c+6p-—- 18 +38" 7 1s 


where as before p? = yi7+y2"+c’. 
The discussion of the results obtained here with geophysical applications will 
be taken up in another communication which is being prepared. 
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RADIATION AND SURFACE CURRENTS FROM A SLOT 
ON AN INFINITE CONDUCTING CYLINDER! 


F. H. NORTHOVER 


ABSTRACT 


The far field produced by an aperture in the surface of an infinite cylindrical 
antenna, whose interior is being electromagnetically excited in some way, has 
been calculated by Papas and others using a “‘saddle-point’’ method of approxi- 
mation to the infinite integrals which express the external field. This method will 
not, however, serve to approximate to the far field near the antenna surface. Also, 
there seems to have been no attempt so far to estimate the order of magnitude of 
the error involved in using formulae derived by such methods. The object of the 
present paper is to remedy these deficiencies and to provide a basis for the 
calculation of the far field in wider regions than has so far been possible. In 
particular, the surface current distribution at large distances from the aperture has 
been worked out. Finally, a justification is given for the Papas procedure of 
integrating his approximation to the Poynting energy flux vector over the whole 
of the portion of the large sphere surrounding the aperture which is outside the 
cylinder itself. 


1. INTRODUCTION 

In recent years the slotted cylindrical antennae have become increasingly 
used-in V.H.F. techniques and it is often a matter of some importance to 
know the characteristic behavior of the field many wavelengths from the 
aperture of the antenna. This far field has been calculated by Papas (1949) 
and others, but methods so far used fail for those far points that are close to 
the surface of the antenna.* Also, there seems to have been no attempt so 
far to estimate the degree of accuracy achieved by such methods; nor has a 
justification been given for the Papas procedure of integrating his approxi- 
mation to the Poynting energy flux vector over the entire surface of that 
portion of a large sphere surrounding the aperture lying in the space outside 
the cylinder. 

We assume that the antenna wall is a perfect conductor. 


2. GENERAL EXPRESSIONS FOR THE FIELD OUTSIDE THE 
APERTURE 

We shall follow closely the treatment given by Silver and Saunders (1950). 
Taking axes relative to aperture and antenna as in their paper (see also their 
figure) and using their notation we have the shape of the slot defined by 
curves $:(z) and @¢2(z) and bounded in the axial direction by planes z = 2; 
and 2 = 22. 

The tangential electric field in the slot will in general have both ¢ and z 
components which we consider to be prescribed functions f1(¢, 2) and f2(¢, 2) 
respectively. In view of the infinite conductivity of the antenna wall, this 


'1Manuscript received November 7, 1957. 
Contribution from the Department of Mathematics, Carleton University, Ottawa, Canada. 
*Hallen’s integral equation approach (1948) seems difficult to apply here. 
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represents a complete knowledge of the tangential electric field over the 
entire surface of the cylinder, namely 


Es = fil, 2) 150 < $2,210 2 <K & 
| = (0 outside the slot, 
(1) | Bs=fel6,2) br SO < os 


0 outside the slot. 


Any field component is given by expressions of the form 


Ea(p, %, 2) = \ é.(p, ¢, h) dh, 


where 


z_Je i dH (Ap) nom y pp [eet sis 
| 6, = | thAy ap — B,H, (Ap) ce, 
aH? (, Ae) cil a 


: é LZ nh 

(2A) ; &= 1217 = 5 ind ( Ap) — iauBy "> 
Sx 
\¢ 


at 
wt 
Il 


seni 4») | oY a 


oe 


where 


(2B) A = (k?—h?)}, k = (w2pe+inow)}, 
and A, and B, are given by 


co 1 ‘iaaliad \ ~~ inB— int 
. A. = 4 ‘7 —}? YH? ( ah a. f2(B, &) é dB, 


22 2(§) 


9C 4 Sie’ “ciscmentaiecte mh iepediiocicaies — inB—iht 7- 
(20) | Bs = Fyoua (eh) (GH (Aa) /0a) J, *) or" ” 
| a ait 4 i. ; sea 

+ FF OH (Aa)/da} J, Spy) a8 


This is equivalent to 


ep2(t) 
E, -f J 6, dt dB, 
o1(t) 
a2( 1) 
(3) _ EK= is &, dt dB, 
o1(0t) 
epe(t) 
E, = 5 a é, dt dB, 
o1(0) 


in which adtd8 is an element of area of the surface of the aperture, dt being 
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parallel to the z direction, and d8 being azimuthal, and the integration is 


over the slot surface, where 


: ‘ 2 (°). @Hi?(Ap) now @ \ ino+ itz 
| 6-3 — pte) — Bat a ETO (Ap) eet ap, 


| (1) 
a) | dem Si f4 2 (ap) — iad, SE Ae pmeestgy 


~ 

a 
n 

II 


= , &e —f°) aH? (Ap) e*** "dg, 


where 
_ fa (8, t) eo itt 
oo ta ‘e—P) (he) ’ 
(3B) 4 —inB—itt e 
hao 


de” 


Sorgen aN 2) fe (8, t) +f1(8, yt. 


In the above f; and fz are as defined in Silver and Saunders’ paper but 
here we take the time factor as e~“‘ in order to obtain (after Papas) results 
in terms of Hankel functions of the first kind. The path of integration is 
along the real axis of the cut plane defined by the cuts 


—3n/2 < arg({—k) < 7/2, 
—n/2 < arg{¢—(—k)} < 32/2. 


It is not necessary to write down the corresponding expressions for H as 
these can be obtained from the above via the fundamental relation 


(4) —ikH = curl E. 


Inspection of the above shows that the field intensity E can be expressed 
in terms of fundamental integrals 


Ie J Hy {(k?— =5ie p} eh Ode 


WH | (R27) 'a} 


Hy? {(kR°- =o} __ gif(2-D gy 


(5) P= J omna —1)'e}/aa) 
Yoo (1) , 
ody 1 He iG ‘oy ites 
| : - HO a 


Hy” {(k°— J cy p} Rs. it (z 
PP Bier “yr 


t) 
at, 


| H = Ode 


as follows* 


*No confusion should result from our different uses for the symbol E. 
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~ 1 Ste ee 2 ae 
E-7abe JJ Lapaz pa oz $1268 #) 


- - Fif (8, »| e " dtdg, 


' = by Eo : 
(6) Es - 4° x ? Je = p Oz a Opdz frlB, t) 


4 Ifa 6 0g, | « "atdB, 
Op 


a0 


* ree, a i Sf mp0, the ™dtdB. 


These fundamental integrals are connected by the equations 


(7) f @G/d2? + RG = E, 
\ 9H /d22 + RH = F. 


3. THE FAR ZONE FIELD 
It is convenient to describe the far-zone field in terms of spherical polar 
coordinates 7, 0, @ and, accordingly, we shall later make the substitutions 


(8) p=rsing, z=rcosé, 


In the work of approximation to the far field the case when the field point 
is near the antenna surface has to be treated separately from the case when 
this is not so. The reason for this will appear later. 


3.1. THE FIELD AT POINTS NOT TOO NEAR THE ANTENNA SURFACE 


This case is dealt with by Papas (1949) and Silver and Saunders (1950) 
by applying the saddle-point method of asymptotic approximation to the 
infinite integrals concerned. The expressions for the field given in the Silver 
and Saunders paper are first approximations only: the analysis necessary to 
estimate the order of the error involved in the use of these expressions is too 
lengthy to include here. It appears in a contemporary D.R.B. report by the 
writer (Northover 1958). The results are: 

(1) when kaé@ is not small the Silver and Saunders expressions must be 
multiplied by a factor of 

1 + O(a?/rx) 
to obtain the true result, 

(2) when ka@ < 1 but not so small as to make this type of approximation 
invalid,* their expressions have to be multiplied by a factor 


1 + Of1/[(Rr)70*}} 


in order to reach the accurate result. 
For small values of ka@ the Silver and Saunders expressions are thus correct 
to order 1/(kr6?) (i.e., the error involved is at most Of (kr)-?6-4}. 


*I.e., we must have 62(kr) > 1 as is shown in the report. 
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APPROXIMATION TO THE FIELD FOR z LARGE AND p MODERATE 
3.2.1. The integrals to be approximated to are E, F, G, and H. We begin by 
considering the integral E. 
Integral E 
This integral is taken along the real axis with an upward indentation at 
¢ = —k and a downward one at ¢ = &, in the cut plane defined by 


— - < arg((—k) < 
(9) . 


a 


- ae <arg({+k) < 57 
We write for brevity z instead of (s—¢) (as the form of the integral is 
thereby unchanged) and assume first that z > 0;* then, owing to the con- 
vergence to zero of the integral over the arcs of the circle |{—k| = R defined 
by —a > arg({—k) > — $a and 0 < arg(¢{—k) < $27, as R- ©, we may 
deform the path of integration away from the real axis into a path passing 
round ¢ = k and going to infinity on opposite sides of the cut which starts 
from this branch point. In fact we deform the path of integration to the 
path shown in the figure: 


As regards the contribution to the integral from the small circular part 
surrounding the branch point k, it is easy to see that this is at most of order 
e, where ¢ is the radius of the indentation. Hence, writing 


¢ = k+in/z 


and noting that kR—¢ = e'"(¢—k), the integral becomes (after using well- 
known properties of Hankel functions and reducing) 


| 2 sus (°° Jnlor) Yaa) — Valor) In(a) 
i? _ Hg J(a,)+ Y,(a1) , dn, 


(oa) 

p 

10A 

(10A) { as) oon vi)! 
a Tbe 


We now require the asymptotic forms of this integral for z large and p 


where 


a, 
> 
ll 


a 
~ 


*The case s < 0 would be dealt with by a corresponding deformation around the branch 
point ¢ = —k. 
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moderate. The analysis falls into two separate cases, namely, » ~ 0 and 
n = 0. We consider these in turn. 
Case 1.—n # 0. For — small and n > 0 
, ; 2” (n—1)! 
(11) 1.) % ==. 
Note that 
H°,(¢) _ H,"(é) 
HE) He(€) 
so that it is only necessary to consider positive values of n. Therefore the 
above integral approximates to 


Qe"? - a 1( e) e 3 x mw = 
do met\a) “\p/S 2 {a—py 
(12) ~— 2s e** 7-"(1kab)" oe , 1 =f : y" 
é z 7 (n—1)! p/ § 


as a first approximation, on taking 
(13) a, = a(2kn/z)* eo *"*. 
Case 2.--n = 0. Integral is 


2 ike r Jo(pr) Yo(a1) — Volos) Jo(ai) 9 
Jo Jo(ai)+ Yo(a1) : ’ 


For & small 
(14) Yo(&) = (2/m) log é. 


Hence the integral approximates* to 


= 4a ik2 J 1 \ ? 
(15) ae e( * ) troe( 455) ; 


Integral F 

Deforming the path of integration into the same contour as that used 
above, the integral becomes equal to 
2 a i Jn(p1) Y,’ (a1) cee Y, (91) Jn’ (a1) a ee dn 


; (T,’(ai)\?+4¥n7(ai)}>— (2kn/z)' e (1 +-ni/2kz)? 


~ 


2 KZ -—-n n ~ 
(16) ~ — See et; + (shad) {1+(¢) (  (n>0), 


~ 
~ 


= = & 
(17) sia? (n = 0). 
on carrying out an analysis similar to the preceding analysis. 


The Integrals of G and H 
It can be seen by inspection that deformation of the path of integration 
into the contour used above is not now permissible since these integrals would 


*Cf. D.R.B. report or Appendix 1, 
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not converge along such a contour owing to the presence of the factor 
1/(k°—¢*) in the integrand. We therefore proceed indirectly, using the 
relations (7). We illustrate the method of dealing with these two integrals 
by giving the analysis in full for the case of G. 

Case n ¥ 0.—We have 


o°G 27 = 2m ike -—m4 ie es j (¢ | 
(18) az + kG srt Crt (Skap) (n—1)! \! . ( i 
Put 
(19) G = ey, 
then 
” . ’ 2a “—n n 1 j a ™) 
(20) x" +2ikx’ ~ — -nzit "(dkap) atntt-(*) ( 
and therefore 
‘ T “—n 1 n J a a | 
oh) X™~ ike’* nl (2kap) e=(2) (> 
so that 
’ T ikz—jnri Ms n J a a} 
(22) G ~ amen ; ni (2Rap) bts {- 
This makes 
(23) x” /kx’| = O(1/kz) K 1, 


thus justifying the neglect of x’ on the left of (20), which led to (22). 
Case n = 0.—-Here 


Me. os. dm ik: 2) /} ( lz yy 
2. ; I~ 4 r % 
(24) ae + k( e€ oe Ss \ log mr { ; 


Similar reasoning to the above gives 


in . - 2e a era. 
(25) G~ ib oe( * \log a Ob {- 


Proceeding similarly for H, we have 


ayn nd T ikz—jnvi _ 1 1 n ) ( a yr 
(27) tl log z (2 = 0). 


~ Qik 
3.2.2. Approximate Evaluation of the Far Field 


Using the above approximate expression, we have 


as Ge estes 1 . = (2 )" 
28 = agit jnri 1 n on 0 . 
oe dz0p —_2iz" , (n—1)! (2kap) \ p f ee) 


2nJ y Pi © ) _ 
(29) ~~ - yl woe (4 ( (ser), 
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0H TW tkz—jnr 1 / n J a ~ 
(30) fae an ge siali (zkap) i+(2) fo 


0z z n( 
‘ dG W oikz—}nri 1 1 n J a ( a , l 
(31) as 2" é n! (hap) ‘fe ? f , 
- aH x wka tez—-yart ly 2. (s ¥ 
(32) Opdz 22" ° n! (3kap) er p me 
eau oF rka 1 a} 1 J (s . 
(33) ap atl (4kap) eS it ( (n > 0), 
(34) ~ O(1/2?) (2:-= 0). 
Writing 
(35) y = JJfe(B, t) dtdg, 


the double integral being taken over the aperture, the first approximation 
to the field for z large and moderate is as follows: 


tk2 


E,~51-;* 
“P — 2xp log Z’ 


(36) { Ey~ of 4) 
| ¥ on «) 
me rz (log Z)° log ( p/’ 


(36A) Z= 


where 


This gives 
y ao 
(37) Eye~ - ee 
7 2mp log Z 

It will be noticed that the field is dominated by the mode u = 0: the 
corresponding expressions for H (as obtained from (4)) are 


of ‘y ’ 
ikz 


Seg 
2p log Z’ 


| anidh ah 


4. CONCLUSION 


| H, 


(38) { Hy~ - 


The main interest in the elucidation of the field near the metal surface of 
the antenna is that this allows the surface currents at long distances from 
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the slot to be investigated. For if H, is the surface component of H and j is 
the conduction current vector on the surface, j is perpendicular to H, and 
4e|j| = |H,|. 

At the ends of a finite antenna the component of the surface conduction 
current perpendicular to the bounding edge is zero. By superposing suitable 
extra terms on the surface current terms appropriate to the infinite antenna 
it may be possible to realize this condition and thus to treat the case of a 
finite antenna. It is hoped to follow up this idea in a later investigation. 

The appendices deal mainly with the question of the range of validity of 
the two types of approximation and with that of the validity of the Papas 
procedure of using the saddle-point type of approximations for E and H 
for the Poynting energy flux vector in an integration of the normal component 
of this vector over an outside surface completely enclosing the radiating 
aperture. 
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APPENDIX 1 
Conditions under Which This Type of Approximation Is Valid 
It is sufficient to consider a typical integral, viz., 


Ho’ {(k’—$")'p} ge 
J ~ He {(k’—¢")*a} eG. 


This is expressible as 


a o tke { Ji (pi) Yo(ai) — Yo(p1) Jo(a1) e ” 1 
zZ oF Jy"(ay) + VYo*(ay) Ca 


We split up the range of integration into 0 < » < K, » > K, where K is 


taken large enough to ensure that 


ey | Kk 
Jk | JG 


but small enough to make the approximation 


T(pi) Yolay) — Yo(pi) Jo(ay) 5 i log(a/p) 
2 7,2 a eee 25 
Jo (ay) + Vo" (a;) wil+(4/n’)flog(a;VYa)} 
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(where a = je~*”, y being Euler’s constant) good throughout the range 
0<1< K. This requires* that 


z > 2kp*K, 


while K has to be taken large enough to make 


| eo | | aK 
| | «K] | 
1 VK | 1 vo | 


We then wish to know under what circumstances the following holds: 


wK =< 
J ae e "dn a oe 
o 1+(1/n°*){log(n/z)— Bir} ~ 14+(1/2°) (log Z+4ix) 
wherein we have used the formula (13) for a, and written 
Z = 2e*1(2/ka*). 


We further split the range 0 to K of the integral {* into the ranges 0 to Z-! 
and Z-! to K; the first integral here is obviously at most of order Z-}(In Z)-? 
and the second approximates to 


—g71/3 _K l 


_ ) Ses 5 : —<5 
oa . 1+(1/r°)(In 2+ 4297)” 


since the integrand in this range approximates to 
7 
4 


I+(1, e*)(In Z+4ir)" 


(with K < Z!, say, Z large). Since K is not small, i. is evidently O(e-* 


and safe criteria for the desired approximations are 
Gs). 2S 2F- 
(ii) A such that e-* « #?/ (In Z)?, 
(ii) s > 2kp?K. 
(i) and (iii) are satisfied for sufficiently large Z and (ii) is satisfied if, for 
example, 
K > 2ln(In 2). 
If we take 
2 > 4e"K, 
Le., Say 


Z > Se” In(In Z), 


then this requires Z > at least 30 (and A at least about 2.6) 
Hence a sufficient condition for this approximation to be valid is that Z > 
about 80. 


* 


mo <I throughO go» CA requires that 2 > 2ep*A 
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If, further, Z > about e**, then 


1 T 


and the simple type approximations we have used in the text are then valid. 


APPENDIX 2 
On the Validity of the Papas Procedure of Integrating the Poynting Energy Flux 


Vector 
Papas uses his expressions for E and H to evaluate 


JJ(EAH).dS 


over the whole of that part of an infinite sphere outside the cylinder and 
containing the radiating aperture, whereas his approximate expressions are 

. + pk —} : ; = 5s as | 
only valid* for 6 > (kr)-?. However, by the radiation conditions, 7/E! and 
r|H| are finite as r — ©, so the error involved is at most of order 


~wVKi/vr 
J Bo do 


a/r 


where B is a constant, and this integral approaches 0 as r > «©. (/E A H.dS! 
will be at most O(r~? 27r? sin 6 d6).) 


APPENDIX 3 
On the Surface Value of H, for Large z 
It has been conjectured that the saddle-point type approximations give 
the correct value for the surface component of H (tangential component of 
H at the surface of the cylinder).t However, this is not the case. 
At the surface, E, = 0 and we consider H,. We have 
dE,  9E, 


Oz Op ’ 
iit. a 2 Oey 
p Op 


tkHH, = 


Therefore 


iky Ps v gt 
ig 00 i re ae 
. 2mp in Z mzp (In Z)” 
_ iky an 


2rpinZ’ 


where 
VS Sffe(a, t) dt dB. 


*Cf. the D.R.B. report (References). 
This surface component is of importance in calculating the surface current. 
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This gives 


ikH, ~ O(1/2?). 
Also 
ikH, = —0E,/d2 = O(1/2°). 


However, the saddle-point type expressions give 


E, ~ ied fo(B, t) = 25 cot 6 ar (H\?(kad)} ~'d6 at 


* ee 
4p In(kap/z) ’ 


, x anf 2 i ey ff 
Ey ~ Qn + x kad filB, 1) dB dt 


27 tar 
Rim is ff fo(B, t) dB dt X & ja 


tkr 


e e ; 
4a In(kap/zs) ’ 


and these give, for the surface value of H, 


H, ~ O(1/r?), 
H, ~ O(1/r’), 


_ tky ao iky ant iky e* 


neh ~ dep In(kap, 3) 4a In(z ka*)  4rainZ’ 


The Papas approximations therefore give half the correct result. Wait has 
suggested the following physical explanation: ‘the Papas result when extra- 
polated to the cylinder surface should be multiplied by two to account for 
the imaging of the conductor”’. 


APPENDIX 4 
Gap in the Range of Applicability of the Approximations 
The theory of this paper shows that neither of the two types of approxima- 
tion to the field developed can be used in the intermediate range 


0.435(kr)- < 6 < (kr)-} 


(taking K about 2.6 in Appendix 1). 

Unfortunately considerable difficulty is encountered when attempts are 
made to close this gap. Fortunately, it is relatively narrow: however, if the 
field was required therein, there would probably be no alternative to heavy 
numerical computation. 











HIGH RESOLUTION RAMAN SPECTROSCOPY OF GASES 
XI. SPECTRA OF CS, AND CO,! 


B. P. STOICHEFF 


ABSTRACT 


The vibrational Raman spectra of CS2, C02, and COs, consisting of the 
strong Fermi diad »;, 2v2, have been photographed with a 21 ft. grating. In the 
spectrum of CS, 12 additional sharp Q branches were observed in the region of the 
diad; three are due to isotopic molecules and the remainder are ‘‘hot’’ bands. 
The rotational structure of the strong v; band was also obtained. These measure- 
ments together with infrared data are used to determine the vibrational con- 
stants w;° and x; of CS. The pure rotational spectrum of CS», with rotational 
lines up to J = 94, yields the constants Booo = 0.10910+0.00005 cm, Dooo = 
1.0X10-§ em, and ro(C==S) = 1.5545+0.0003 A. For COs, the rotational 
structure of the diad was analyzed and the results are in agreement with recent 
infrared data. 


A. INTRODUCTION 


One of the vibrational Raman spectra most frequently studied in the past is 
that of the CS, molecule. Interest in this spectrum was stimulated primarily 
by its unexpected complexities, for example: the Fermi resonance between »; 
and 2v,. (Placzek 1934), the appearance of the forbidden bands v2 and v3 (Wood 
and Collins 1932; Venkateswaran 1932; Evans and Bernstein 1956), and the 
appearance of “hot’’ bands (Langseth, Sorensen, and Nielsen 1934) and of 
vibrational bands due to isotopic molecules (Giulotto and Caldirola 1941). All 
of these investigations were carried out with liquid CS.. The Raman spectrum 
of solid CS. has been studied by Sirkar (1936), and that of the vapor by 
Imanishi (1935), who observed the »; and 22 Q branches. 

In order to unravel the complex vibrational spectrum, Langseth, Sorensen, 
and Nielsen (1934) and Giulotto and Caldirola (1941) used spectrographs of 
relatively large dispersion (7 A/mm.) and were successful in interpreting most 
of the observed details. However, as will be evident from the present results, 
much more information can be obtained from the Raman spectrum of the 
vapor particularly when the spectrum can be photographed at high resolution. 

The infrared absorption spectrum of CS, vapor was studied by Bailey and 
Cassie (1933) and the rotational structure of one of the infrared bands was 
first resolved by Sanderson (1935). The ultraviolet electronic spectrum was 
photographed under high resolution by Liebermann (1941) and yielded accu- 
rate values of the rotational constants and of the vibrational interval 2v.. More 
recently, several of the overtone and combination bands have been observed 
at high resolution by Gailar and Plyler (1952) and Allen, Plyler, and Blaine 
(1956). These investigations provide further information about the rotational 
and vibrational constants. 

\Manuscript received October 30, 1957. 
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The vibrational Raman spectrum of CO, has also been studied frequently, 
and as in CSs, the Fermi diad »;, 2v2 has been the most interesting feature in 
the spectrum. But in contrast to CS», the infrared spectrum of CQO, has also 
been the subject of a large number of investigations. One of the most extensive 
of these is the recent study by Courtoy (1957). During his research, further 
data on the diad »;, 2v2 were desirable and one of the possible means of 
obtaining this information was by photographing the Raman spectrum at 
high resolution. 

Up to the present time, this series of papers has dealt with the pure rotational 
spectra of molecules and not with the less intense rotation—vibration spectra. 
In contrast, the primary purpose of the present investigation was to attempt 
to photograph the vibrational spectra of CS. and CO: with a 21 ft. grating and 
if possible to observe the rotational structure of the most intense bands. 
Fortunately, the vibrational spectra of these molecules are relatively intense 
and several plates of the rotational structure of the vibrational bands were 
obtained in reasonable exposure times. These have vielded useful information, 
particularly about the role of Fermi resonance. Also in the CS: spectrum, 
sufficient data were obtained from the “‘hot’’ bands to permit the calculation of 
vibrational frequencies and anharmonic constants. In addition, the pure 
rotational spectrum of CS. was photographed and has given precise values of 
the rotational constant and internuclear distance for the ground state. 


B. EXPERIMENTAL PROCEDURE 


The apparatus, including the mirror-type Raman tube, mercury lamps, 
and 21 ft. grating spectrograph, has been described in the first paper of this 
series (Stoicheff 1954). For the pure rotational spectrum of CSe, the cell was 
filled to a pressure of 10 cm. Hg. Exposure times of 4 hours were used with a 
slit width of 0.02 mm. For the vibrational spectrum of CS., pressures up to 1 
atmosphere were used and the cell was kept at a temperature of 60° C. A filter 
of NaNOsz solution was circulated around the cell. Since the vibrational spectra 
are much weaker than pure rotational spectra in the Raman effect, a cylindrical 
lens was placed in the spectrograph, in front of the photographic plate as 
described in paper VIII (Callomon and Stoicheff 1957). This reduced the 
exposure times by a factor of about 15. Several plates of the vibrational spec- 
trum were obtained with exposure times of 4 to 20 hours, using slit widths of 
0.02 and 0.05 mm. 

The spectra of COs and (of a 60°% enriched sample) of COs, consisting ot 
the sharp Q branches »;, 2v2, were photographed with gas pressures of | atmos- 
phere, using the cylindrical lens, in exposure times of 4+ to 8 hours. Exposures 
lasting 80 hours were necessary to photograph the rotational structure of the 
vy, and 2p. bands of COs. 

All of the spectra were photographed in the second order of the N.R.C. 21 ft. 
grating spectrograph (dispersion of 1.25 A mm.). Second- and third-order iron 
lines were used as standards. The plates were measured with a photoelectric 
scanning comparator built in this laboratory, after the design of Tomkins and 
red (1951). 
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C. ROTATIONAL SPECTRUM OF CS, 


A photograph of the pure rotational spectrum of CS, is shown in Fig. 1(qa). 
The Stokes and anti-Stokes branches extend to 40 cm~ from the exciting line 
and include lines from J = 10 to J = 94. The line spacing is 0.87 cm~!. Two 
plates were evaluated and the average values of the wavenumber shifts, 
accurate to +0.02 cm— or better, are listed in Table I. 


TABLE I 


WAVENUMBER SHIFTS (UNITS CM~!) IN THE ROTATIONAL 
RAMAN SPECTRUM OF CS, 











Ap a 











J Ave I 
10 4.958 50 22.469 
12 5. 868 52 23.343 
14 6.761 54 24.217 
16 7.637 56 25.091 
18 8.497 58 25.959 
20 —_—_—* 60 26.826 
22 10.247* 62 27.695 
24 11.125 64 28.560 
26 12.002 | 66 29.423 
28 12.871 68 30.299 
30 13.745 70 31.176 
32 14.615 72 32.037 
34 15.483 74 32.909 
36 16.355 76 33.790 
38 17.230 78 34.657 
40 18.109 80 35.529 
42 19.003* 82 36.401 
44 19.839 84 37.271 
46 20.723 86 38.151 
48 21.589 88 38.997 
| 90 39. 862 
92 40.742 
| 94 41.594 





2All values are averages of the Stokes and anti-Stokes 
lines measured on two plates. 
*Lines blended with grating ghosts. 


Since the CS, molecule has a center of symmetry (point group D,,,) and the 
most abundant isotope, S*, has a nuclear spin of zero (Bose statistics), alter- 
nate rotational levels are missing. Since the ground electronic state is 'Z}, it is 
the odd rotational levels which are absent. The selection rule for rotational 
transitions is AJ = 2 and leads to the expression 


|Av| = 4Booo(J+#) —8Doo0(J +3)’, 


for the displacements of the rotational lines from the unshifted Rayleigh line 
(AJ = 0). The rotational and centrifugal distortion constants were evaluated 
in the usual way by plotting a graph of |Av|/(J +3) against (J +3)? as shown 
in Fig. 2. The values* obtained are 


FRE Bp. 5-05 


*At the Symposium on Molecular Structure and Spectroscopy, Ohio, 1957, A. H. Guenther, 
T. A. Wiggins, and D. H. Rank reported the results of their investigations of the 3»; and 
vo+3v3—ve2 infrared bands. They obtained Bo = 0.10910+0.00002 and Do = (1.10.3) 
X 10-8 cm, in good agreement with the present values determined from the Raman spectrum, 
which were also reported at this symposium. 
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Fic. 1. Raman spectra of CS: and CO, photographed with a 21 ft. grating. (a) Pure rota- 
tional spectrum of CS». (6) Fermi diad »;, 2v2 of CS, photographed in exposure times of 4 hours 
and 20 hours. (c) Q branch lines of the Fermi diad »;, 2v2 of C?0O2 and C¥O.. (d) Rotational 
structure in the diad »;, 2v2 of C?Oz. 
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Booo = (0.10910+0.00005 cm=! and Dooo = 1.0X10-° cm. 
The effective moment of inertia is therefore 


27 gR¢ ‘i 
he= re"? x 10-" = (256.54+0.12) 107" g. em? 
000 


and gives 


ro(C=S) = 1.5545+0.0003 A. 


0.4360 





(J+3/2)? 


Fic. 2. |Av|/(J+3/2) graph for the ground state of CS. (X represents rotational lines 
blended with grating ghosts). 


The present values of the constants are considered to be an improvement 
over the recent infrared values given by Allen, Plyler, and Blaine (1956). Their 
value of Do = 4.341078 is high and is reflected in their high value of By = 
0.10927 cm~'. The relation Dy = 4By'/w,;? should give a reliable value of Do. 
With By = 0.109 and w, = 659 cm~', one obtains Dy = 1.2 10-8 cm—, which 
is almost the same as the value obtained in the present research. 

In order to determine B, and hence the equilibrium internuclear distance 
in CS, the quantity }(a@a:+2a2.+a;) must be known. Its value can be obtained 
experimentally from the convergence of rotational lines in the band v1 +272+ 3, 
or in three bands which include the vibrations v1, v2, v3. For CSe, the infrared 
band v;+2v2.+v;3 at 2961.7 cm has been investigated by Allen, Plyler, and 
Blaine (1956), who report a value of a;+2a.+a; = 0.000460 cm~. Un- 
fortunately, there appears to be an error in the measurement of this band 
since it gives high values of By and Do; there is therefore some uncertainty 
(about 10%) in their value. Consequently larger errors than those found for 
By and ro are assigned to the values of B, and r;: 


B, = 0.10933+0.00008 cm-', r, = 1.5529+0.0005 A. 


D. VIBRATIONAL SPECTRUM OF CS:z 


Photographs of the vibrational spectrum of CS, are shown in Fig. 1(b). The 
lower strip shows the strong vibrational lines photographed in 4 hours, and 
the upper strip shows the rotational structure of the »; band as well as other 
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weaker vibrational Q branches photographed in 20 hours. The vibrational 
spectrum is indeed a complex one consisting of the prominent Fermi diad 1, 
2v. and including three Q branches of isotopic molecules as well as nine ‘‘hot”’ 
bands. 


(1) The v; Band 

The most intense line in the spectrum is the Q branch of the »; band. The 
weak rotational structure on either side of this line forms the S and O branches 
of this band. The wavenumber shifts are listed in Table II. 


TABLE II 
WAVENUMBER SHIFTS® (UNITS CM7!) OF THE LINES IN THE yy BAND OF CS. 


J S(J) OW) x S(J) OV) 

8 - 654.90 42 676.67 639.62 
10 662.92 653.77 44 677 .68 638.77 
12 663.79 652.99 46 678.42 637.76 
14 664.72 652.13 48 679.13 636 .94 
16 665.57 651.29 50 679.99 636.04 
18 666.25 52 680.83 635.04 
20 667 .22 54 681.65 634.21 
22 668.16 649.11 56 682.56 633.41 
24 668.99 58 683.33 632.49 
26 669.86 646.80 60 684.17 631.74 
28 670.68 645.79 62 684.96 630.58 
30 671.52 645.02 64 629.62 
32 672.39 644.11 66 628.49 
34 673.26 643.21 68 627 .86 
36 674.13 642.33 70 626.87 
38 675.00 641.29 77 626.00 
10 676.67 640.62 74 625 .22 








“Average values of measurements of two plates. 


The formulae of the rotation—vibration lines are given by the well-known 


equations 
S(J) = vo +6B’ + (5B’ — BY’) J+(B’—B") J’, J=0,2,4,6..., 
O(J) = »9+2B’ — (3B’+B")J —(B’— BY") SJ’, FBG os 6x 


These can be written as a single equation: 
vy = v1 —2B"+(B'+3B")m+ (B'—B”)m’, 


where m = J+2 for the S branch, and m = —J+1 for the O branch. A graph 
of vy versus m gave the relative numbering of the lines in the O and S branches, 
and the absolute numbering was then determined from the known value of 
B” and the assumption B’ = B’’. The band origin vp and (B’—B’’) were 
determined from a graph (Fig. 3) of the combination sum 


31S(J) +O(J)] = vo+4B’+(B’—B")J(J+1). 


The values are vp = 657.98 cm~'! and (B’—B”) = —0.000155 cm~. This 
value of B’— B” gives —a; toa very good approximation (see part (2) of this 
section), since the Fermi resonance of the (10°0) and (02°0) levels is weak. 
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This fact also explains the value D’— D” = 0 observed for this band, while in 
COs, where the resonance is very strong, values of D’ # D” have been ob- 
served for many of the interacting levels (Courtoy and Herzberg 1955). 


cm-! 
S(J)+0(J) ede ae 
658.40 SOOO) Uy + 48 +(B'-B") J (J+!) 
658.20 
658.00 





° 1000 2000 3000 4000 


J(J+1) 
Fic. 3. A graph of 3[S(J)+O(J)] for the »; band of CSs. 


(2) The Vibrational Spectrum and Vibrational Constants 

The wavenumber shifts and assignments of the vibrational lines are given 
in Table 111. The assignment of the isotopic bands was carried out with the aid 
of known isotopic abundances (CS.” = 90%, CS®S* = 1.4%, CS®S* = 8%, 
C8S,.22 = 1%) and with the aid of the expected isotopic shifts Av;(CS®S* = 
—5 cm, CS®S* = —9.7 cm—'!, and CS,” = —2.5 cm —'). For the assignment 
of the “hot’’ bands, the Boltzmann factor e~”’/*? and the known degeneracy 
of the initial levels must be taken into account. All of the assignments given 
here obey these requirements except for the band assigned as (271+ v2) — (1 
+v2), which appears to have almost double the expected intensity. At the 
bottom of Table III are listed the known infrared fundamentals and overtone 
bands. In Fig. 4, all of the observed transitions are represented in an energy 
level diagram. 

The vibrational levels of a linear triatomic molecule (with respect to the 
lowest vibrational state) can be represented by 


(1) Go(v1, Vo, V3) = 1°01 Fw 2Vo+ 3 °V3 HV +X 2202? +N 33037 FX 120102 
+. 130103 +X 2302034 Zool 


if no resonance occurs. In Eq. (1), /2 is the quantum number of the vibrational 
angular momentum about the symmetry axis (/2 = v2, v.—2,..., 1, or 0). 
Since Eq. (1) has 10 unknown constants, w;, «ix, and g22, 10 of the Raman and 
infrared bands, marked with an asterisk (*) in Table III, were chosen for 
calculating these constants. The values obtained are 





w," = 659.23, Nig = — 12:25, Xy = —7.74, 
(2 we? = 396.55, Kee = 2.25; Xo; = —6.67T, 
ws" = 1537 .637, X33 = —5.13f, x13 = —4.96, 

2 = —2.00. 


Only the values marked with a dagger (f) are not affected by Fermi resonance. 
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Fic. 4. Vibrational levels of the ground electronic state of CS». The vertical solid lines 
connecting various vibrational levels (horizontal lines) represent observed Raman transitions: 
the broken lines represent cbserved infrared transitions. 


The constants were then used to calculate the remaining nine bands. Fair 
agreement is obtained as shown in column three of Table III. The deviations 
are due in part to the use of Q branch measurements instead of band origins 
and to the limited accuracy with which some of the vibrational frequencies 
are known. But the greatest part of the deviations are presumably due to the 
neglect of Fermi resonance amongst many of the vibrational levels. 

The occurrence of Fermi resonance between the vibrational levels (10°0) 
and (02°0) in CS. is most clearly shown by the appearance of a relatively 
strong Raman line at almost twice the frequency of the fundamental v2 band. 
To be sure, even without Fermi resonance, a weak line is expected at this 
frequency in the Raman spectrum; however, its intensity would most likely 
be much less than 1% of the »; band and would only be observed with long 
exposure times. In contrast to this expectation, the observed intensity of the 
2v2 band is between 5 and 10% of the intensity of the »; band. For this reason, 
and since a satisfactory analysis of the vibrational spectrum is then possible, 
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TABLE III 
VIBRATIONAL BANDS OBSERVED IN THE RAMAN AND INFRARED 












vy (cm~!) Assignment O—C (cm) 
Raman bands 
637.76 (vy +3p2!) —3p2! 3.0 
641.29 (yi tre!) —py! CSS34 
643.64 (vy +2y2?) — 2v-? 1.14 
645 .02 (v1 +2y2°) —2y2° 2.52 
648.01 (2v; + v2!) — (vi +2!) 0.26 
648 . 37 v, CS2S# 0.07 
650.24 (vy +2!) —ve! 0 * 
652.99 v, CS8S33 0 
655.49 2vi—v1 0 > 
657 .98¢ Vy 0 * 
802.11 22" Go * 
811.12 3yy' — ve! Oo * 
818.15 Av.” —2v2° —1.95 
819. 18 4vo? —2p.? —0.92 
Infrared bands 
396.82." va! = © 
1532. 5° « V3 oe * 
2185 . 522.4 vitv: 0 
2833. 142.4 2vitvs —2.90 
2961 . 732.4 vp t+2vo+v3 —1.49 
4566. 694: 3y3 q 
4546.74. (vo! +3y3) — ve! 0 * 
0.5 


4527 .22.°¢ (2v2 +33) —2ve 





*Bands used to calculate vibrational and anharmonic constants 
listed in (2). 

“Band origins; other values of » represent measurements of sharp Q 
branches. 

’Measurements by Bailey and Cassie (1933). 

‘Measurements by Gailar and Plyler (1952). 

4Measurements by Allen, Plyler, and Blaine (1956). 


the band at 802.11 cm™! is considered to arise from the transition (02°0)— 
(00°0), which gains intensity through the interaction of the upper level (02°0) 
with the nearby level (10°0). 

The Fermi resonance in CS» is, however, not strong, as is evident from the 
large difference in intensity of the »; and 2v2 bands. A good estimate of the 
positions of the unperturbed levels and of the interaction energy can be 
readily obtained from the positions of the four Raman bands », 2v2 and 
(votv2—v2), 3v2—ve. The spacing 2v2—%, gives (e#+4W?)}, and (3ve—v2) — 
(vy; tve—ve) gives (®+8W2)), if anharmonicity is neglected. In these expres- 
sions, 6 is the separation of the unperturbed levels and W is the matrix element 
of the cubic term aj.2 in the potential energy which gives rise to the perturba- 
tion. One obtains 6 = 126.0 cm and |W) = 35.4 cm“. 

If anharmonicity is taken into account, one finds 


[22° — v,°|? = 4W7*+6, 
(3) [3vo! — (vi tvs!) )? = 8W2+(6+4x22—X12)’, 
[4vo? — (vi +2ve*)]? = 12W?+ (6+8x22— 2x12)", 
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where 6 = 2w»°+4xo2—w,"—x11, as before. On substituting the observed fre- 
quency differences, 144.13, 160.88, and 175.54 cm—', respectively, one obtains 


6 = 124.80 cm™! and |W| = 36.05 cm“, 


very close to the values obtained when anharmonicity is neglected. The small- 
ness of the effect of anharmonicity on 6 and |W’! is due to the smallness of the 
value of 4x%.2—x,}2 = —0.355 cm~!. The shift in the vibrational levels due to 
Fermi resonance is 9.67 cm~! up for (02°O) and 9.67 cm~! down for (10°0). 
According to the above definition, 6 is a positive quantity, but only the 
magnitude of W can be determined from Eq. (3) and not its absolute value. 
The latter can be found in the following way: The rotation vibrational coupling 
constant a» is given by the formula (Dennison 1940) 
B.’, 4Bew: 8B Te" arm 


ty 2 ee 
) 


2 3 
w. (w3 —we ww. he 


On the right-hand side, the last term (containing a2) is the predominant one, 
and can be given in terms of W (see Herzberg 1945, p. 218). From these equa- 
tions and the experimental value of ag = —0.00023 (Allen, Plyler, and Blaine 
1956), it follows that W = —31.8cm~ in good agreement with the value deter- 
mined from the Raman data alone, and with a negative sign. 

The separation of the unperturbed levels, given by 6, is almost 90% of the 
observed separation (Av) and shows that the perturbation of the two levels is 
not strong. The mixing of the two levels affects their eigenfunctions and, of 
course, affects the rotational constants of the interacting levels. For example, 
the observed rotational constant Boo for the 1,0,0 level in terms of the un- 
perturbed constants Byoo and Bozo is given by the equation 


D0 9p? 
Byoo = Biot bBorw, 


where the coefficient a? = (Av+6)/2Av and b? = (Av—6)/2Av. In the present 
case a = 0.93 and 6? = 0.07. Therefore if one uses the experimental value of 
Bioo reported in Section D(1), one obtains Byoo = 0.10892 cm-! (or a, = 
0.00018 cm~') for the unperturbed value of the rotational constant of the 
1,0,0 level. 

The vibrational constants and anharmonic coefficients (except for x13) can 
now be determined without neglecting the effect of Fermi resonance. The sums 
of the wavenumbers of a resonance multiplet (diad, triad... ) are given by the 
same formula as without Fermi resonance. For example, »;+2v2° = w)°+x11 
+2w.°+4x22 = 1460.09 cm. Similarly, (vi+v2!)+3y.! = 2254.96 cm, 2p," 
+ (vi +2y2°) +429 = 4380.86cm™!, (vy) + 22? — 2v2?) + (4v2? — 2ve") = 1462.82 cm. 
Other information used was 6 = 124.80 cm~ and 4xo2.—x12 = —0.355 cm7! 
evaluated from Eq. (3), and vz = 396.8 cm—!. The constants w3°, X33, and x23 
remain as in (2) since they were calculated from the bands v2, v3, 3v3, and 
vo+3v3—v2, which are not affected by Fermi resonance. Also since the available 
data are not sufficient to determine x3, the assumption is made that the diad 
(vitys, 2v2+3), of which only v1; +»; is observed, can be described by the same 
values of 6 and W found for the diad (»;, 2v2). This assumption leads to a 
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shift of 9.67 cm~! to lower frequency for the v;+ v3 band from its unperturbed 
position. The constants (in cm~!) obtained from the preceding analysis are 


as follows: 


wy," = 668 . 69, 1 = —1 054, X12 >= 0.86, 
wo® = 395.97, Xe = 0.126, Xo3 = —6.67, 
( 4) w3? = 1537 . 64, +33 = —5. 135, X13 = —4 9s, 


124.80, Z22 = 0.67;, 
W = —36.05. 
The frequencies (in cm~') for infinitesimal amplitudes, w;, obtained from 
the above w,° and x;, (Herzberg 1945, p. 206) are 


& 
ll 


w, = 671.36, we = 398.62, w3; = 1551.92. 


Apart from errors in measurement or assignment, any uncertainty in the 
constants given here derives largely from the unknown accuracy of the value 
$xo2—X,2 = —0.355 cm~'. This quantity was determined from the difference 
of two large numbers (see Eq. (3)) and hence may contain a considerable error. 

Unfortunately, all of the available infrared and Raman vibrational data 
have been used in the determination of the constants (4) so that a check of their 
accuracy is not possible at present. For this purpose and for the improvement 
of the present set of constants, it would be desirable to investigate the v2 
fundamental and other overtone and combination infrared bands at high 
resolution. 

E. VIBRATIONAL SPECTRA OF C®0, AND C0, 

The vibrational Raman spectrum of CO: is characterized by two intense 
and sharp Q branches, of almost equal intensity, which form the diad 7, 2v2. 
Adjacent to each of these lines are the much weaker Q branches of the related 
diad v;+v2—v2, 3ve—v2, in which the lower level is 0,1',0. The Q branches of 
these diads have been photographed for both C’Oz and C¥QO:, and are shown 
in Fig. 1(c). The »; Q branch of CO, is sharp and symmetrical; in fact this 
Raman line clearly shows the hyperfine components of the Hg 4358 exciting 
line. On the other hand, the 2v2 Q branch is slightly broader and is degraded 
towards the v; band. For CQsz, the »; and 2v2 Q branches are degraded towards 
each other. Also the v; band of CQO, is about three times the intensity of the 
2v. band, while in CO, the corresponding intensity ratio is approximately 
2: 1. The wavenumber shifts of these lines are given in Table IV. For CO, the 
average measurements of six plates are shown, the deviation from the average 
being at most +0.06 cm~'. For C¥QO., the average measurement of three plates 
is shown, except for the line »;+»2—v2, which was only measured on one plate. 


TABLE IV 
WAVENUMBER SHIFTS (UNITS CM~!) OF THE Q BRANCHES OBSERVED 
IN THE VIBRATIONAL SPECTRA OF C#Q2 AND C®O. 


3vo! — ve! 2v.° v9 019 +p2! — vo! 


C®O, 1265.17 1285.51 1388.15 1409.44 
C80, 1248 .23 1266.03 1369.90 1388 . 47 
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The separation of the Q branches in each of the two diads was used to deter- 
mine 6 and W, by the method described in Section D(2). The values (in cm~') 
obtained are: 

for C®O.: 6 = 15.8, W = —50.7; 
for C¥O.: 6 = 46.4, W = —47.1. 


These lead to shifts of 43.4 cm~! up for »; and down for 2y2° in COs, and to 
shifts of 28.7 cm! in CO ,. The values of 6 and W given here are necessarily 
approximate since they are based on Q branch measurements only and since 
higher order terms discussed by Taylor, Benedict, and Strong (1952) and 
Courtoy (1957) have been omitted. Nevertheless these values show that W for 
both isotopic molecules is nearly equal, as expected. Also, the resonance is 
not as strong in CQ, as it is in C"’Oe, as shown by the fact that the separation 
of the unperturbed levels in CO, is three times the separation in C?O.. This, 
of course, explains the much larger intensity ratio of »,/2v2 in C®QO, evident in 
the spectrum shown in Fig. 1(c). 

Several plates showing the weak rotational O and S branches of the C?O, 
diad »;, 2v2. were obtained and a photograph is shown in Fig. 1(d). The average 
wavenumber shifts measured on two plates are given in Table V. Unfortunately 
even an exposure of 30 hours (with a cylindrical lens in front of the photo- 
graphic plate) produced only a weak spectrum; only rotational lines up to 
J = 40 were observed, and these were difficult to measure. The analysis of 
the bands was carried out as outlined in Section D(1). 


TABLE V 
WAVENUMBER SHIFTS® (UNITS CM!) OF THE LINES IN THE v; AND 2po 
BANDS OF CQ, 


2v. Band 


v,; Band 


J S(J) OV) S(J) OJ) 

0 

2 1393.58 

4 1396.56 1294.02 

6 1399.85 1379.28 1297.16 1276.65 

8 1403 .00 1376.42 1300.25 1272.94 
10 1406.10 1373.30 1303.42 1270.55 
12 1409.42 1369.88 1306.50 1267 .46 
14 1412.32 1366. 24 1309.63 1264.43 
16 1415.46 1363.96 1312.84 1261.27 
18 1418.55 1360.77 1315.88 1258.16 
20 1421.69 1357.72 1319.04 1254.99 
22 1424.78 1354.60 1322.20 1251.72 
24 1427.91 1351.71 1325.30 1248.83 
26 1431.10 1348.26 1328.48 1245.72 
28 1434.19 1345.28 1331.16 1242.67 
30 1437.16 1342.14 1334.91 1239.57 
32 1440.31 1339.04 1337.93 1236.45 
34 1448 .33 1335.92 1341.10 1233.14 
36 1446.60 1332.83 1344.23 1230.12 
38 1449.71 1329.78 1348 . 26 1227 .01 
40 1326.49 1350.82 
42 —_— 1353.77 
44 1320. 392 


* Average values of measurements of two plates. 





b 
F 
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The constants (in cm~') obtained are: 


Vi00 = 1388.15, Byoo— Booo = — 0.00003, 
Vo2°9 = 1285.40, Bo2%9— Booo = 0.00019. 


All of the values are in good agreement with those reported by Rossman, 
France, Rao, and Nielsen (1956) and by Courtoy (1957). However, the smaller, 
yet significant effects, such as D’ # D’’, observed by these authors were not 
evident in the present analysis. 

The smallness of the values of B’— B” immediately explains the sharpness of 
the Q branches of the »; and 2v2 bands; the larger value for the 2v2 band is in 
agreement with the observed broader line, and its + sign agrees with the 
observed degradation towards longer wavelengths. Furthermore, the small 
values of B’— B” confirm the identical values of the Q branch measurement 
and band origin of the »; band and the very close agreement of the corres- 
ponding values for the 2v2. band. Since the resonance in CO, is strong, the 
values of B’— B” given here do not even represent approximately the coupling 
constants a; and ag; for a discussion of the effect of resonance on the rotational 
constants of CO, the reader is referred to the detailed account by Courtoy 
(1957). 

F. CONCLUSIONS 

Although the rotational structure of some vibrational bands in CS» and 
CO, has been photographed with a 21 ft. grating, the spectra are rather weak. 
It is an unfortunate fact that the low intrinsic intensity of rotation—vibrational 
Raman spectra still presents formidable difficulties and one must resort to the 
use of high gas pressures and long exposure times. Yet in spite of these diffi- 
culties, the investigation of vibrational Raman spectra at high resolution 
should be pursued, since much useful information can be obtained in this way. 
Two examples of the type of information obtainable are shown by the present 
study. When rotational structure is resolved, the band origins and rotational 
constants can be determined; when no rotational structure is observed, but 
only sharp Q branches, the vibrational constants can be evaluated with 
relatively good accuracy. 

One of the most interesting features of the present investigation has been the 
appearance of many “‘‘hot”’ bands in the Raman spectrum of CS,. As shown in 
Section D(2), these have led to the determination of the effect of Fermi reso- 
nance on the v; and 2v2 bands, and to the evaluation of a set of vibrational 
frequencies and anharmonic constants. Clearly, this situation is not peculiar 
to the spectrum of CS». ‘‘Hot”’ bands may occur in the spectrum of any mole- 
cule having a low-lying vibrational frequency (<400 cm~'), and similar studies 
of their Raman spectra at high resolution may prove useful in determining 
vibrational and anharmonic constants. 

It is well known that Fermi resonance is of frequent occurrence in the spectra 
of polyatomic molecules, and the present results on CS2 once again demonstrate 
the importance of determining the effects of the resonance, even when the 
interaction is weak. For example, in CS» the unperturbed vibrational levels 
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are separated by 125 cm~!. Yet each level is shifted by 9.7 cm~ through their 
mutual interaction. As a consequence of this shift, there is a significant effect 
on the vibrational frequencies and anharmonic constants obtained from an 
analysis of the vibrational spectrum which includes resonance as compared 
with the values obtained neglecting resonance. Thus some care is necessary in 
evaluating the accuracy or usefulness of vibrational constants (and in some 
cases, of rotational constants) of molecules in which Fermi resonance is known 
or suspected to occur but has not been adequately treated. 
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A SEARCH FOR a-PARTICLES FROM THE DECAY OF Bi?” ! 


E. P. Hincks anp C. H. MILLAR 


With an appendix by G. C. HANNA 


ABSTRACT 


A search for a-particles from Bi?° has been carried out in Ilford C-2 bismuth- 
loaded nuclear emulsions maintained at 5° C. for 3 years before development. 
Contrary to the findings of Faraggi and Berthelot (1951), Riezler and Porschen 
(1952), and Porschen and Riezler (1956), no evidence was obtained of a-particles 
in the 3 Mev. region which could be unambiguously attributed to bismuth, the 
experimental lower limit for the half-life thus determined being 71/2 > 2X10" 
years. A detailed study of the energy available for the a-disintegration of Bi? by 
the methed of reaction cycles gives a value of 2.93+0.11 Mev. For this disinte- 
gration energy the half-life of Bi?®® predicted on the basis of a-decay theory and 
systematics is between 5X10!® and 5X10” years. 


INTRODUCTION 


Bi?®*, which occurs in nature, has long been known to be unstable against 
a-particle emission by about 3 Mev., and several attempts have been made to 
observe this activity. Jenkner and Broda (1949), using bismuth-impregnated 
emulsions, found no evidence for the a-activity of bismuth and concluded that 
the half-life is greater than 3X10" years. Faraggi and Berthelot (1951), using 
bismuth-loaded emulsions, observed some a-particle tracks which they attri- 
buted to bismuth and concluded that Bi?®* emits 3.15 Mev. a-particles with a 
half-life of 2.710!" years. The present authors carried out a similar experi- 
ment the following year and reported (Hincks and Millar 1952) that they found 
no a-particle tracks that could be unambiguously attributed to bismuth. 
Shortly thereafter Riezler and Porschen (1952) published the results of a study 
of bismuth-impregnated emulsions in which they attributed seven a-particle 
tracks to Bi?®® and deduced an a-particle energy of 2.9 Mev. and a half-life 
of 2X10!" years, in rough agreement with the work of Faraggi and Berthelot. 
In a later comprehensive paper Porschen and Riezler (1956) reported five 
additional a-particle tracks presumed due to Bi?’ and revised their estimate 
of the a-particle energy to 3.0+0.2 Mev. 

The present authors delayed formal publication of their results pending a 
proposed confirmatory experiment with counters; however, since this experi- 
ment has not been done and since no more conclusive results have appeared 
in the literature, the authors feel that the publication of their original data, 
re-evaluated in the light of the latest available information, will constitute a 
useful contribution to the question of the natural radioactivity of bismuth. 


EXPERIMENTAL PROCEDURE 


Type C-2 bismuth-loaded nuclear emulsions which were manufactured by 
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Ilford, Ltd. on February 10, 1949, were kept in a refrigerator at 5° C. until 
March, 1952, at which time they were examined and it was found that the 
emulsion had become detached from the glass backing. It was possible, how- 
ever, to develop and fix the loose emulsion and, by means of a final wash in a 
weak gelatin solution, to reattach the developed emulsion to a glass backing 
with less than 1% distortion in the lateral dimensions. For comparison pur- 
poses, an unloaded Ilford C-2 emulsion, 50 microns thick, manufactured in 
November, 1949, and stored with the bismuth-loaded emulsion, was developed 
in the normal manner. An area of 92 mm.? of the bismuth-loaded emulsion was 
searched for single a-particle tracks; this corresponds to a volume of 0.014 
cm.*, since the emulsion was actually 150 microns in thickness instead of the 
nominal 100 microns. The emulsion was studied with a Cooke, Troughton, and 
Simms nuclear research microscope, the initial scan being carried out at a 
nominal magnification of 405 (45X6X1.5) and the track measurements at a 
magnification of 1425 (95X10X1.5). In the background plate an area of 
68 mm.” (~0.0034 cm.*) was searched for single a-particle tracks. 


AUXILIARY EXPERIMENTS 
It was realized that there were several factors that could lead to the mis- 
interpretation of any tracks observed, so a series of test experiments and 
checks was carried out to eliminate such errors as far as possible. These 
auxiliary experiments will be discussed first, since the interpretation of the 
observations in the main experiment is dependent on them. 

The fact that bismuth-loaded emulsion differs significantly in composition 
from normal emulsion gives rise to the possibility that the usual range-energy 
relationship for a-particles may not be valid. The emulsion compositions as 
supplied by Ilford Limited are shown in Table I. 


TABLE I 
ATOMIC COMPOSITION OF ILFORD C-—2 EMULSION 


Atomic composition (g./cm.*) 
Emulsion . ae een ee 





type Ag Br I c H O N S Na Bi Total 
Normal 1.85 1.34 0.052 0.27 0.056 0.27 0.067 0.01 —_ - 3.915 
Bi-loaded 1.39 1.01 0.039 0.33 0.047 0.48 0.062 0.002 0.06 0.27 3.640 











Yagoda (1949, p. 83) suggests a method of range correction based on the 
Bragg and Kleeman (1905) observation that for a-particles atomic stopping 
power varies directly as the square root of the atomic weight. The application 
of Yagoda’s formula would lead to a correction factor of 1.0165 to convert 
ranges in bismuth-loaded emulsion to ranges in standard emulsion. It seems 
doubtful, however, that this empirical ‘‘law”’ is valid to a sufficient degree to 
justify using it to determine a correction factor to better than a few per cent. 
In addition Rotblat (1950) notes that in the low-energy region a-particle 
ranges in emulsion under ‘‘atmospheric”’ conditions are about 3°% greater than 
in “‘vacuum-dry” emulsion. Since the standard range—energy relationships are 
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given for ‘‘vacuum-dry”’ emulsions (Lattes, Fowler, and Ciier 1947; Rotblat 
1951) some adjustment must be made to take into account the fact that our 
“‘refrigerator-desiccated’’ emulsion presumably had a moisture content some- 
where between “‘vacuum-dry”’ and ‘‘atmospheric”’ during its sensitive lifetime. 

Since there would appear to be some uncertainty in making each of these 
corrections, we have chosen to calibrate the emulsion using a-particles of 
known energy, and obtain an experimental total correction factor in this way. 

In the first calibration experiment both faces of a stripped, bismuth-loaded 
emulsion were exposed to a-particles from a thin U*** source. The emulsion 
was then developed, reaffixed to a glass backing, and the lengths of a group 
of surface tracks measured. The mean range of 44 of these tracks was 19.7+0.3 
microns, which agrees with the value given by the standard range-energy 
relationship for U*** a-particles (mean energy 4.808 Mev.) and indicates that 
ranges in the bismuth-loaded emulsion are within +1.5% of those in standard 
dry emulsion. 

The second check experiment gave information on several aspects of the 
main experiment and thus merits detailed consideration. A general survey of 
the bismuth-loaded plate showed that the great majority of radioactive 
“stars’’ from successive a-disintegrations of naturally radioactive contamina- 
tion in the emulsion originated from Ra**®. These four-pronged stars were 
easily distinguishable from the rare five-pronged stars from thorium, and a 
considerable number of them showed the characteristic separation by a few 
microns of the Ra®”® a-particle track from the remaining three prongs of the 
star—a well-known phenomenon caused by the diffusion of Em*” during its 
3.8-day half-life. These separated Ra**® tracks provided an unambiguous set of 
monoenergetic a-particle tracks deposited throughout the emulsion which were 
presumably typical of the tracks from any long-lived a-particle-emitting iso- 
tope contained in the emulsion. Careful measurements of these Ra”® a-particle 
tracks were carried out, and in addition the tracks were divided into three 
categories—‘‘dense’’, ‘‘medium’”’, and “‘faint’’—depending on their observed 
grain density. The results of this survey are shown graphically in the histograms 
of Fig. 1 and numerically in Table II. 

First it may be noted that the mean range of the dense tracks is 19.63+0.2 
microns, which is 1.3+1.0% greater than the 19.4; micron range predicted 
for the 4.766 Mev. (mean energy) a-particles of Ra®® from the standard 
range-energy relationship for dry emulsion. This is in good agreement with the 
result obtained with U*** a-particles. We are thus confident that when the 
ranges of low-energy a-particle tracks in our bismuth-loaded emulsion are 
multiplied by 0.99 to correct for atomic composition and moisture content, 
a-particle energies may be reliably determined from the relationships of Lattes, 
Fowler, and Ciier (1947) or of Rotblat (1951). 

The range straggling of the tracks may be conveniently measured by com- 
puting the r.m.s. deviations about the mean of the observed values. These 
r.m.s. deviations are shown in Table II, and it will be noted that the values 
increase by a factor of two from the dense tracks of small dip angle to the 
faint tracks of all angles. The increase in the r.m.s. deviation caused by adding 
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the steep dense tracks—arbitrarily defined as those with an observed dip of 
more than 3 microns—to the flat dense tracks is indicative of the large effect 
of the dip measurement inaccuracy in these special plates, whose shrinkage 
factor is as high as 4.6. The increase of r.m.s. deviation with decreasing track 
density is indicative of the increasing probability that some of the tracks will 
be shortened by fading of the end grains. 
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Fic. 1. Histograms of various classes of Ra®* @-particles measured in the bismuth-loaded 
emulsion in which the search for a-particles from Bi?®? was carried out. Predicted mean range 
for standard dry emulsion is shown dotted. Solid curve in (A) is smoothed distribution fitted 
to histogram. 


TABLE II 
MEAN RANGES AND R.M.S. DEVIATIONS OBSERVED FOR VARIOUS CATEGORIES OF 
Ra226 @-PARTICLE TRACKS IN BISMUTH-LOADED C-2 EMULSION 








r.m.s. 
Mean range = deviation 

Type of track No. observed (microns) (%) 

Dense (excluding steep tracks) 18 19.6.+0.2 1.9 

Dense (all tracks) 25 19.6;+0.2 6.1 

Medium (all tracks) 11 18.6,;+0.5 3.2 

Faint (all tracks) 43 18.1,;+0.3 9.9 


All tracks observed 79 18.7 +0 


|} NS 
~ 


These r.m.s. deviations are considerably greater than the value 2.9% which 
Rotblat (1950) obtains for a-particles of this energy in normal emulsion; 
however, our value of 4.9% for dense ‘“‘flat’’ tracks is comparable to the value 








HINCKS AND MILLAR: SEARCH FOR a-ACTIVITY IN Bi20 235 


of 4.7% (our calculation) obtained by Faraggi and Berthelot (1951) for the 
tracks from Th?” in a bismuth-loaded emulsion similar to ours. This greater 
apparent range straggling in bismuth-loaded emulsion is presumably due to 
the loading, which causes the silver bromide grains to be less closely packed, 
as well as to the large shrinkage factor. 

It will be noted from Fig. 1 that the addition of the ‘“‘medium”’ and ‘‘faint’”’ 
tracks to the ‘‘dense’’ ones not only lowers the value of the mean track length 
and increases the r.m.s. deviation, but also produces a skew distribution with a 
tail on the low-energy side. This is seen clearly in Fig. 1A, where a smooth curve 
is fitted to the experimental points and thus represents the empirical distri- 
bution function for a-particles of a single energy emitted at a constant rate 
throughout the sensitive lifetime of the emulsion. 

Table I] shows also that the mean range for the dense tracks does not change 
with the addition of the steep track measurements to those of the flatter tracks, 
indicating that our shrinkage corrections are accurate despite the very large 
value of shrinkage factor. 

A third experiment was carried out to check on the fading of tracks produced 
uniformly throughout the emulsion during its sensitive lifetime. This consisted 
of searching a large area (2.05 cm.*) of the emulsion for the tracks of y-mesons 
which came to rest therein and comparing the number found with the number 
predicted. Identification of these tracks was of necessity based on a visual 
estimate that the scattering near the end of the track was characteristic of 
u-mesons, it being impossible because of fading to make use of the grain 
density. This made it difficult to distinguish short or steep meson tracks from 
proton tracks; however, each track that entered and stopped in the emulsion 
was classified visually as either a meson or a proton, and the “‘mesons’’ thus 
selected were further classified as either “certain” or ‘probable’. 

From the emulsion composition, the meson intensity, and the time of expo- 
sure it was calculated that between 50 and 100 y-mesons should have stopped 
in the area scanned. A more precise estimate of this number was not possible 
because of the lack of accurate knowledge of atomic stopping powers at low 
velocities corresponding to u-mesons of ranges ~100 microns and less, and also 
because the absorbing material surrounding the emulsion during its storage 
period was neither controlled nor recorded. 

We observed 16 “‘certain’’ and 37 ‘‘probable’’ meson tracks terminating in 
this area and, by assuming that roughly as many protons were counted as 
‘‘probable’’ mesons as mesons were misidentified as protons, we estimate 
that a total of about 50 true meson tracks were observed. Thus the product of 
the ‘‘fading factor’ and our scanning efficiency for this type of track is in the 
range 0.5 to 1 and on this basis we feel sure that for the more densely ionizing 
a-particles the value of this over-all efficiency factor is at least 0.5. 

Another qualitative observation on the fading of tracks in the emulsion was 
provided by the radium stars for which a search was conducted in connection 
with the range-energy calibration. Because of their characteristic appearance 
the efficiency for observing faded stars should be greater than for similarly 
faded single a-particle tracks. However, no stars were observed in which the 
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individual a-particle tracks were faded to a degree such that if the tracks had 
been separate they would have been missed by the observer. This result adds 
confirmation to our previous conclusion that in our experiment fading of the 
latent image did not cause a serious loss in observing efficiency. 

Finally, to check that the emulsion really contained bismuth, samples were 
ashed, and the residue, when analyzed qualitatively, was found to contain a 
large percentage of bismuth. We were thus satisfied that there had been no 
mix-up of emulsions nor large loss of the bismuth component during either 
manufacture or storage. 

On the basis of these check experiments, we feel confident that for a- 
particle tracks in bismuth-loaded emulsions we know the range-energy rela- 
tionship with considerable accuracy, as well as the range distribution caused 
by straggling, fading, and measuring error. We also know that the over-all 
loss of tracks due to both fading of the latent image and observer inefficiency 


was less than 50%. 


MAIN EXPERIMENT RESULTS 


The histogram in Fig. 2A shows the corrected ranges of the 212 single a- 
particle tracks (excluding Ra*® tracks) found in the bismuth-loaded plate. 
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Fic. 2. Histograms of short single a-particle tracks found in (A) bismuth-loaded plate 
and (B) unloaded “background” plate. Solid curves show fitted ‘‘background”’ curve (see text). 
Dotted curve in (A) shows size and position predicted for peak due to Bi? if its a-particle 
energy were 3 Mev. and its half-life 210! years. Track lengths have been corrected for 
emulsion composition and moisture content. The a-particle energy scale is from Rotblat (1951). 
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Standard deviations are shown for the higher intensity points. The range- 
energy scale shown is from Rotblat (1951). 

Fig. 2B shows the range distribution of the 114 single a-particle tracks found 
in the background plate. There would appear to be no a priori reason why two 
batches of different types of emulsion should have the same amount of radio- 
active contamination and this is borne out by the observation that in the 
background plate the number of single a-particle tracks (per unit volume per 
unit time of ‘‘exposure’’) was 2.6 times that in the bismuth-loaded plate. It 
would seem that some method other than the subtraction of a “background” 
distribution of normalized intensity must be sought in order to determine 
whether the loaded plate contains an extra group of tracks which might be 
attributable to the loading. Our method of analysis is based on the following 
considerations. 

Single a-particle tracks in nuclear emulsion will most probably arise from 
contaminants from one of the three naturally occurring radioactive series; 
however, short tracks could also arise from the a-decay of Sm"? and the 
(n, p) reaction from the capture of thermal neutrons in the nitrogen of the 
emulsion. The characteristics of these ‘‘contamination tracks’ are shown in 


Table III. 


TABLE III 
POSSIBLE NATURAL SOURCES OF SHORT, SINGLE ‘“‘CONTAMINATION TRACKS” 
IN NUCLEAR EMULSION 


Mean 





a-particle Track length 
Source Isotope energy (Mev.) (microns) 
(4n+2) series U238 4.18 16.1 
Um 4.75 19.3 
Th 4.67 18.9 
Po?!0 5.30 22.6 
(4n+3) series 235 4.41 17.4 
Pa! 4.95 20.5 
(4n) series Th? 4.00 15.2 
Samarium Sm'7 2.18 7.3 
Neutron capture in N'4(n, p)C' _ 6.6 
nitrogen (El Bedewi 1951) 





We know from the abundance ratios and half-lives of U** and U* that in 
natural uranium contamination the activities of the members of the (4%+3) 
series will be only 4.5% that of the members of the (4n+2) series, assuming 
radioactive equilibrium exists within each series. We may also infer from the 
relative scarcity of five-pronged thorium stars as compared with the four- 
pronged radium stars that uranium contamination was a much stronger source 
of background tracks than thorium contamination in the emulsion studied. 
Thus, to a first approximation, we may attribute all single a-particles from 
contamination by elements of naturally occurring radioactive series to the four 
members of the (42+2) series listed in Table III. These a-particles have mean 
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ranges of from 16.1 to 22.6 microns and their tracks should thus be quite well 
separated from any bismuth a-particle tracks, which are expected to have a 
range of 10-11 microns, and from tracks due to either samarium a-particles 
or N'4(n, p)C" reactions, each of which has a range of about 7 microns. 

Considering in detail the members of the (4n+2) series that have long-lived 
daughter products so that their a-particles may appear singly in nuclear emul- 
sion, it is evident that if radioactive equilibrium has been established, a- 
particles from U**, U4, and Th? will appear in equal numbers unless there 
has been some chance separation of the few parts per hundred thousand of 
Th®° that exists in equilibrium with its uranium parents. Po*!° a-particle 
tracks will not necessarily be equally abundant, since in this case the source 
is probably the Em*? normally present in the atmosphere rather than the 
uranium contained in the emulsion. 

On these premises a semiempirical curve to compare with the observed distri- 
bution of a-particle tracks in our background plate was obtained in the follow- 
ing manner. Three smoothed distributions as obtained in Fig. 1A of equal 
amplitude and of mean ranges corresponding to the energies of a-particles 
from U**, U**, and Th*8° were added to a fourth such distribution of mean 
range corresponding to Po?!® a-particles. The relative amplitude of the Po?!® 
distribution was adjusted so that the mean range of the combined distribution 
corresponded to the mean range for all tracks observed between 10 and 28 
microns, the range limits of the semiempirical distribution. This combined 
distribution, normalized to the total number of tracks observed within its 
range limits, is shown as a solid curve in Fig. 2B, and it will be seen that it 
provides a statistically good fit to the observed points. 

In applying this same technique to the tracks observed in the bismuth- 
loaded emulsion it was noted that the mean range of tracks between 10 and 
28 microns was 19.3 microns as compared with 19.9 microns for the back- 
ground plate. These differ by nearly twice the standard error of the difference 
of the two means (0.35 microns) indicating either a lower relative abundance 
of Po*® or, possibly, the presence of a new group of low-energy particles. A 
semiempirical distribution fitted on the basis of the first of these possibilities 
does, indeed, give a statistically good fit to the observed distribution. However, 
if a new group of short tracks is present, this fitting process is not justified. To 
check this second possibility we have shown in Fig. 2A the curve obtained by 
renormalizing the semiempirical distribution obtained for the background 
plate. It will be seen that on the lower side of the peak the curve provides a 
statistically adequate fit to the observed points and there is little indication 
of any excess of tracks which might be attributed to bismuth a-particles. On 
the other hand, the experimental points in the region of the Po”!® peak fall 
below the curve, which is some evidence that the shorter mean range of tracks 
in the bismuth-loaded plate is due to a relatively lower abundance of Po?!® 
rather than to a new group of short-range tracks. 

In neither observed distribution is there definite evidence for the 4.00 Mev. 
a-particles from Th*” nor the 2.18 Mev. a-particles from Sm'7. Porschen and 
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Riezler (1956) have suggested that our technique was insensitive to low-energy 
a-particles and cite our non-observation of tracks in the 2-3.2 Mev. region as 
proof. In their own work they find tracks in this region in their background 
plates which they attribute to Sm" and to protons from N'4(, p)C" reactions. 

From the known flux of cosmic-ray neutrons we estimate that of the order 
of 10 N'(n, p)C" reactions would have occurred in the area searched in our 
bismuth-loaded emulsion during its sensitive life. However, for the proton 
tracks from this reaction, we expect a greater loss from fading than for the denser 
a-particle tracks, and the extreme shortness of these tracks (6.6 microns) 
presumably leads to a greater loss through non-observation than for the longer- 
a-particle tracks studied in the check experiments. Thus we might well have 
for the short and lightly-ionizing proton tracks an over-all observation effi- 
ciency less than 0.5, which means that we might have expected to observe 
less than five tracks from neutron capture in the nitrogen. The fact that only 
one track was observed in this range region in the bismuth-loaded emulsion 
(and two others in the background plate) is not incompatible with the above 
predictions, but it does not give firm proof that our technique was, indeed, 
sensitive to low-energy a-particles. However, on the basis of our efficiency for 
observing lightly-ionizing y-meson tracks, of the consistency of our results in 
measuring a-particles of 4-5 Mev., and of our previous experience in studying 
extremely short a-particle tracks (e.g. Millar and Cameron 1953), we feel 
confident that our technique was not anomalously insensitive to a-particle 
tracks in the 3 Mev. region. 

The dotted curve in Fig. 2A shows the predicted distribution for the 42 
a-particles that we would expect to have observed if the half-life of Bi?°® were 
2X10"? years and the a-particle energy 3.0 Mev. as suggested by Porschen 
and Riezler (1956), and if our observation efficiency were only 0.5. Bearing in 
mind that the lower limit of our detection efficiency has been used, and that 
no cosmic-ray background has been subtracted, it would appear from our 
observations that the upper limit of the specific activity of Bi?®® is about an 
order of magnitude lower than that represented by the dotted curve—that is 
the specific activity is less than 3X 10~5 disintegrations per second per gram. 
The corresponding half-life limit is 74,2(Bi?®) > 2X10'*® years. The experi- 
mental limit of the half-life will vary, but not very rapidly, with the a-particle 
energy chosen. We estimate that it should be approximately a factor of two 
smaller at 3.5 Mev. and a factor of two greater at 2.5 Mev. As will be shown 
in the next section, reaction-cycle calculations indicate that the energy avail- 
able for a-particle emission from Bi? is probably < 3.0 Mev. 


THEORETICAL CONSIDERATIONS 
Disintegration Energy 
Neither semiempirical mass formulae nor the methods of a-decay systematics 
(e.g. Perlman, Ghiorso, and Seaborg 1950) are adequate to predict accurately 
the a-disintegration energy of Bi?®*, which has a magic number of neutrons 
(126) and one proton more than magic (83). Mass spectroscopic data are also 
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not accurate enough for our purpose, but a fairly good estimate of the available 
energy can be made on the basis of reaction cycles. Fig. 3 is a chart of the 
nuclides involved in pertinent closed cycles for which experimental data are 
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Fic. 3. Diagram showing pertinent decay and reaction energies used in the determination 
of the energy available for the a-disintegration of Bi?°® by the method of reaction cycles. 
Arrows indicate direction of decay or of reaction on which values have been based. Numbers 
in parentheses are experimental limiting values. Numbers in square brackets are derived from 
reaction cycles. Error limits are shown beneath the appropriate figure. 


On this diagram are shown the natural a- and 8-decays and the artificially 
induced reactions (n, y), (y, 2), (p, 2”), etc., with the direction of the decay 
or reaction indicated. The experimental total energies—disintegration energies, 
neutron binding energies, or Q-values—are given, the lower figures showing the 
uncertainties in these values. Where experimental data give only a limiting 
value for the energy difference, this is shown in parentheses, while values 
derived indirectly from reaction cycles are shown in square brackets. The 
sources of these experimental energies and of their uncertainties are given in 
Appendix I. 

The energy available for the a-decay of Bi?, E,.(Bi?®), may be obtained 
without reference to any derived energy values from the cycle 


(1) Bi29 — Po268 — Ph2 — T1204 — T1298, 
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The energies involved (see Appendix I) are: 


Q(Bi2°*(p, 2n)) = —9.65+0.09 Mev. 


E,(Po*") = 5.21+0.01 Mev. 

— E3(T) = —0.76+0.01 Mev. 

E, (TI?) = 7.35+0.12 Mev. 

M,-M, = 0.78+0.00 Mev. 

therefore E,(Bi?°) = 2.93+0.15 Mev. 


There are two other possible cycles that are independent of (1). These are 


(2) Bi2° 2 Pb?99 ad Pb?2°8 as Pb?97 = Pb?°6 a Pb29 os T12% 
and 
(3) Bi2°9 <. Bi?! ce T12°6 = Pb2°6 a Pb2 id TI2%, 


but both of these unfortunately involve a neutron binding energy for which 
only a lower limit is directly measured, that is E,(Pb?°%) in cycle (2) and 
E,,(Bi?!”) in cycle (3). 

In the case of cycle (2), it is possible to trace an alternative route from 
Pb?°? to Pb? via the 4n+1 and 4n+2 radioactive series for which all energy 
values are known from experiment, thus avoiding the uncertainty introduced 
by including £,(Pb*°’) in the calculation. Referring to Fig. 3 this alternative 
route is Pb2°9 ate Pa? = Np? as U237 = U238 na Pb?2°S, 


The energies involved in cycle (2) via this alternative route are: 


— E,(Bi?) = —0.63+0.01 Mev. 

— E6a.3g(Pa?**) —40.38+0.04 Mev. 

— E,(Np”3’) —4.96+0.01 Mev. 

— E,(U*87) = —0.51+0.01 Mev. 

E,(U™**) 5.82+0.10 Mev. 

Ega.6g(U?4*) 51.64+0.06 Mev. 

— E,, (Pb?) —8.09+0.09 Mev. 

Ex-cap-(Pb*®) = 0.05+0.03 Mev. 

therefore E,(Bi?®) = 2.94+0.16 Mev. 


ll 


ll 


It is not possible to derive a value for £,(Bi!°) which is independent of 
cycles (1) and (2), so that cycle (3) cannot be used to provide an independent 
value of E,(Bi?®’). 

Combining the values from cycles (1) and (2) we obtain a best estimate for 
E,(Bi?) of 2.93+0.11 Mev. 

It is of interest to note that the best derived value of E,(Bi?!°) is a few 
hundred kilovolts above the value obtained from (mn, y) and (d, p) reactions in 
agreement with arguments advanced by Pryce (1952) and Harvey (1953). 

Similarly the best derived value of E,(Pb*°’) is greater than that obtained 
from a (d, p) reaction which is consistent with the argument that, since it is 
necessary to add a neutron with four or five units of angular momentum to 
Pb?°§ to produce Pb? in the ground state, it is unlikely that the ground state 
transition would have been observed in a (d, p) reaction. 
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Lifetime 

Assuming a value for the Bi?®’ a-disintegration energy enables its half-life 
to be estimated by analogy with the decay characteristics of neighboring 
bismuth isotopes. Such an estimate has been carried out by G. C. Hanna (see 
Appendix II) and his prediction of the half-life as a function of the a-dis- 
integration energy is shown by the solid curve in Fig. 4. In particular, we note 
that the estimated half-life for E, = 2.93 Mev. is 1.5107! years, although 
the uncertainties in the theory would allow it to be as low as 3 X 10°° or as high 
as 6X 10” years. 
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Fic. 4. Diagram showing results of various experiments on the a-decay of Bi?®* together 


with the prediction from a-decay systematics of the half-life as a function of disintegration 
energy (see Appendix IT). 


CONCLUSION 

Ina search of 0.014 cm.’ of bismuth-loaded emulsion that had been ‘‘exposed”’ 
for 3.1 years we have found no evidence for a group of a-particle tracks that 
could be attributed to Bi?®’. The disintegration energy obtained from reaction 
cycles using the best available experimental data is 2.93+0.11 Mev., corres- 
ponding to an a-particle energy of 2.87+0.11 Mev. We have used this result 
as a basis on which to obtain from our experimental data the best estimate of 
the lower limit set on the Bi?® half-life. This limit is 2 10!§ years and is well 
below (see Fig. +) the best estimate based on a-decay systematics of 5X 10° 
to 5X10” years for the range of disintegration energy quoted above. 

Our result is in disagreement with previous results of Faraggi and Berthelot 
(1951) and Porschen and Riezler (1956), who found an apparent a-activity 
due to bismuth with a half-life of ~2X 10" vears and an a-particle energy in 
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the region of 3 Mev. This half-life is inconsistent with the predictions of a-decay 
systematics unless the a-disintegration energy is at least 3.2 Mev., a value that 
is rather improbable on the basis of reaction-cycle analysis. 
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APPENDIX I. ENERGY VALUES FOR REACTION-CYCLE CALCULATIONS 
All energy values and their uncertainties in this Appendix are quoted to 
the nearest 10 kev. 
(a) a-Disintegration Energies 
E,(Po**) = 5.21+0.01 Mev. (Perlman and Rasmussen 1956) 


E, (Bi?) = 5.03+0.02 Mev. (Levy and Perlman 1954) 
E.(Np***7) = 4.96+0.01 Mev. (Magnusson et al. 1955) 


(b) B-Disintegration Energies 


E;(TP") = 0.76+0.01 Mev. (King 1954) 
E;(TP"*) 1.51+0.01 Mev. (King 1954) 
E;(Pb?"’) 0.63+0.01 Mev. (King 1954) 
E;(Th**) = 1.23+0.01 Mev. (King 1954) 
E;(U™7) = 0.51+0.01 Mev. (King 1954) 
Ey cap (Pb?) = 0.05+0.03 Mev. (see below) 


The value for E;- c¢ap.(Pb?°*) is based on the observation of Huizenga and 
Wing (1956) that Pb®®® decays by Z-capture with a half-life of 5X10" years 
and on the conclusion by Herber et a/. (1956) that decay by A-capture must 
have a half-life > 10'° years to account for the non-observation of Tl AK X-rays. 
We have assumed that the decay energy lies between the AK- and L-capture 
thresholds and have chosen an energy and uncertainty value to encompass 
this range. 


(c) (p, 2n) Reaction Energy 


Bi?’*(p, 2n)Po?°> Q = 9.65+0.09 Mev. (Andre et al. 1956) 


Andre et al. give +80 kev. as their cyclotron energy-calibration uncertainty. 
We have added +50 kev. in quadrature to account for a possible departure 
from the linear theory upon which their extrapolation is based. 


(d) Neutron Binding Energies 

Several neutron binding energies are based on the threshold for the (y, 7) 
reaction, and many of these measurements make use of the Cu®(y, 2)Cu® 
threshold energy for calibration purposes. The four most recent and nominally 
most precise measurements of this value are 
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E, (Cu®) 10.73+0.03 Mev. (Robinson et al. 1956) 
= 10.78+0.05 Mev. (Bendel et al. 1956) 
= 10.54+0.04 Mev. (de Souza Santos et al. 1955) 


10.61+0.05 Mev. (Birnbaum 1954). 


ll 


| 


ll 


Since it is evident that the differences between these four values are quite 
inconsistent with the errors indicated, we can only assume that each one is 
equally likely to be affected by a considerably larger systematic error and 
take the simple mean, 10.67 Mev., as the best value. We have arbitrarily 
assigned an uncertainty of 100 kev. to this value so that it is not in disagree- 
ment with any one of the four values quoted, that is 


E,(Cu®) = 10.67+0.10 Mev. 


This 100 kev. uncertainty has been added in quadrature to the quoted un- 
certainty of any (vy, ”) threshold energy referred to the Cu®(y, 7) threshold. 


E,(TE*) = 7.35-£0.12 Mev. 


This value is obtained from the three following determinations of the (y, 7) 
threshold: 

(1) Parsons and Collie (1950) give an “‘absolute’’ value of 7.3+0.25 Mev. 

(2) Sher et al. (1951) give a value of 7.55+0.20 Mev., but this is based on 
an energy calibration that gives the Cu®(y, 2) threshold as 10.85 Mev. Thus 
their value should be adjusted to 7.42+0.22 Mev. 

(3) Hanson et al. (1949) give a value of 7.48+0.15 Mev. based on a Cu® 
(y, n) threshold of 10.9 Mev. The adjusted value is therefore 7.32+0.18 Mev. 

The weighted mean value of these three corrected determinations is 7.35+ 
0.12 Mev. 


E,(TE°*) > 6.20+0.03 Mev. 
= |6.53+0.10 Mev.| 


The first value is based on the (n, y) reaction studied by Bartholomew and 
Kinsey (1953) and the (d, p) reaction studied by Harvey (1951). However, 
neither study proves conclusively that this represents a ground-state reaction, 
and indeed the indication of a 6.54 Mev. neutron binding energy ascribed by 
both to TP? on an intensity basis could possibly be properly ascribed to TH, 
or could mask a TI?°* binding energy in this region. Because of this uncertainty 
we believe that the best value is the one obtained from the reaction cycle 
T1295 — Pb26 — Pb? — T]?%, The value of £,,(TI?°*) thus derived is 6.53+0.10 
Mev. 


E,,(Pb*°*) = 8.09+0.09 Mev. (Harvey 1953; Palevsky and Hanson 1950) 


‘ 


Harvey’s value of 8.09+0.10 Mev. has been determined from the (d, ¢) 
reaction, using the Pb?°%(d, ¢)Pb?°? reaction for energy-calibration purposes. 
The latter is accurately known (see below). Palevsky and Hanson obtain a 
value of 8.25+0.10 Mev. from the (y, 2) threshold, but this value is based on a 
value of 10.9 Mev. for the Cu®(y, 2) threshold. Adjusting to 10.67 Mev. for 
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this latter value, we obtain E,(Pb*°*) = 8.08+0.14 Mev. Combining these two 
values and their errors, we obtain the value 8.09+0.09 Mev. 


E,,(Pb*°’) = 6.73+0.01 Mev. (Kinsey et al. 1951) 


The accurate value of the highest-energy y-ray of neutron capture is used 
here since confirmation that this is the transition to the ground state has been 
found in the (d, p) reaction (Harvey 1951), the (d, ft) reaction (Harvey 1953), 
and the (y, 2) reaction (Parsons and Collie 1950; Sher et a/. 1951; Palevsky 
and Hanson 1950). (The values in the last two references need to be adjusted 
to the proper Cu®(y, ) threshold.) 


E,,(Pb***) = 7.38+0.01 Mev. (Kinsey et al. 1951) 
The same arguments and references as for E,(Pb?”) apply here. 


E,,(Pb*?"’) > (3.87+0.05 Mev.) (Harvey 1951) 
= [4.14+0.11 Mev.] 


The value shown in parentheses is obtained from the highest-energy proton 
group observed in the (d, p) reaction and therefore can set only a lower limit 
on the neutron binding energy. This is confirmed by the decay cycle Pb?’ — 
Th? — Th?3 — Pa*3 — Pb?°°, which gives a value of >3.76+0.21 Mev., again 
only a lower limit since the only experimental value of /,(Th**) has been 
obtained from a (d, p) reaction. 

Absolute values of £,(Pb?°*) may be obtained from two different cycles 
which are independent except for the Pb?°* — Pb?°’ — Pb?°S chain whose energy 
separations are accurately known. 

The cycle Pb?°® — Pa — Np? — U?37 — U238 — Ph? — Pb?” — Pb?°8 gives a 
value E, (Pb?) = 4.14+0.12 Mev. 

The cycle Pb?2° = Bi2°9 a Po0?8 = Pb? aa TI? TH. Pb2% Pb2°6 bi Pb?°7 saad 
Pb2°8 gives a value /,,(Pb?°’) = 4.15+0.18 Mev. 

The weighted mean value of these two independent determinations is 
4.14+0.11 Mev. 


E,, (Bi?) > 4.17+0.02 Mev. 
= [4.44+0.11 Mev.] 


Here again both the (7, y) reaction (Kinsey et a/. 1950) and the (d, p) reaction 
(Harvey 1951; Wall 1953) give only a minimum value (>4.17+0.02 Mev.) for 
the neutron binding energy, while theoretical considerations indicate that the 
true value is probably somewhat higher (Pryce 1952; Harvey 1953). Here again 
two independent estimates may be made of £,(Bi*!®) but these involve cycles 
similar to those used for estimating £,,(Pb?°*). We therefore obtain a best value 
for E,(Bi##°) by using the derived value for E,(Pb*°’) in the cycle Bi*?® - 
Pb?°9 — Pb?°8 — Pb?0?7 — Pb?°6 — T12% — Bi#!°. The value thus obtained is E,(Bi?!®) = 
4.44+0.11 Mev. 


E,(Th***) > 4.88+0.20 Mev. (Harvey 1951) 
= (5.26+0.12 Mev.] 





246 CANADIAN JOURNAL OF PHYSICS. VOL. 36, 1958 


arvey’s experimental value is based on a (d, p) reaction an erefore sets 
H J I tal val based 1 t d theref et 
only a lower limit. 

The derived value consistent with other experimentally measured values 
may be obtained from the cycle Th? — Pb?°’ — Pb?99 — Pa? —- Np?8* and is 
5.26+0.12 Mev. 


E, (U*) = 5.82+0.10 Mev. 


From the (y, 2) reaction Parsons and Collie (1950) obtained an absolute 
value of 5.8+0.15 Mev. for E,(U*8*). From the same reaction Huizenga et al. 
(1951) obtained 5.97+0.10 Mev., but this was based on a Cu®(y, m) threshold 
value of 10.9 Mev. Correcting this latter value, we obtain 5.84+0.14 Mev. 





The weighted mean of these two values is 5.82+0.10 Mev. 


(e) Natural Radioactive Series Decay Energies 


(i) ‘‘4n’”’ series 


E,(Th??) = 4.08+0.01 Mev. 
E3(Ra®8) = 0.03+0.02 Mev. 
E3(Ac”’) = 2.24+0.02 Mev. 
E,(Th2’) = 5.52+0.00 Mev. 
E,(Ra™4) = 5.78+0.00 Mev. 
E,(Em”°) = 6.40+0.00 Mev. 
E,(Po?!*) = 6.90+0.00 Mev. 
E3(Pb?”) = 0.58+0.01 Mev. 
E,3(Bi2”) = 2.25+0.02 Mev. 
E,(Po?”) = 8.95+0.00 Mev. 








Total = 42.73+0.04 Mev. 


Th? (6a, 43 )Pb?°* 


E = 42.73+0.04 Mev. 


(Harvey et al. 1957) 
(see below) 

(Kyles et a/. 1953) 
(Asaro et al. 1953) 
(Briggs 1954) 
(Briggs 1954) 
(Briggs 1954) 
(King 1954) 

(King 1954) 
(Briggs 1954) 


The decay scheme and energy of Ra** are not well known; however, the 
value chosen has sufficient range to encompass most measured values (see King 


(1954) for references). 


(ii) “‘“4a+1” series 


Ea(Pa**) = 0.56+0.02 Mev. 
E,(U?*) = 4.90+0.00 Mev. 
E.(Th”*) = 5.11+0.01 Mev. 
E,(Ra™5) = 0.36+0.03 Mev. 
E.(Ac”*) = 5.92+0.01 Mev. 


E,(Fr2") = 6.45+0.01 Mev. 
E,(At?!?) = 7.18+0.01 Mev. 
E3(Bi?") = 1.39+40.01 Mev. 
E,(Po"*) = 8.51+0.01 Mev. 





Total = 40.38+0.04 Mev. 


Pa?*8 (6a, 38)Pb?°° 


E = 40.38+0.04 Mev. 


(King 1954) 

(Gol’din et al. 1955) 
(Hagemann et al. 1950) 
(Perlman et al. 1955) 
(Stephens 1955) 
(Stephens 1955) 
(Stephens 1955) 

(King 1954) 

(Stephens 1955) 
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(iii) ‘‘42+2” series U*55(8a, 63)Pb?°* E = 51.64+0.06 Mev. 
E.(U*) = 4.27+0.01 Mev. (Harvey et al. 1957) 
E,(Th?*4, UX) 0.19+0.01 Mev. (Stoker et al. 1953)? 
E,(Pa*4, UX.) = 2.30+0.01 Mev. (Stoker et al. 1953)? 


E,(U**4) = 4.85+0.00 Mev. (Harvey et al. 1957) 
E.(Th®®) = 4.77+0.01 Mev. (Briggs 1954) 

E.(Ra™®) = 4.86+0.01 Mev. (Briggs 1954) 

E,(Em””) = 5.59+0.00 Mev. (Briggs 1954) 

E..(Po*"*) = 6.11+0.00 Mev. (Briggs 1954) 

E3(Pb?"*) = 0.99+0.04 Mev. (King 1954) 

E3(Bi*"*) = 3.26+0.03 Mev. (Daniel and Nierhaus 1956) 
E.(Po*!*) = 7.83+0.00 Mev. (Briggs 1954) 


E,(Pb?", RaD) = 0.06+0.00 Mev. (King 1954) 
E,(Bi?!®, RaE) = 1.16+0.01 Mev. (Plassmann and Langer 1954) 
E,.(Po?!®) = 5.40+0.01 Mev. (Perlman and Rasmussen 1956) 


Total = 51.64+0.06 Mev. 
APPENDIX II. THE EXPECTED a-HALF-LIFE OF Bi? 
By G. C. HANNA 


(a) General 

In this appendix the a-half-life of Bi?®® is estimated, for a range of values of 
energy release near 3 Mev., by using a relation between lifetime and energy 
that is derived from a study of the experimental data on different types of 
a-emitting nuclei. 

The experimentally observed regularities in a-decay can be presented in 
several ways; we have chosen here the one originated by Winslow (1955). 
This approach recognizes that the dependence of lifetime on energy is con- 
trolled by the barrier penetration, an effect that is calculable if spherical sym- 
metry can be assumed. Thus for a (spherical) a-emitter of known energy and 
lifetime, the amplitude of the a-particle wave can be evaluated at any radius, 
and with the further assumption of a sharp nuclear boundary of known radius, 
at the nuclear surface itself. In this way the experimental data are used to 
give information on the specifically nuclear features of a-decay, that is on the 
probability of forming an a-particle out of nuclear matter. 

It is reasonable to expect that formation probabilities, or 
amplitudes”, will be similar for a-emitters of similar type. In particular it 
should be possible to predict, with some assurance, the lifetime of an unknown 
species if there exists a series of similar nuclei all having similar values of 


“Surface wave 


surface wave amplitude. 
There are essentially two uncertainties. First, especially near a closed shell, 
the value of either the surface wave amplitude or the nuclear radius of the 


2Error values quoted are our estimate based on the Fermi-Kurie plots shown in the reference 
quoted. 
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unknown species may differ from that of the rest of the series. Secondly, it is 
not possible to distinguish clearly between the inside and outside of the 
nucleus.’ In the present context this difficulty reduces to that of choosing an 
appropriate ‘‘effective’’ value for the radius, R; fortunately the results are not 
very sensitive to R. 

(b) Calculation of the Surface Wave Amplitude 

The only a transitions of interest for the present purpose happen to be ones 
in which LZ, the angular momentum of the a-particle, is uniquely defined. This 
is fortunate, since only in such cases can the surface wave amplitudes be cal- 
culated without additional information on the relative importance of the 
different L waves. 

The parameter to be calculated is the dimensionless product R|¢(2) |*, where 
¢(R) is R times the radial wave function; this product will henceforward be 
abbreviated to S. It can be shown (Winslow 1955) that 

S = R|¢(R)|? = MGR) = XE G(R)Qx 
where d is the decay constant, v is the velocity of the a-particle relative to the 
recoil nucleus, and G,;(R) is the irregular Coulomb function. G(R) is that part 
of G,(R) not containing L, that is 


— a \ 1/4 ‘ 9 : ; 
G(R) = =) exp) ee (are cos 2x12 


f2(L+})* a2 val 
a oer 
A more complete expression for Q; has been used by Winslow and the results 
given graphically. In these expressions 
B, the barrier height = 2Ze?/R, Z being the atomic number of the daughter 
nucleus, 
x = E/B, E being the a-disintegration energy (particle 
plus recoil), 
and k = Mv/h = (2ME)'/h, M being the reduced mass of the 
a-particle = (MaMpr)/(Mat+Mp). 


Numerically for R = 9.3 fermis these expressions reduce to 


—log S = 21.033+log 71/2(sec.) +43 logl/E(1+4/A)]—log Q,—2 log G(R) 


and Q. = exp 


with 


i ] =f) 1.0946 Z 1/2 1/2 \1/2 
9 rG = —} y ———————7-;5 larc COS X —Xx —x 
2 log G(R) a +i 4/A)|7 [are cos x x ‘"(1—x) "J, 


Xx 





8’Winslow formally avoids this difficulty by considering the wave amplitude near the nucleus, 
that is, at a radius large enough for the use of pure Coulomb wave functions. However, the 
whole procedure is useful only if it achieves a largely complete separation of the barrier pene- 
tration term. In fact Winslow uses a value of 9.3 fermis for the radius, which is certainly too 
small for the calculation to be accurate, but which is a reasonable estimate of the “effective” 
nuclear surface. 
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where A is the mass number of the daughter, E is in Mev., and all logarithms 
are to the base 10. 

These expressions refer to a bare nucleus. The effect of electron screening is 
taken into account by using a value for E larger than the observed disintegra- 
tion energy by an amount £,. Perlman and Rasmussen (1956) give 


E, = 65.3Z,7" —80Z,?” electron volts, 


where Z, is the atomic number of the parent nucleus. 
For bismuth a-decay E, = 31 kev. 


(c) Observed Regularities: Prediction for Bi°° 

The even-even a-emitters with a neutron number NV > 128 all have similar 
values of S, and lifetimes can therefore be predicted rather accurately. In other 
nuclei a-decay is generally less probable, especially for even-even a-emitters 
with N < 126 and the bismuth isotopes. Table AI lists the values of S for 
various relevant categories of a-emitters. Only ground-to-ground-state transi- 
tions are considered, since it seems probable that in Bi*®’ decay, transitions to 
excited states will be of minor importance (see below). Allowance has been 
made for the effect of spin-change, that is, the ‘‘centrifugal barrier’ term Q;; 
the values of Z were obtained from shell theory, except for Bi?” and Bi?!’ decay, 
where the values adopted by Perlman and Rasmussen (1956) have been used.‘ 


TABLE AI 
a-DECAY SURFACE WAVE AMPLITUDES NEAR THE CLOSED SHELLS Z = 82, N = 126 


—_ logioS 
L me 
Individual 


Nucleon 
— _ configuration 


a-Emitter 
Average (proton: neutron) 


Category 


Z= 


Example 
84, Po!8 .40 

J > 128 Po?é .43 (6h)?, (6g)? 
Po?! 42 ; 
Po? .60 


Po?!0 .85 
Po?98 .69 2.66 (6h)?, (4p)? 
Po? .54 


Po#4!(AcC’) f a. 37 3. (6h)?, 6g 4p 


Bi23 f 3.07 
B22 : 3.49 6h 3s, (6g)? 
Bit! 5 3.24 


Bi? ? f 6h 3s, (4p)? 


Bi(RaE) 22 5.55 6h 3s, 6g 4p 
Bi2!°(long) 4 7.42 


The last column of Table AI gives the shell model assignments (see e.g. 
Klinkenberg 1952) of the protons and neutrons from which the a-particle has 
to be assembled. Following the ideas of Perlman and Rasmussen (1956), it will 


‘The measured spin of Bi®®’ is 9/2, and that of TP is 1/2, so that, independent of shell model 


considerations, which fix the parity, it is known that L = 4 or 5. 
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be seen that there is a correlation between the value of S and the overlap of the 
nucleon wave-functions. For example, in the case of the polonium isotopes the 
considerable reduction in S with N < 126 is reasonably ascribed to the poorer 
overlap of the wave functions of the 64 protons and 4 neutrons, compared 
with the 6g neutrons above the shell. On the other hand the 126-neutron effect 
is not expected to be so large for the bismuth isotopes because of the broken 
proton configuration (one low and one high orbital), and it will be assumed 
here that the value of S for Bi?°’ decay is the same as for the decay of the 
heavier bismuth isotopes. 

Regarding a-decay to the first excited state of TP? (at 205 kev., probably ad3/2 
state), the nucleon overlap might be somewhat larger (6h 4d (4p)*) and the 
spin change smaller (L = 3), but these effects should be much less important 
than that of the lower energy release which corresponds to a factor of about 
700 in barrier penetration. This mode of decay should therefore be quite 
unimportant. 

A complication arises from the possibility of a significant difference in nuclear 
radius between Bi*®’ decay and the ‘“‘comparison”’ decays. It will be noted that 
S is rather smaller for Po?” decay than for the heavier polonium isotopes 
(and these have somewhat lower values of S than the emanation and radium 
a-emitters). This behavior might possibly indicate a shrinkage in radius as 
the closed shell is approached. However, this effect appears to be less serious 
than the other uncertainties. 

With the assumption that logioS is —3.25, values of a-half-life have been 
calculated for a range of disintegration energies. The results are given in 
Table AII, and plotted graphically in Fig. 4. 


TABLE AIlI 


CALCULATED VALUES OF HALF-LIFE FOR Bi? a-DECAY 


a-Disintegration Half-life 

energy (Mev.) (years) 
20s 0X 10?! 
2.82 3.5X 107 
2.92 5X10! 
3.02 5X10" 
3.12 4.3 10'8 
3.22 om 1G" 


Note: The disintegration energies are the ‘“‘observed”’ 
values, i.e. before the 31 kev. electron screening cor- 
rection is added. 


It is difficult to give a reliable estimate of the uncertainty. The error due to 
poor choice of R is not expected to exceed a factor of 2 (corresponding to a 
change of 25% in Rk). The value assumed for S is probably less reliable. How- 
ever, Bi?°* decay would certainly be expected to be no more probable than 
Po?!” decay, for which S is roughly 2.5 times greater; on the other hand a finite 
126-neutron effect would increase the predicted lifetime, by a factor of 20 if 
it were as important as for the polonium isotopes. 

Rather conservatively then, the half-life is expected to be not less than 0.2 
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times and not greater than 40 times as large as the value given, for any energy, 


in Table AII. 
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LETTERS TO THE EDITOR 





Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words and, for any issue, should be submitted not later than six weeks 
previous to the first day of the month of issue. No proof will be sent to the authors. 


Polarization of Electrons in Forbidden 8-Decay* 


As a special case of a lengthy calculation of the angular distributions of the electrons and 
neutrinos and the polarization of the electrons emitted in the arbitrarily forbidden B-decay 
of an oriented source, an expression for the longitudinal polarization of the electrons from a 
randomly oriented source has been obtained. Because of the importance and the simplicity 
of the result, it is believed that immediate publication of this special case may be useful. 

The longitudinal polarization is expressed as the quotient of a ‘‘modified correction-factor”’ 
and the usual correction-factor. The modifications are certain changes in the combinations 
of coupling constants and in the functions L,, M,, etc. of Greuling (1942) and Pursey (1951). 
Let us make the definition Dyy = Re(CxCy*+Cx’Cy’*). X or Y stands for one of the five 
interactions SVTA P; Cx and Cx’ are coupling constants as defined by Lee and Yang (1956). 
The usual correction-factor depends on the coupling constants only through combinations 
of the type Dxyy. For the modified correction-factor we shall requireGxyy = Re(Cx’Cy*+CxCy’*) 
and I'yy = Im(Cy’Cy*+CxCy’*). The modifications in the coupling constant combinations 
are as follows: 


(i) If X and Y both belong to STP, replace Dxy by Gxy. 

(ii) If X and Y both belong to VA, replace Dyy by —Gyy. 

(iii) If X belongs to VA and Y to STP, replace Dxy by Ixy. 

(iv) If X belongs to STP and Y to VA, replace Dxy by —Txy. 
Because [xy = —TI yx, cases (iii) and (iv) are equivalent. It is of interest to note that with 
Pauli’s (1957) new 8-decay interaction Dxy is replaced by the real part of the invariant Kyy, 
Gxy and yy by the real and imaginary parts of the invariant Lyy. 

In addition to the coupling constant changes we must replace the combinations of electron 
radial wave-functions according to the scheme given in Table I. 


TABLE I 
TRANSFORMATION OF L,, M,, ETC. 


| 
| 


eS nn —(p? Fop?”)"g_y_1 fai sin Q 


M,— (p? Fop?”*?) Np asf oat sin Q 


Ny =e (2p? Fop?’*)-"( gy 418-1 —frsif—v-1) sin Q 


Py — (p?Fop?”)-1g_y-1fr41 cos Q 


Q— (p?Fop??**)—19741 f_yn1 cos 2 


R,y—> — (2p? Fop?”*!)-"(2y412-»—1 —fossf—y-1) cos Q 


y? = (v+1)?—a2Z?, 





Except for Q, the notation of Table I is that of Rose, Perry, and Dismuke (1953); 
0 = 6+—6-, 5* and 5~ being the phase-shifts for angular momentum j = v+} and for the 
values (7+3) and —(j+4) respectively of the quantum number x. 


*Issued as A.E.C.L. No. 533. 
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A useful first approximation is to treat the nuclear charge distribution as a point charge 
and to evaluate the wave-functions to the lowest order in pp. The expressions in the right- 
hand column of Table I correspond to this approximation. When these formulae are applied 
to allowed transitions, agreement with the results of Jackson, Treiman, and Wyld (1957) is 
obtained. If in addition the two-component neutrino theory (Lee and Yang 1957) is assumed, 
the longitudinal polarization is +v/c for interactions drawn from either STP or VA. If the 
ratio Cy’/Cx = +1 has the same value for each of STP and has the opposite value for each 
of VA, the polarization is again +v/c. By demonstrating that the approximate relations of 
the last column of Table I fail when one retains terms of the next higher order of pp in the 
radial wave-functions, one sees that the formula +v2/c for the polarization is not exact. 

The results given so far apply to negatron decay. For positron decay one must, in addition 
to changing the sign of Z in Table I, carry out a suitable transformation of the coupling 
constants. If Cx+ and Cx’* refer to positron emission, it is sufficient to take Cyt = Cxy* when 
X is V or T, Cx+ = —Cyx* when X is S, A, or P, and Cy’t/Cxt = —Cy’*/Cx* for all X. 
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The Viscosity of Liquid Helium between 0.78° K. and 1.1° K. 


The theory of the viscosity of liquid helium IT (Landau and Khalatnikov 1949) gives two 
expressions for the variation of the viscosity » with temperature T 


Io K.:< P< te R., n = A+BT-V2eA/kT, 
T < 0.8°K.,, n = A+(CT¥%e-4/kT 4. DTS)—, 


A, B, C, and D are constants which can only be estimated, and 4/k has the value 8.8° K. 
(Wiebes et al. 1957). The viscosity measurements above 1.1° K. made by Heikkila and Hollis 
Hallett (1955) using a rotating cylinder viscometer were consistent with the theory, but an 
adequate test can only be made from measurements at lower temperatures. 

Another viscometer of the type described by Heikkila and Hollis Hallett has been used to 
measure the viscosity of liquid helium down to 0.78° K. in a booster cryostat (Atkins et al. 
1955). The radius and length of the inner cylinder were 0.995 cm. and 3.03 cm. respectively, 
and the width of the annular gap between the cylinders was 0.60 mm. The suspension for the 
inner cylinder was immersed in the liquid helium and the deflections of the inner cylinder 
were determined with the aid of two specially shaped condenser plates, one fixed, and the other 
attached to the inner cylinder. The measured change in capacity of this condenser was pro- 
portional to the angular deflection of the inner cylinder, and the constant of proportionality 
was determined from viscosity determinations above 1.1° K. assuming that the viscosity 

values were known (Heikkila and Hollis Hallett). The viscometer was attached to the appendix 
of a metal dewar of such a size that approximately 300 cm.? of liquid were cooled to 0.78° K. 
by pumping. Temperatures were measured using a carbon thermometer, which was calibrated 
against a paramagnetic salt, and are believed to be accurate to 0.01° K. The velocity of the 
outer cylinder was less than the critical value for this apparatus. 

Interpretation of the results was somewhat complicated by the fact that the effective 
torsion constant of the fiber suddenly increased somewhere between 0.95° K. and 1.00° K. 
This effect was due in some way to the superconducting state of the aluminum alloy of which 
the viscometer was made, because the torsion constant, as measured by the period of free 
oscillation of the inner cylinder, returned to its original value when the superconductivity 
was destroyed by a magnetic field. 

The results are shown in Fig. 1 and include measurements taken with the apparatus in its 
normal state (i.e. ina magnetic field), and in its superconducting state, corrections being made 
for the observed change in torsion constant. The curve which has been drawn represents the 
theoretical variation of » with temperature using the constants A, B, and A/k given by Heik- 
kila and Hollis Hallett. (The known slight variation of B with temperature (Khalatnikov 
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1956) and use of the latest value of A/k (with corresponding changes in A and B) do not change 
the shape of the curve significantly. ) 
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The experimental results agree, at least qualitatively, with the predictions of the theory. 
Quantitatively, the value of 140+30 micropoises for the viscosity at 0.78° K. compares 
favorably with the value of 173 micropoises at 0.8° K. quoted by Khalatnikov. The measured 
value at 1.01° K. is 3747 micropoises. It should be pointed out that the theory does not give 
specific predictions for the range 0.8° K. < T < 1.0° K.; this region can only be covered by 
interpolation (i.e. the dashed curve in the figure). 

These results also agree with the indirect determinations made by Zinov'eva (1956) from 
measurements of the absorption of second sound down to 0.83° K. 
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